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It is well known that classical ECM techniques against an incoming active missile 
attack is effective only at distances where the jamming signal to missile return signal 
ratio is high. There is no guarantee that it will work well against all missiles. As the 
missiles get closer to the target, the effectiveness of classical ECM techniques is 
reduced. Hence, there is a need to design a short range ECM system that is effective at 
medium range and increases its countermeasure capabilities as the range shrinks. The 
US Army is considering a high power, mobile ECM system capable of delivering up to 
5 GW of peak power at 10 GHz. The narrow beamwidth of this system creates a very 
high peak power density of the order of 10-20 W/cm 2 at 10 Km distance. With this 
amount of concentrated energy, it is capable of damaging or degrading the highly 
sensitive sensors of the missile and causing the missile to lose track of its target. 

The purpose of this project is first, to evaluate the various possible array types 
which must meet a far field peak power density requirement subject to a constraint on 
maximum aperture E-field strength. At the same time, the sidelobe and backlobe levels 
should not cause radiation hazard to the operators of the system. After an array 
design is chosen, the whole vehicular structure is modeled and near field patterns are 
examined to determine the likely radiation hazard zone. 

This project is based on the usage of the Numerical Electromagnetic Code - Basic 
Scattering Code (XEC-BSC) from the Ohio State University. Four additional software 
programs are written to support the code. The evaluation shows that Uniform - 
Cosine tapering offered the lowest sidelobe and backlobe level in the near field and yet 
met the far field requirement and maximum aperture E-field constraint. With the 
vehicular structure included, the system has a safety margin of more than 27 dB below 
the recommended radiation hazard limit near ground level. 
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I. INTRODUCTION 



A. GENERAL BACKGROUND 

The classical electronic countermeasure (ECM) techniques against incoming 
active missiles attack is to employ noise and deceptive jammers to capture the missile 
sensor's range gate and angular gate and steer it toward some false target generated by 
passive ECM devices. Such an ECM technique works well when the distance between 
missile and jammer is large and the jamming signal to missile return signal ratio is high. 
On the other hand, modem missiles are designed to counter classical ECM techniques, 
such as the use of a monopulse tracking technique which makes angular gate stealing 
very difficult, and the home-on-jam technique to home onto the jammer when the 
jamming signal is very high. As the missile closes in on the jammer, the jamming 
signal to missile return signal ratio decreases and the jamming effectiveness drops. By 
the time the distance between the missile and the jammer is reduced to less than 
10 Km, there is even less chance that the jammer will achieve its mission. Hence, there 
is a need to design another form of ECM technique that is effective at the short range. 

One solution to the short range ECM problem is to use a high power narrow • 
beam device, concentrating very high energy towards the missile at ranges up to 
10 Km. As most sensors are designed to have high sensitivity for target acquisition at 
far distances, they usually do not have very wide dynamic range for handling very high 
input power density. If the power density produced by the ECM device at 10 Km 
exceeds the breakdown limit of the missile sensor, the sensor may be damaged or 
malfunction, causing the missile to lose tracking ability. 

B. SYSTEM SPECIFICATIONS 

The Flyswatter Array is a high power, vehicle mounted, ECM system capable of 
delivering up to 5 GW of output peak power within a very narrow pencil beam. This 
creates a very high power density of the order of 10-20 W/cm 2 at a distance of 10 Km. 
This amount of energy when radiated onto a seeker head of an incoming missile will 
damage or degrade the highly sensitive sensor of the missile, thereby causing it to lose 
track of its target. The specifications of such a system are listed below : 



Frequency 

Total output peak power 
Maximum aperture E-field 



10 GHz 
5 GW 
10.6 KV/cm 
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Reflector size 
Array size 

Array mounting plate size 
Array type 
Element orientation 



5.1m x 2.7m x 0.1m 
2m x 2.5m 
3m x 4.5m 

endfeed broadside array 
length in parallel with the 
long side of the plate 



Full gain elevation coverage 

Azimuth coverage 

Pulse width 

Pulse repetition rate 

Required power density at 10 Km 

Overall structure 



0° to 45° 



10 - 20 W/cm 2 
(see Fig. 1.1) 



180° 

100ns 

lOpps 



The above specifications were derived from a simplified model consisting of uniformly 
distributed elements. More detailed description of such a simplified model will be 
described in Chapter IV. 

C. PROJECT PURPOSE 

Given that we have such a high power device, one can ask "what is the radiation 
intensity in the vicinity of the system? Will it cause a radiation hazard to the operators 
of the vehicle?" The purpose of this project is to simulate the ECM system using a 
more realistic TE 10 mode slot array rather than uniformly distributed elements. The 
far field patterns are recorded to ensure that the power density requirement of 
10-20 W/cm 2 at 10 Km is met. Then the near field radiation patterns are plotted to 
determine the radiation hazard regions, if any. 



13 




Figure 1.1 Flyswatter Array Mounted on an Armored Vehicle. 



14 



II. RADIATION HAZARD EVALUATION PROCEDURE 



A. PROBLEM AREA IDENTIFICATION 

The project can be divided into two main areas. The first one is to determine the 
array configuration to meet system requirements and the second is to determine any 
radiation hazard zones. 

1. Array Configuration Determination 

The system specifications listed in Chapter I, section B do not define the array 
configuration. The only restrictions are on the array type, element orientation and 
overall array dimension. Other parameters that must be chosen are as follows: 

• the dimension of element, 

• the spacing between elements, 

• the types of field tapering and 

• the number of elements. 

By varying all these parameters, we hope to design a TE 10 mode slot array 
that will meet the system requirements; namely, deliver up to 5 GW of output power to 
achieve 10-20 W/cm 2 of power density at 10 Km distance without exceeding the 
maximum aperture E-field breakdown limit of 10.63 KV/cm. At the same time, the 
sidelobes and backlobes must be kept small so that the near field radiation level outside 
of the main beam is low enough not to cause any hazard to the operators of the 
system. 

2. Radiation Hazard Regions Determination 

After the array configuration has been defined and the system specifications 
are met, the near field patterns of the ECM system are then computed at various 
reflector tilting angles and array azimuthal pointing directions. The radiation hazard 
regions are then recorded. 

B. EVALUATION PROCEDURE 

The evaluation procedure used to accomplish the task is as follows : 

Step 1 : Compute the far field patterns of a simple plate mounted array with fewer 
elements and no reflector plate. The purpose is to determine the element 
spacing and element sizes so as to reduce unwanted grating lobes. 
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Step 2 : Compute the far field and near field patterns of the array with the reflector 
plate, with different tapering types. This will determine the type of tapering 
that will give high gain, narrow main beam beamwidth and low sidelobe 
and backlobe levels. 

Step 3 : Compute the far field patterns of the flyswatter array with full input power 
of 5 GW. By adjusting the number of elements and the type of tapering, 
one can determine if the required power density at 10 Km distance is 
achieved. At the same time the maximum aperture E-field must also be 
checked to ensure that it is within the limit specified. 

Step 4 : Compute the far field patterns of the array in the presence of the modeled 
vehicle and determine if the requirements are still met at various elevation 
and azimuthal pointing angles. Vary the number of elements and tapering 
type if necessary. 

Step 5 : Compute the near field patterns of the array in the presence of the modeled 
vehicle at various elevation and azimuthal pointing angles and determine 
the radiation hazard zones. 

C. RADIATION EXPOSURE LIMITS 

The radiation exposure limits of radio frequency electromagnetic spectrum for 
humans used in this report is based on the recommendations made by the National 
Council on Radiation Protection and Measurements (NCRP) and American National 
Standards Institute (ANSI). The NCRP report [Ref. 1] suggested two different limits 
on the far field power density, averaged over 6 min. and 30 min. respectively as shown 
in Fig. 2.1(a). The average-over-6-min. curve is meant for occupational exposure 
criteria which is similar to the limits given in ANSI C95. 1-1982 protection guide 
[Ref. 2]. The average-over- 30-min. curve is one-fifth that of the occupational criterion 
and is the criterion for the general population. These limits apply to both continuous 
wave and pulsed signals. No specific peak power limitation was given. 

The ANSI C95. 1-1982 protection guide suggested that the radiation intensity be 
measured at a distance of 5 cm or more from any object in the field and averaged over 
a 6 minute duration as shown in Fig. 2.1(b). Hence, it can be applied to both near 
field and far field radiation conditions. At 10 GHz, the power density limit is 
5 mW/cm 2 , the same as that given in NRCP report. Therefore, we will use this as the 
occupational exposure limit. However, the Flyswatter Array has a duty cycle of 10*^, 
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hence the limit on peak power density based on an average of 5 mW/cm* is equivalent 
to 5 KW.cnr or 77 dB \V,m“. 




1000 

Frequency, HHz 



(a) NCRP Proposed Limits. 



Frequency (/) 
range 


Equivalent 

power 

density* 


(Electric field) 1 


(Magnetic 

field) 1 


MHz 


mW/cm’ 


VVm* 


A 1 /®’ 


0.3-3 


100 


4 X 10‘ 


2.5 


3.0-30 


900// 1 


4 X 10 3 (900//*) 


0.025(900//*) 


30-300 


1.0 


4 X 10* 


0.025 


300-1500 


//300 


4 X 10 3 (//300) 


0.025 (//300) 


1500-100.000 


5.0 


2 X 10‘ 


0.125 



•From ANSI (1932). 

b measured 5 cm or greater from any object in the field and averaged for any 0.1 h (6 
min). 

* (Electric Field) 2 / 1 200ir or 12r(magnetic field) 1 , whichever is greater. 

(b) ANSI C95. 1-1982 Protection Guides. 



Figure 2.1 Radiation Exposure Limits. 
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III. SOFTWARE SUPPORT 



The simulation of the Flyswatter Array and vehicle requires a software code that 
is capable of modeling the w'hole system structure and computing the near field and far 
field radiation patterns at X-band. The Basic Scattering Code of the Numerical 
Electromagnetic Code (NEC-BSC) series is chosen for its capabilities in computing 
radiation patterns of complex structures using the Geometrical Theory of Diffraction 
(GTD) method. This code provides the required calculation and analysis power and 
outputs the data in tabular form. Additional routines were then written to support the 
code. A brief description of each of these codes is included in this chapter. The names 
enclosed in the parenthesis are the file names of the codes. The flow charts and listings 
of the four supporting software codes are attached in Appendix A. A detailed 
description of the Basic Scattering Code and its applications can be found in the 
references [Refs. 3.4]. 

A. BASIC SCATTERING CODE (GTD2P51) 

1. Capabilities. 

• Frequency Range : 

* UHF and above. 

• Modeling : 

* Infinite Ground Plane : Perfect conductor or dielectric filled plate. 

* Plates : Perfect conductor plates or thin dielectric slabs. 

* Cylinders : Elliptical shape with slanted end caps. 

* Antenna : 

Dipole or slot. 

Single element or array. 

Uniform, piecewise sinusoid or TE Q1 distribution. 

Electric or magnetic excitation. 

* Receiver : 

Dipole or slot. 

Single element or array. 

Uniform, piecewise sinusoid or TE 01 distribution. 

Electric or magnetic excitation. 
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• Computation : 

* Far field : E-field in either 9 or <p direction or versus frequency. 

* Near field : E-field, H-field and Power density computations with one of 
the following options : 

Spherical coordinates. 

Rectangular coordinates. 

Frequency sweeping. 

* Antenna Coupling : Magnitude, Phase and dB values of reaction between 
the antennas with one of the following options : 

Mutual impedance between the two antennas. 

Ratio of the power-out over power-in using a modified Frii's 
transmission formula for conjugate matched loads. 

Maximum gains between the antennas using the Linville method. 

2. Limitations. 

• Plates should have edges at least a wavelength long. 

• Thin dielectric slabs should be at least a wavelength away from any source. 

• Incident fields should not strike the slab too close to grazing. 

• Each antenna element should be at least a wavelength from all edges. 

• The code delivers up to a second order answer, ie., it can include fields for ray 
paths that interact with at most two scattering centers. Hence, the plots are 
accurate for the first 30 dB only. 

• Double diffracted fields are not included. 

In many cases, the one wavelength limits may be reduced to a quarter wavelength for 
engineering purposes. 

3. Modifications. 

The GTD2P51 code has been modified to suit our requirements for this 
project. The modifications are carried out largely in the main program and subroutine 
OUTPUT which are listed below : 

• Added a CT card to read in the title of input data set. This title will be printed 
out together with the data set and passed over to the plotting routine to be 
included as the heading of the plots. 
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• Modify the subroutine OUTPUT to print out the title, abbreviated codes and 
output data at specific formats for output unit No. 8. The abbreviated code 
consists of 5 characters : 

* Character 1 : F, N for Far field or Near field. 

* Character 2 : F, S for Fixed or Scanned frequency. 

* Character 3 : R, S for Rectangular or Spherical coordinates. 

* Character 4 : E, H, P for E-field, H-field or Power density. 

* Character 5 : 1, 2, 3 defines the variable parameters; 

X, Y or Z for rectangular coordinates, 

R, 6 or (p for spherical coordinates. 

• Maximum dimensions for plates, cylinders and sources have been increased as 
follows : 

* Maximum number of plates = 50. 

* Maximum number of cylinders = 100. 

* Maximum number of sources = 10100. 

B. INTERFACING CODE (GTDRWSEL) 

The output from the Basic Scattering Code consists of a string of abbreviated 
control codes followed by several columns of data. The interfacing code (GTD Read, 
Write, SELect) reads in these abbreviated codes and converts them into useful 
information. It then displays this information for the user to select the type of plots 
that he wants. These selections include, 

• The type of field, ie. E-field, H-field or power density in either X, Y, Z or R, 
0, <p directions or the total field. 

• The type of plots, ie. magnitude plot or phase plot, linear or dB scale for Y 
axis values, linear or log plot. 

After the selection, the GTDRWSEL program will then generate another set of 
control codes to be used by the plotting program. These codes and the labels for each 
plot are placed in a control file using output unit No. 3. The program also chooses the 
correct data sets and places them in the output data file using output unit No. 2. 
These two files will be used by the plotting program (GTDPLOT). 
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C. OUTPUT DATA PLOTTING CODE (GTDPLOT) 

This program is modified from an existing program called PLOTNF (Near Field 
Plot). Additional features were added to allow for all options in the interfacing code 
(GTDRWSEL). It reads in the control codes from the control file and obtains the 
data from output data file and plots them using DISSPLA subroutines. 

D. GEOMETRIES PLOTTING CODE (GTDGEOMP) 

This program reads in the same set of input data that is to be fed into the N’EC- 
BSC code. It then selects only the geometrical data and plots out its 3-D geometries 
using DISSPLA subroutines. It is used as an aid to check the geometrical data of the 
input data set. 

The GTDGEOMP code selects the geometrical data and plots the following 

items: 

• Infinite ground plane, 

• Plates, 

• Cylinders and 

• Sources. 

* 

It also reads in and executes the following commands that are used in GTD2P51 : 

• CT Print title. 

• NC Remove all cylinders. 

• NG Remove ground plane. 

• NP Remove all plates. 

• NS Remove all sources. 

• NX Start next set of data. 

• RT Rotate and translate coordinates. 

The user may choose to view the object from any viewing angle in order to give a 
better perspective of the structure. This code, however, does not have the capability to 
blank out hidden lines. 

E. SA CARD GENERATING CODE (GTDSA) 

The code SA stands for Source-Array. This is an interactive program where the 
user is asked to enter a series of data regarding the geometries of the array, the type of 
element used, the tapering type and the total input power. It then computes and 
generates a data file in the format specified by the NEC-BSC code for the SA card. 
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This program is capable of generating a slot antenna array with a rectangular 
grid arranged on the X [t Z t plane, as used in NEC-BSC code. All elements of the array 
are assumed to have the same dimensions and orientation. It can simulate either a 
TE 10 mode or Uniform distribution for each element. There are 5 different tapering 
types to choose from, for each of the directions (along X t or Z t axis). Namely, 
Uniform , Cosine , Cosine-square , Triangular and Binomial . One can also choose to 
add a pedestal height to the tapering. Note that even without any pedestal (ie. 
pedestal height equal zero), the outer most elements will still have some excitation 
because the algorithm always includes one more element at each side of the array, 
having the same spacing as the rest, to be the one with zero excitation. 

The distribution used in each element is assumed to be perfect so that the 
relationship between the E-field and the output power for TE 10 mode excitation is 



abE 2 




A 

where E is the maximum E-field across the width b of the element, a is the length of 
the element and rj is the intrinsic wave impedence. For an array of N elements with 
total output power P tota p the above equation becomes 

P t°tal " ^f n ti E n ’ 

A 

with E n being the maximum E-field of the n* element. Knowing the output power, 
element dimensions and the tapering types, the excitation for each element E n 2 can 
be computed. 

For a Uniform distribution, the relationship between the input power and the 
excitations is given below : 



P 



total 



ab 

n 




Note that given the same input power, element dimensions, number of elements and 
tapering types, the element excitation of each element for uniform distribution is only 
one half of that of the TE J0 mode. Hence, a Uniform distribution will have much 
lower aperture E-field than TE 10 mode. 



IV. ARRAY CONFIGURATION ANALYSIS AND RADIATION HAZARD 

EVALUATION 

The evaluation of radiation hazard is carried out according to the 5 steps spelled 
out in Chapter II under Evaluation Procedure. 

A. ELEMENT SPACING AND SIZE VARIATION 

In order not to complicate the evaluation, we assume that the elements are 
arranged on a rectangular grid and all elements are of the same rectangular size and 
orientation as shown in Fig. 4.1. The axes X and Z £ orientate the array. (The array 
has no fixed relationship with the global co-ordinates X, Y and Z until specified.) 
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Figure 4.1 Plate mounted rectangular grid array. 
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There are three main factors that impose restrictions on the element spacing and 
sizes. Namely, supporting the TE J0 mode only, suppressing grating lobes and sidelobes 
and reducing maximum aperture E-field. 

1. TE10 mode Propagation 

In order to support the TE 10 mode, the waveguide elements must have a 
length of at least half a wavelength (X/2). Similarly, for the TE Q1 mode to propagate, 
the width of the element must be at least half a wavelength, too. Hence, we have to 
restrict the width to less than half a wavelength. As the length increases, and exceeds 
one wavelength, the TE 20 mode will begin to propagate. Therefore, in order to support 
only the TE 10 mode, the element size must be restricted to the following values : 

1.0X > element length (a) £ 0.5X 
0.5X > element width (b) 

2. Grating Lobes and Sidelobes Suppression 

Grating lobes will appear whenever the spacing between elements is greater 
than or equal to one wavelength. This can be seen in Fig. B1 of Appendix B. These 
three plots show the far field patterns of a simple plate mounted array with fixed 
spacing of I. OX, but different element lengths of 0.5X, 0.75X and l.OX. There is no 
reflector plate used and the other parameters that are common to all three plots are as 
follows : 



Array mounting plate size 


4.5m x 3m 


Number of elements in X t direction 


10 


Number of elements in direction 


10 


Element width 


0.25X 


Amplitude Tapering in X f direction 


Uniform 


Amplitude Tapering in Z ( direction 


Uniform 


Output pow r er 


100 W 


Frequency 


10 GHz 


Distance 


10 Km 



The array is placed on the X-Y plane with the X t axis parallel to the Y-axis. 

These figures show that grating lobes always occur at 0 = 90° and 270° 
irrespective of the element length. The main difference between the three plots is the 
gain at 0 = 0°. It is highest for element length = 1.0X and reduces as the length 
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decreases. Hence, it may be desired to have the length of elements as large as possible 
to ensure higher gain, but not too large to create grating lobes. 

Fig. B2 shows the effects when the spacing between elements in the Z t 
direction is reduced. Fig. B2(a) and Fig. B2(b) are for element spacing of 0.75X while 
Fig. B2(c) shows the radiation pattern plot when the spacing is further reduced to 0.5X. 
Note that at 0.5X spacing, the antenna gain pattern at 0 = 90° and 270° is totally 
cancelled as one would expect. It is very important to ensure that the radiation level 
near the region where 0 = 90° and 270° be kept to its minimum because these 
directions could be in line with the operator. 

Similar deductions can be drawn for the width of the element and the spacing 
in X t direction. Fig. B3 shows 3 plots of fixed element length and spacing in both the 
X t and Z t directions, but with different element widths of 0.5X, 0.375X and 0.25X. The 
w'ave in the element is still propagating in the TE 10 mode and hence the E-field 
strength is uniform in the X t direction. This results in broader beamwidth and even at 
a spacing of 0.75?i in the X direction, the grating lobes appear to be quite significant. 
As the spacing in the X t direction is reduced to 0.5X, as shown in Fig. B4, the grating 
lobes disappear. Further reduction in the spacing causes the radiation level at 0 = 90° 
to increase as shown in Fig. B5. Hence, the spacing in the X direction should be fixed 
at 0.5X. 

3. Maximum Aperture E-field Radiation • 

The equation below shows the relationship between the input pow r er and the 
maximum E-field of the array, assuming a perfect TE 10 mode of propagation. 



abE' 



total = 



max 



4Z r 



^ ^ ratio^ n ) 



From this equation, it is obvious that by increasing the element size a and b, the 
maximum E-field E max will be reduced since E rat j 0 (n) is only a function of N and the 
tapering type. Hence, by making the elements as large as possible, without affecting 
the mode of propagation and the grating lobes, the maximum aperture E-field may be 
reduced. 

4. Effects of varying the element spacing and size 

From the above investigation, the effects due to the variation of the element 
spacing and size are summarized below'. 
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• Spacing between elements in the Z t direction should be less than one 
wavelength in order to reduce grating lobes. It may be desired to operate 
with spacing equal to 0.5X to reduce the radiation at regions near 
0 = 90° and 270°. 

• Spacing between elements in X t direction must be 0.5X. or less to reduce 
grating lobes. At 0.5X spacing, the radiation level at 0 ^ 90° is very low 
and is preferred. 

• Element size affects the gain of the main lobe as well as the maximum power 
flux at array output for a fixed number of elements. The larger the element 
size, the higher the gain and the lower the maximum power flux. 

From the above observations, one realizes that the spacing between elements 
in both Z t and X t directions are restricted to 1.0X and 0.5X respectively in order to 
support only the TE 10 mode of propagation and at the same time suppress the grating 
lobes. Hence, by operating with minimum gaps between elements in both directions so 
that the element size is as large as possible, one can achieve high array gain at low 
aperture E-field strength. Alternately, if radiation hazard is a more critical issue, then 
one would prefer to operate with element size of 0.5X x 0.25X. and with element spacing 
of 0.5X in both directions to ensure minimum direct radiation for 0 between 90° and 
270°. 

Before we conclude this section, let us consider the case where the array area 
is fixed and we want to fill the area with as many elements as possible so that the peak 
excitation is reduced. Fig. B6 to B9 show that a fixed array area of 7.5X x 3.75X is to 
be filled by elements of length of either 0.75X or 0.5X and width of either 0.375X or 
0.25X. The plots in Fig. B6 and B7 are for Uniform tapering and they show identical 
radiation pattern in both Eg and E^ cuts. Even the maximum excitations are all equal 
to 0.2442 x 10 4 V/m with 100W of input power. The plots in Fig. B8 and B9 are for 
Cosine tapering and again the radiation patterns are almost identical except for some 
variation in the sidelobe level. However, the maximum excitation varies with the plots, 
with 0.4349 x 10 4 V/m for 0.75X x 0.375X, 0.4462 x 10 4 V/m for both 0.75X x 0.25X 
and 0.5X x 0.375X. and 0.4578 x 10 4 V/m for 0.5X. x 0.25X. Hence, the peak excitation 
for Cosine tapering increases as the element size is reduced, while that of the Uniform 
tapering remains unchanged. 
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B. TAPERING TYPES VARIATION 

Five commonly used tapering types are selected for comparison. These are : 

• Uniform tapering 

• Cosine tapering 

• Cosine-Square tapering 

• Triangular tapering and 

• Binomial tapering 

Fig. BIO to B14 show the far field and near field patterns of the array with 
different taperings for each figure. The tapering in X t and Z t directions for each plot is 
the same to simplify the evaluation. The array is placed on the X-Y plane with the X t 
axis parallel to the Y-axis and the reflector plate is tilted at 45° to give maximum 
radiation pointing towards 0 = 90°. The other parameters chosen are as follows : • 



Fig 4.2 shows the geometry of the array. 

1. Far Field Pattern Comparison 

For far field radiation patterns, we are interested in the gain and beamwidth. 
Hence, the plots are taken with ±45° coverage about the boresight for cuts in both 
elevation and azimuth. The origin of the coordinate system is located at the center of 
the reflector plate and the plots are computed at a distance of 10 Km. Table 1 show r s 
the comparison of gain and beamwidth. 

From Table 1, we have the following relationship : 

Gain : Uniform > Cosine > Triangular > Cosine-Square > Binomial 
Beamwidth : Binomial < Cosine-Square < Triangular < Cosine < Uniform 



Reflector plate (single plate) 

Array Mounting plate 
Number of elements in X t direction 
Number of elements in Z t direction 
Element length 
Element width 

Element spacing in X t direction 
Element spacing in Z t direction 
Output power 
Frequency 



5.1m x 2.7m 
4.5m x 3m 



10 

10 

0.5X 
0.25X 
0.5X 
0.5X 
100 W 



10 GHz 
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Figure 4.2 Simplified Array Geometry. 

Note that the reduction in gain and the increase in bearmvidth for Binomial tapering as 
compared to the others are quite significant. Hence the Binomial tapering may not be 
suitable. 

2. Near Field Pattern Comparison 

For near field comparisons, we are interested in radiation levels in all 
directions. Hence, 0° to 360° coverage in both vertical and horizontal cuts are 
obtained. These patterns are computed at 5m distance from the center of the reflector 
plate. 

The plots show that with Uniform tapering, the sidelobe levels are very high 
and they occupy quite a wide area, hence, they are likely to cause a radiation hazard 
in a wider region. With Cosine tapering, the mainlobe beamwidth is slightly broadened 
but the sidelobe level is reduced by about 10 dB. The mainlobe beamwidth is even 
wider for Cosine-Square and Triangular tapering. However, the sidelobe level is also 
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TABLE 1 

COMPARISON OF GAIN AND BEAMWIDTH 


Taping 

Parameter 


Uniform 


Cosine 


Cosine 2 


Triangular 


Binomial 


Max. E-field 
in dB V/m 


-27.15 


-28.29 


-29.90 


-28.95 


-32.47 


BW in 

vertical cut 
(deg) 


10.25 


12.50 


15.50 


13.50 


21.00 


BW in 

horizontal 
cut (deg) 


10.00 


12.50 


15.00 


13.50 


20.00 



proportionally reduced. As for Binomial tapering, the mainlobe bearmvidth is the 
broadest and theoretically, it should offer the lowest sidelobe level among the 5 
tapering types. But due to the broad beamwidth and the limited reflector plate size, 
there is spillover at 0 near 40°, making the sidelobe level appear to be much higher 
than for Cosine-Square and Triangular tapering. However, as the number of elements 
increases, the beamwidth will be reduced and Binomial tapering should offer the lowest 
sidelobe level. 

Notice that in the vertical cut near field plots, there is a considerable amount 
of radiation at 0 near to 180° and 300° for Uniform, Cosine and Triangular tapering. 
These additional sidelobes also appear in the far field pattern but are not shown in our 
plots due to the limited angular coverage. The first lobe near 180° is part of the 
radiation pattern continuing from the main lobe. The sudden drop in radiation level 
prior to this lobe is due to the presence of the array mounting plate. By simply 
increasing the plate size, this lobe can be eliminated. The second lobe near 0 = 300° is 
the result of one of the sidelobes generated by the array itself. The presence of the 
reflector causes distortion on this sidelobe. For Cosine-Square and Binomial tapering, 
the sidelobe levels of the array are so low that no significant sidelobes appear at the 
near field patterns beyond 0 = 180°. 
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3. Preferred Tapering Type 

The far field pattern comparison shows that Uniform tapering gives the 
highest gain and narrowest mainlobe beamwidth. On the other hand, the near field 
pattern comparison suggests that Cosine-Square or Binomial tapering offer the lowest 
sidelobe level. A compromise between gain and sidelobe levels would make Cosine- 
Square tapering the preferred choice. 

C. OVERALL ARRAY CONFIGURATION ANALYSIS 

In this section, we will attempt to determine the number of elements in both the 
X t and Z t directions, element sizes, as well as the tapering type in order to meet the far 
field requirement of 10-20 W/cm 2 of power density at 10 Km, subject to the following 
constraints as spelled out in the system specifications. 

• Maximum total output power of 5GW. 

• Maximum aperture E-field of 10.6 KV/cm. 

• Maximum array area of 2.5m x 2.0m. 

1. Array Model Simplification 

The array area is given as 2.5m x 2.0m. In order not to exceed the maximum 
aperture E-field constraint, one would expect to have the whole array area filled with 
array elements. Since we are operating at 10 GHz, and the element spacing is 
restricted to less than one wavelength in Z t direction and half a wavelength in X t 
direction, the minimum number of elements needed will be N = 84 x 133 = 11,172 as 
shown in Fig. 4.3(a). With such a large number of elements, it is impossible to execute 
the GTD2P5 program. Hence, a simplified array model is needed. 

The GTD2P5 code can model the TE 10 mode or uniformly distributed 
elements regardless of the element length and width. If we assume that tapering along 
the X t direction is Uniform, and there is no gap between elements in X t direction, then 
we may use a simplified array that has the length of each element unchanged but the 
width extended to include all elements in X t direction. That is, each element's width is 
equal to 2m to cover the whole width of the array area as shown in Fig. 4.3(b). Such a 
simplified array model will produce exactly the same solution as one using the actual 
model. 

Alternatively, if we wish to vary the tapering type along X t direction, we may 
replace the elements along the Z t direction with a single element of equivalent length 
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Figure 4.3 Array Model Simplification. 
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while the width remains unchanged. This configuration is shown in Fig. 4.3(c). 
Because the GTD2P5 code can only simulate either Uniform or Cosine (TE 10 ) 
distribution within the element, we are restricted to using either Uniform or Cosine 
tapering along Z ( direction. Together with the use of either a uniform distribution or a 
TE 10 mode excitation for each element, we have the following four configurations. 

(a) Uniformly tapered TE 10 linear array simulated by a single uniformly 
distributed element. 

(b) Cosine tapered TE 10 linear array simulated by a single TE I0 element. 

(c) Cosine tapered uniformly distributed linear array simulated by a single TE 10 
element. 

(d) Uniformly tapered, uniformly distributed linear array simulated by a single 
uniformly distributed element. 

The aperture E-field distribution of these configurations are shown in Fig. 4.4. 

Except for the last configuration, the simplified model will not yield an exact 
solution. Fortunately, the near and far field radiation patterns for the simplified 
models are very similar to those of the actual arrays. The relationship of the gain and 
maximum aperture E-field between the actual and simplified models are derived in 
Appendix C. These derivations allow us to accurately predict the performance of the 
actual array from the data obtained from the simplified model. 

To illustrate the similarity between the actual array and the simplified model, 
5 diagrams consisting of far field and near field pattern plots are shown in Fig. B15 to 

B19. These plots are measured at 10 Km away from the center of the reflector. The 

distance of 10 Km ensures that the field generated by any element within the array area 
of 2.5m x 2.0m will arrive almost in phase, so that direct addition of the individual field 
patterns is possible. The input power is 25W, the mounting plate is 4.5m x 3m and the 
reflector plate is titled at 45°. The elements are arranged in a straight line with a 
constant width of 0.25X along Z t axis. The other parameters of each plot are listed 
below: 

• Model 1 : Actual array for configuration (a) (Fig. B15). 

* Tapering type = Uniform 

* Aperture field distribution = TE 10 

* No. of elements = 20 

* Element size = 0.5X x 0.25X. 
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Figure 4.4 Aperture E Field (E a ) for Both the Actual Array 
and the Simplified Model. 
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* Element spacing = 0.5X 

• Model 2 : Simplified model for configuration (a) (Fig. B16). 

* Aperture field distribution = Uniform 

* No. of elements = 1 

* Element size = 10.0X x 0.25X 

• Model 3 : Actual array for configuration (b) (Fig. B17). 

* Tapering type = Cosine 

* Aperture field distribution = TE 10 

* No. of elements = 20 

* Element size = 0.5X x 0.25X 

* Element spacing = 0.5X 

• Model 4 : Simplified model for configuration (b) & (c) (Fig. B18). 

* Aperture field distribution = TE 10 

* No. of elements = 1 

* Element size = 10.5X x 0.25X. 

• Model 5 : Actual array for configuration (c) (Fig. B19). 

* Tapering type = Cosine 

* Aperture field distribution = Uniform 

* No. of elements = 20 

* Element size = 0.5X x 0.25X 

* Element spacing = 0.5X 

Fig. B15 and B16 allow the comparison of radiation patterns between the 
actual array and the simplified model of configuration (a). Fig. B17, B18 and B19‘ 
show the radiation patterns of the configurations (b) and (c). Except for the absolute 
magnitude, both the near and far field radiation patterns of the simplified models 
appear to be very similar to those of the actual array. Table 2 lists the maximum 
aperture E-field, boresight far field strength and boresight near field strength of each 
model. It also shows computed ratio of the values of the actual array to those of the 
simplified model. The computed ratios of the maximum aperture E-field and the 
boresight far field strength turn out to be very close to the derived value. Even the 
ratios for the near field field strength also appear to be closely related to the derived 
values for far field. Hence, we can use the data obtained from the simplified model to 
predict the actual array performance accurately. 
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TABLE 2 

MAXIMUM APERTURE E-FIELD AND BORESIGHT FIELD STRENGTH COMPARISON 





Max. Aperture 
E Field 


Boresight Far 
Field Field 
Strength 


Boresight Near 
Field Field 
Strength 


Model 1 


4095.0 V/m 


0.0192 V/m 


23.30 V/m 


Model 2 


2047.0 V/m 


0.0151 V/m 


18.30 V/m 


Model 1 

Ratio= 

Model 2 


2.0005 


1.2715 


1.273 


Derived Result 
for Con- 
figuration (a) 


2.0000 


1.2732 




Model 3 


5635.0 V/m 


0.0176 V/m 


21.40 V/m 


Model 4 


3996.0 V/m 


0.0196 V/m 


23.80 V/m 


Model 3 

Ratio= 

Model 4 


1.4102 


0.8980 


0.899 


Derived Result 
for Con- 
figuration (b) 


1.4142 


0.8986 




Model 5 


2818.0 V/m 


0.0139 V/m 


16.80 V/m 


Model 5 

Ratio= 

Model 4 


0.7051 


0.7092 


0.706 


Derive result 
for Con- 
figuration (c) 


0.7071 


0.7058 





2. Uniform Field Realization 

The system specifications are derived based on the simplified model of the 
array. It assumes that the field distribution within each element is uniform (not TE l0 ) 
and there is no gap between adjacent elements. The area is then occupied by 96 
elements each of length 0.8687>w and width 2m with the tapering along the Z ( axis being 
Triangular (linear). This simplified model generates 10.47 KV/cm of maximum 
aperture E-field. The far field patterns generated by this array are attached in Fig. B20. 
and the peak E-field of this array is about 78.81 dB V/m or equivalent to 20.2 W/cm 2 
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of power density. The near field patterns are attached in Fig. B21. Fig. B21(a) shows 
the vertical cut radiation pattern measured at 3.5m. It shows very high radiation level 
in both forward lobes and backlobes. Fortunately, the azimuthal beamwidth is very 
narrow as shown in Fig. B21(b) and the hazardous region is confined to within ± 10° 
of the boresight, and is much less for the backlobe since the peaks are much lower. The 
array will be mounted on a vehicle 2m high, hence, it is not possible for the mainlobe 
or the backlobe to point directly at any operator. Fig. B21(c) and B21(d) shows the 
horizontal radiation plots of the array at two sidelobe peaks near the antenna 
mounting plate where the operators could be located. At the front (cp = 90°) and back 
((p = 270°) of the array, the radiation is still above the hazardous level and hence 
operators should be kept away from the front and back of the array. 

3. TE10 mode Realization 

Although a uniform distribution is easier to simulate and it also simplifies the 
computations, it is not realizable in a slot antenna. Instead, a perfect TE J0 mode is 
chosen to make the simulation more realistic, bearing in mind that the peak E-field in 
each element of the TE ]0 mode excitation is much higher than that generated by the 
uniform distribution. 

In order to minimize the aperture 'field strength, we will fill the 2.5m x 2.0m 
array area with elements close to each other. Using the knowledge we gained in 
Section A of this chapter, we select the element size of 0.8687k x 0.5k so as to reduce 
the aperture field strength w’hile at the same time keeping the sidelobes as low r as 
possible. Once the element size is chosen, the number of elements is then fixed at 
96 x 133. With this configuration, we can then use the two suggested simplified model 
methods to determine the type of tapering that will meet the requirements. 

After several trials with TE 10 mode excitations, we find out that by using 
5 GW of input pow'er, the far field criterion can be met by most tapering types. But at 
the same time, the aperture field strength constraint is exceeded. Hence, the input 
power must be lowered. It turns out that when the input power is lowered, only 3 
tapering types will satisfy both the far field criterion and near field constraint. These 
are Uniform - Uniform tapering, Cosine - Uniform tapering and Uniform - Cosine 
tapering. (With the tapering type along Z t direction positioned before the hyphen and 
the tapering type along X t direction placed after it.) 
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Table 3 lists the maximum aperture field strength and the far field boresight 
field strength computed at 10 Km for the 3 tapering types mentioned above. The far 
field field strength is also converted to power density for ease of comparison. Two sets 
of data are obtained for each tapering type. One with 5 GW of input power and the 
other with 1.88 GW of input pow r er to meet the requirements. As can be seen from the 
table, Uniform - Uniform tapering offers the highest far field power density and lowest 
maximum aperture field strength. The Cosine - Uniform and Uniform - Cosine 
taperings provide approximately the same performance. 



TABLE 3 

COMPARISON OF MAXIMUM APERTURE FIELD STRENGTH AND 
FAR FIELD POWER DENSITY FOR DIFFERENT TAPERING TYPE 



Tapering 
• Types 


Input 

Power 

(GW) 


Max. Aperture 
Field Strength 
(V/m) 


Far Field 
Boresight 
Field Strength 
at 10 Km (V/m) 


Far Field 
Boresight 
Power Density , 
at 10 Km (W/cm") 


Uniform- 


5.00 


0.123E7 


0.128E5 


43.49 


Uniform 


1.88 


0.753E6 


0.799E4 


16.95 


Cosine- 


5.00 


0.173E7 


0.118E5 


36.96 


Uniform 


1.88 


0.106E7 


0.723E4 


13.88 


Uniform- 


5.00 


0.173E7 


0.118E5 


36.96 


Cosine 


1.88 


0.106E7 


b . 724E4 


13.91 



Fig. B22, B23 and B24 show the near field radiation patterns of the array 
using Uniform - Uniform tapering, Cosine - Uniform tapering and Uniform - Cosine 
tapering respectively. These plots are measured at 3.5m aw r ay from the center of the 
reflector plate with 1.88 GW of input power. Note that the NEC-BSC code suggested 
the output is only accurate for the first 30 dB of the pattern, due to the absence of the 
higher order diffraction/reflection terms. For these 3 plots, however, we are only using 
2 flat plates, the antenna mounting plate and the reflector plate, with the reflector plate 
slanting at an angle of 45° from the antenna mounting plate. Furthermore, most of 
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the energy is concentrated towards the center of the plates and then reflected 
horizontally out of the plates, causing very little interaction within the structure. 
Therefore, we can safely say that the accuracy range limit of 30 dB does not apply to 
these plots. Unfortunately, it is not within the scope of this project to construct an 
actual model for direct measurement and comparison. We will bear in mind the 
accuracy problem and place additional safety margins on the radiation hazard level. 

For Uniform - Uniform tapering, the radiation level is very high at both the 
forward sidelobes and backlobes as can be seen in Fig. B22(a). There may be a 
radiation hazard at 0 around 140° and 230°, assuming that the array is placed on the 
vehicle 2m above ground. The radiation level at the sides of the array (tp = 0° & 180°) 
is also quite high as shown in Fig. B22(b), but it drops down to below 70 dB W/m 2 for 
0 between 140° and 230°. Fig. B22(c) shows the horizontal radiation pattern at 
0 = 90°. It shows that the main lobe beamwidth is very narrow and is less than 1° at 
the -10 dB point. Similarly, for the backlobe at <p = 270°, its beamwidth is as narrow 
as the mainlobe. Hence, the radiation hazard regions at 0 around 140° and 230° can 
be expected to be very narrow. For Cosine - Uniform tapering, the problem remains at 
the front and backlobes as in the case of Uniform - Uniform tapering. Its radiation 
level at the sides of the array is very low due to the Cosine tapering along Z t direction. 
For Uniform - Cosine tapering, because Cosine tapering is used along X t direction, the 
vertical cut radiation pattern at (p = 90° shows very low sidelobe and backlobes. 
Hence, there is no possibility of causing a radiation hazard outside of the mainlobe. 
The radiation level along the sides of the array is quite high due to Uniform tapering 
along the Z t direction, but it decreases faster than that of the Uniform - Uniform 
tapering and is below 60 dB W/m 2 for 0 between 120° and 240°, therefore not causing 
any radiation hazard. 

From the above comparisons, we see that with Uniform tapering along the X t 
direction, there is very likely to be a radiation hazard due to high sidelobe and 
backlobe levels. Uniform tapering along the Z t direction also causes high radiation 
lobes at the sides of the array. But because the array area is 0.5m (or 16.7 wavelength) 
longer along the Z t direction than along the X t direction, its sidelobe level decreases 
faster and hence is less likely to cause a radiation hazard. Therefore Uniform - Cosine 
tapering type is preferred. 

Fig. B25 and B26 show the far field radiation patterns of the array using 
Uniform - Cosine tapering with the reflector plate tilted at 45° and 22.5° respectively. 
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Although the patterns are not quite as one would expect, the maximum E-field of both 
plots are equal to 0.724E4 V/m. Hence, the tilting of the reflector plate does not affect 
the gain of the array and the far field criterion is still met. 

D. VEHICLE MODELING AND FAR FIELD REQUIREMENT 

Fig. 1.1 shows the drawing of an armored vehicle carrying the Flyswatter Array. 
We need to model this vehicle to determine the effects of the presence of the vehicle on 
the radiation pattern of the array and then determine the radiation hazard region cause 
by the entire system. Fig. 4.5 and 4.6 show the computer model input to the NEC- 
BSC code, drawn by GTDGEOMP program. Note that the wheels of the vehicle are 
not shown on the model because they do not affect the radiation pattern significantly. 
The supporting poles and pivots for the array reflector plate are formed by cylinders 
while the rest of the structure is formed by plates of different shapes. Fig. 4.5 is not 
the true plot from GTDGEOMP. The hidden structures have been removed manually 
for easy comparison between the computer generated model and the actual vehicle as 
shown in Fig. 1.1. Fig. 4.6 on the other hand, shows the actual GTDGEOMP plot 
with the array rotated to face the front of the vehicle. The darkened region is the array 
area filled by small elements. 

•The full listing of the input data for Fig. 4.6 with far field pattern plot commands 
is attached in Appendix D, and the far field radiation plots are shown in Fig. B27. 
Fig. B28 shows the far field radiation plots of the same system but with the reflector 
plate tilted at 22.5°. Comparing them with Fig. B25 and B26, suggests that the 
mainlobe and the near in sidelobes have not been affected by the presence of the 
vehicular structure, and the gain, also remain unchanged. By rotating the array at 
various azimuth angles and changing the tilt angle of the reflector plate within the full 
gain range of 22.5° to 45°, the gain of the model remain unaltered. Hence, the model 
has met all the requirements specified. 

E. NEAR FIELD RADIATION HAZARD 

Before we evaluate the near field patterns, we need to consider again the accuracy 
of these plots. Unlike previous discussions, we now have many plates and cylinders 
being placed around the array. The vehicular structure plates are placed mostly below 
the array mounting plate. W'ith most of the energy radiating in a narrow beam out of 
the structure, we can be quite sure that there is not much reflected and diffracted field 
from these plates. The cylinders on the other hand are placed above the antenna 
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Figure 4.5 Computer Generated System Model with Hidden Structures Removed. 

mounting plate and are subjected to direct radiation. From the previous plot 
Fig. B24(b), we see that there is quite a considerable amount of energy radiating along 
the sides of the array. These radiated energy' will be dispersed by the cylinders. Hence, 
we can expect more higher order diffraction and reflection terms to appear and the 
accuracy of the plot is reduced. Using the NEC-BSC code guide line of 30 dB accuracy 
range, w-e will need to apply more safety margin to the radiation hazard evaluation. 

Fig. B29 and B30 show the near field radiation patterns of the full model with its 
array reflector plate tilted at 45 s and 22.5° respectively. The array radiates towards the 
front of the vehicle. From these plots, it is obvious that most of the energy is 
concentrated near the main beam. The only region that the radiation is beyond 
hazardous level is within ± 5° in the elevation angle and ± 20° in the azimuth angle 
near boresight. With the center of the array reflector erected at about 5m above the 
ground, the safety margin for operator standing at ground level is more than 27 dB. 
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Figure 4.6 Computer Generated System Model. 

Fig. B31 and B32 show the near Field radiation patterns of the Full model with the array 
rotated to point at (p = 45° and 90° respectively. The reflector plate is tilted at 45° For 
, these plots. Again, there is no sign of radiation hazard near the ground level. Hence, 
for a Uniform - Cosine tapered array, there is a safety margin of 30 dB below the 
allowable radiation hazard level. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



From the preceding investigation, we see that although Uniform tapering 
provides high gain and narrow mainbeam beamwidth, it generates very high sidelobe 
and backlobe levels and is likely to cause a radiation hazard. Cosine tapering on the 
other hand, reduces the sidelobe and backlobe levels by a considerable amount at the 
expense of lower gain and broader mainbeam beamwidth. In order to meet the far 
field requirement of producing 10-20 W/cm 2 of power density at 10 Km distance and at 
the same time be subjected to the constraint of 10.6 KV/cm of maximum aperture E- 
field strength, we are then restricted to only 3 tapering types, namely, Uniform - 
Uniform, Cosine - Uniform and Uniform - Cosine. After comparing the near field 
radiation patterns, we see that Uniform - Uniform tapering produces the highest 
sidelobe and backlobe levels and is very likely to cause a radiation hazard to the 
operators working nearby. Cosine - Uniform tapering reduces the radiation level 
towards the sides of the array. But the radiation near the front and back of the array 
is still quite high and is also likely to cause a radiation hazard. A Uniform - Cosine 
tapered array operating at 1.88 GW of input power, is capable of producing 
13.91 W/cm 2 of far field power density at 10 Km range, with a maximum aperture E- 
field strength of 10.6 KV/cm. The near field patterns show that the radiation level near 
ground level is very much lower than with the other two tapering types. With the 
vehicular structure included, we see that there is a safety margin of more than 27 dB- 
below the recommended radiation hazard limit near ground level. Although the NEC- 
BSC code suggested that the accuracy range for the radiation pattern is only up to 
30 dB from the peak value, the contributions from the higher order diffracted and 
reflected terms cannot bring the near field level up by 27 dB. Hence, we can safely say 
that there is no radiation hazard to the operator operating nearby when the Uniform - 
Cosine tapering type is used. 

From this evaluation, we see that the NEC-BSC code is very useful for 
electromagnetic radiation analysis. However, its capability in dealing with complex 
structures is very restricted by its limitation in handling only up to second order 
reflected and diffracted terms excluding doubly diffracted field. It will be able to yield 
more accurate results if higher order terms can be included. The supporting software 
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programs, due to the limited time frame, are written mainly for this project. Further 
improvement will make them fully compatible with the NEC-BSC code. These 
modifications are listed below : 

• Modify NEC-BSC code to output phase and antenna coupling values. 

• Modify GTDRWSEL and GTDPLOT codes to handle phase and antenna 
coupling output. 

• Improve the cylinder plotting algorithm of the GTDGEOMP code. 

• Improve the GTDGEOMP code to remove hidden objects. 

• Include a phase variation subroutine in the GTDSA code. 



43 



APPENDIX A 

PROGRAM FLOW CHARTS AND LISTINGS 
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b. Program Listing 



********************** 

* GTDRWSEL * 

* BY KOH WEE JIN * 

* ON MAY 26 1986 * 

* AT NPS * 

'k'k'k'k'k&'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

This program reads in the abbreviated codes and data created 
by GTD2P5 

It then selects the plots needed, convert the abbreviated codes 
into useful information and put them into a Control file. 

It also selects the appropriate data set and put them into a 
Data file 

These two file will be used by GTDPLOT for plotting 

REAL DATAIN(1000,5) 

INTEGER IPLOT (4) 

CHARACTER* 1 C0ND(4) 

CHARACTER*3 CONTRL 
CHARACTER*6 YHEAD ( 4 ) 

10 CONTINUE 

DO 30 1=1,5 

C+++ 

C+++ Initialize data set 
C+++ 

DO 20 J=l,1000 
DATAIN( J , I )=0 . 

20 CONTINUE 
30 CONTINUE 

DO 40 1=1,4 
IPLOT ( I )=0 
40 CONTINUE 
C+++ 

C+++ Read in control and comment data and 
C+++ Determine if job end is encountered 
C+++ 

CALL READC (COND,IXPARM,NPTS) 

CONTRL = COND ( 1 ) / / COND ( 2 ) / / COND ( 3 ) 

IF (CONTRL .EQ. 'END') GOTO 1000 

C+++ 

C+++ Convert abbreviated code to appropriate variables 
C+++ 

CALL CONVRT ( COND, I XP ARM, I COORD , YHEAD) 

C+++ 

C+++ Read in plotting data 
C+++ 

CALL READD (DATAIN, IXPARM,NPTS) 

50 CONTINUE 
C+++ 

C+++ Select the type of plot needed 
C+++ 

CALL CSTART (COND , I COORD , IPLOT , YHEAD) 

CALL SELECT (YHEAD , IYPARM, IPLOT , COND) 

IF (IYPARM .EQ. 5) GOTO 10 

C+++ 

C+++ Write control data to control file 
C+++ 

CALL WRITEC (NPTS,COND) 

C+++ 

C+++ Write output data to gtdpdata file 
C+++ 

CALL WRITED (DATAIN, COND ,NPTS , IYPARM) 

GOTO 50 

1000 WRITE (3, 1 (A3) ' ) CONTRL 
STOP 
END 



C 

C 

C 

C 

C 

C 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 
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c 

c 

c 

c 

c 

c 

C+++ 

C+++ 



************************* 

* SUBROUTINE CONVRT * 
************************* 

This subroutine converts the abbreviated codes 



C+++ to appropriate variables 
C+++ 

SUBROUTINE CONVRT (COND , IXPARM, ICOORD , YHEAD) 

CHARACTER*! COND(4) 

CHARACTER*3 YDIM,XAXIS ,YAXIS 

CHARACTERS COORD(3 ,2)/ 1 X 1 , ' Y' , ' Z ' , ' R 1 , 'THETA' , 'PHI'/ 
CHARACTER*6 XQ , YQ , XPARM , YPARM , YHEAD ( 4 ) 

COMMON /DPARMS/XQ , YQ . XPARM , YPARM , YDIM, XAXIS , YAXIS 
IF (COND (3) .EQ. ' R ' ) ICOORD = 1 
IF (COND (3) .EQ. 'S') ICOORD = 2 
IF ( COND ( 2 ) .EQ. 'S') THEN 
XQ = 1 F ' 

XPARM = 'FREQ' 

IXPARM = 4 
END IF 

IF (COND (2) .EQ. ' F ' ) THEN 
XPARM = 'DIST' 

IF (COND(3) .EQ. 'S' .AND. IXPARM .EQ. 2) XPARM = 'ANGLE' 

IF ( COND ( 3 ) .EQ. 'S' .AND. IXPARM .EQ. 3) XPARM = 'ANGLE' 

XQ = COORD (I XPARM, ICOORD) 

ENDIF 

DO 10 I = 1,4 

IF (I .EQ. 4) YHEAD ( I ) = COND (4)// 'PEAK' 

IF (I .NE. 4) YHEAD ( I ) = COND (4) //COORD (I , ICOORD) 

10 CONTINUE 
RETURN 
END 



C 
C 
C 
C 
C 

C+++ 

C+++ 

C+++ 



************************* 

* SUBROUTINE CSTART * 
************************* 

This subroutine displays the comment statements 
as a start for the selections that follows 



C+++ 

SUBROUTINE CSTART (COND , ICOORD , IPLOT , YHEAD) 

INTEGER IPL0T(4) 

CHARACTER*1 COND (4) 

CHARACTER* 5 COORD(3,2)/ 'X' ,'Y','Z','R' , 'THETA' , ' PHI '/ 
CHARACTER*6 YHEAD (4) 

CHARACTER*60 LABEL1 , LABEL2 , LABEL3 , LABEL4 , CNTRL 
COMMON / START/ FSTART , RXS , TYS , PZS 
COMMON /LABELS/LABEL1 , LABEL2 , LABEL3 , LABEL4 , CNTRL 
CALL FRTCMS( ' CLRSCRN ') 

PRINT * 

PRINT *, ' LABEL1 = ', LABEL 1 
PRINT * 

IF ( COND ( 1 ) .EQ. ' F ' ) PRINT *, 'FAR FIELD PATTERN PLOT' 

IF (COND(l) .EQ. 'N') PRINT *, 'NEAR FIELD PATTERN PLOT' 

PRINT * 

PRINT *, 'INITIAL CONDITIONS OF PLOTS : ' 

PRINT *, ' FREQ. = ', FSTART,' ' ,COORD(l , ICOORD) , ' = ' ,RXS 
PRINT*, ' ' ,COORD(2, ICOORD) , ' = ',TYS,' ', COORD (3 , ICOORD) , 

*' = ' ,PZS 
PRINT * 

DO 10 I = 1,4 

IF (IPLOT ( I ) .EQ. 1) PRINT *, YHEAD ( I ) , ' IS BEING PROCESSED' 
10 CONTINUE 
PRINT * 

RETURN 
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END 



C 

C 

C 

C 

C 



C+++ 

C+++ 

C+++ 



100 

110 

1000 



************************ 

* SUBROUTINE RE ADC * 

************************ 



This subroutine reads in the comments and control cards 

SUBROUTINE READC (COND, IXPARM, NPTS) 

CHARACTERS COND (4) 

CHARACTER’S CT 

CHARACTER’ 1 ' 60 LABEL1 , LABEL2 , LABEL3 , LABEL4 , CNTRL 
COMMON / START/ FSTART ,RXS , TYS ,PZS 
COMMON /LABELS/LABEL1 , LABEL2 , LABEL3 , LABEL4 , CNTRL 
READ (1,100, END=1000 ) CT,LABEL1 
FORMAT (IX, A2 , IX ,A60 ) 

IF (CT .NE. ’CT 1 ) LABEL1= 1 $ 1 

READ (1,110) COND, I XP ARM, NPTS, FSTART, RXS, TYS, PZS 
FORMAT ( 1X,4A1 , II , 14 , F7 .3 , 3F9 . 3 ) 

GOTO 20 
CONTINUE 



C+++ 

C+++ Indicate the end of the input data set 
C+++ 

COND ( 1 ) = 'E 1 
COND(2) = ' N 1 
COND ( 3 ) = ' D 1 
20 CONTINUE 
RETURN 
END 

C 

C 
C 
C 
C 

C+++ 

C+++ 

C+++ 



IF (IXPARM .EQ. 2 .AND. X(2) .EQ. 180.) ICHANG = 1 
DATAIN (1,1) = X(IXPARM) 

10 CONTINUE 
RETURN 
END 



************************ 

* SUBROUTINE READD * 

************************ 

This subroutine reads in the data 

SUBROUTINE READD (DATAIN, IXPARM, NPTS) 

REAL X(4) , DATAIN (1000, 5) 

ICHANG = 0 
DO 10 I = 1 ,NPTS 

READ (1,100) X(4),X(1),X(2),X(3),( DATAIN ( I , J ) , J=2 , 5 ) 
FORMAT ( IX , F7 .3 , 3F9 . 3 ,4E10 .3) 

IF (ICHANG .EO. 1) X(2) = 360. - X(2) 



C 
C 
C 
C 
C 

C+++ 

C+++ 

C+++ 

C+++ 



************************* 

* SUBROUTINE SELECT * 
**£*** ******* ************ 



This subroutine displays the various options 
and reads in the choices 



SUBROUTINE SELECT ( YHEAD , I YPARM , IPLOT , COND ) 

INTEGER IPL0T(4) 

CHARACTER*! COND (4) 

CHARACTERS YDIM,XAXIS, YAXIS 
CHARACTERS XQ , YQ , XPARM , YPARM , YHEAD ( 4 ) 

CHARACTER’S YUNITS (2 , 3 ) / 1 V/M ' ,'DBV/M ' , 'A/M ' 

| 'DBA/M ','W/SQM 1 ,'DBW/SQM'/ 

COMMON /DPARMS/XQ , YQ , XPARM , YPARM , YDIM , XAXIS , YAXIS 
10 CONTINUE 
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C+++ 

C+++ Select the type of plots 
C+++ 

PRINT *, 'SELECT ONE OF THE FOLLOWINGS : 1 
PRINT *, '1) PLOT ' , YHEAD(l) 

PRINT *, ’2) PLOT YHEAD (2) 

PRINT *, '3) PLOT YHE AD ( 3 ) 

PRINT *, '4) PLOT YHE AD (4) 

PRINT *, '5) NO PLOT 1 

PRINT *, 'ENTER THE NUMBER OF THE PLOT.' 

READ (5, '(II)') IYPARM 

IF (IYPARM .LT. 1 .OR. IYPARM .GT. 5) GO TO 10 
IF (COND(l) .EQ. ' F ' .AND. IYPARM .EQ. 1) THEN 
CALL FRTCMS ( ' CLRSCRN ' ) 

PRINT *, 'THERE IS NO ' , YHEAD ( 1 ) , ' FOR FAR FIELD PLOT' 

PRINT *, 'RE-ENTER PLOT SELECTION' 

PRINT * 

GOTO 10 
END IF 

IF (IYPARM .EQ. 5) GOTO 1000 
YQ=YHEAD( IYPARM) 

IF (IPLOT( IYPARM) .EQ. 1) THEN 
CALL FRTCMS ( 'CLRSCRN ') 

DO 15 1=1,4 

IF (IPLOT(I) .EQ. 1) PRINT *, YHEAD ( I ) , ' HAS BEEN PROCESSED' 
15 CONTINUE 

PRINT *, 'RE-ENTER PLOT SELECTION' 

PRINT * 

GOTO 10 
END IF 

I PLOT (IYPARM) = 1 
20 CONTINUE 
C+++ 

C+++ Select magnitude or phase plot 

C+++ (Presently, this code only caters for magnitude plot) 

C+++ 

CALL FRTCMS ('CLRSCRN ') 

PRINT * 

PRINT * 

PRINT *, 'SELECT ONE OF THE FOLLOWINGS : ' 

PRINT *, '1) MAGNITUDE PLOT' 

PRINT *, '2) PHASE PLOT' 

PRINT *, 'ENTER THE NUMBER OF THIS PLOT.' 

ISELEC = 0 

READ (5, ' (II) ' ) ISELEC 

IF (ISELEC .LT. 1 .OR. ISELEC .GT. 2) GOTO 20 
IF (ISELEC .EQ. 1) YPARM = 'MAG' 

IF (ISELEC .EQ. 2) YPARM = 'PHASE' 

C+++ 

C+++ Select the Y-axis' unit. 

C+++ 



IF 


(CONDI 


14) 


.EQ. 


' E ' ] 


1 IYUNIT = 1 


IF 


(CONDI 


4) 


.EQ. 


1 H ' 


1 IYUNIT = 2 


IF 


(CONDI 


4) 


.EQ. 


■ p ■ 3 


1 IYUNIT = 3 



30 CONTINUE 

CALL FRTCMS ( 'CLRSCRN ') 

PRINT * 

PRINT * 

PRINT *, 'SELECT THE UNITS TO BE USED.' 
PRINT*, '1) ' ,YUNITS(1, IYUNIT) 

PRINT*, '2) ' ,YUNITS( 2, IYUNIT) 

PRINT *, 'ENTER THE NUMBER OF YOUR CHOICE.' 
ISELEC = 0 

READ (5, ' (II) ' ) ISELEC 

IF (ISELEC .LT. 1 .OR. ISELEC .GT. 2) GOTO 30 
IF (ISELEC .EQ. 1) YDIM = 'LIN' 

IF (ISELEC .EQ. 2) YDIM = ' DB ' 

40 CONTINUE 
C+++ 
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C+++ Select linear or log scale plot 
C+++ 



CALL FRTCMS ( 
PRINT * 



CLRSCRN ') 



PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
ISELEC = 
READ ( 5 , 1 



'SELECT 
■1) X s 
■2) X : 
•3) X : 
'4) X : 
' ENTER 
0 

(ID 1 ) 



IF (ISELEC .LT. 
XAXIS = 'LIN' 
YAXIS = 'LIN' 



1000 



IF 

IF 



(ISELEC 

(ISELEC 



CONTINUE 

RETURN 

END 



THE TYPES 
LINEAR, Y 
LINEAR, Y 
LOG, Y 

LOG, Y 

THE NUMBER 



OF PLOT : 
: LINEAR' 
: LOG 1 
: LINEAR' 
: LOG 
OF YOUR 



CHOICE. ' 



ISELEC 
1 .OR. 



ISELEC .GT. 4) GOTO 40 



.EQ. 

.EQ. 



.OR. 

.OR. 



ISELEC 

ISELEC 



.EQ. 

.EQ. 



YAXIS = 
YAXIS = 



'LOG' 

'LOG' 



C 
C 
C 
C 
C 

C+++ 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 



************************* 

* SUBROUTINE- WRITEC * 
************************* 



This subroutine prepares and writes the various labels 
and fields needed to create the control file 

SUBROUTINE WRITEC (NPTS,COND) 

CHARACTER*! C0ND(4) 

CHARACTER*3 XAXIS , YAXIS ,YDIM 
CHARACTER*4 UNITS$ 

CHARACTER* 5 NPT$ 

CHARACTER*6 RXC$ , TYC$ , PZC$ , XQ , YQ , XPARM , YPARM , FIELD$ 
CHARACTER*? FS$ 

CHARACTER*9 RX$,TY$,PZ$ 

CHARACTER*60 LABEL1 , LABEL2 , LABEL3 , LABEL4 , CNTRL 
COMMON /START/ FSTART , RXS , TYS , PZS 

COMMON /DP ARMS / XQ , YQ , XPARM , YPARM , YD IM , XAXI S , YAXIS 
COMMON /LABELS/LABEL1 , LABEL2 , LABEL3 , LABEL4 , CNTRL 

Converts numerical variables to characters 



• C+++ 

OPEN (4 , STATUS= ' NEW ' ) 

WRITE (4,100) NPTS, FSTART, RXS, TYS, PZS 
100 FORMAT ( IX, 14 , IX, F7 .3 , 3F9 .3) 

CLOSE (4) 

OPEN (4 , STATUS= ' OLD ' ) 

READ (4,110) NPT$ ,FS$ , RX$ ,TY$ , PZ$ 

110 FORMAT ( IX, A5,A7,3A9) 

CLOSE ( 4 , STATUS= 1 DELETE ' ) 

C+++ 

C+++ Create the control statement 
C+++ 

CNTRL = 'NEW ' //NPT$//XQ// ' '//XPARM//' '//YQ//' '//YPARM//' 
*YDIM// ' '//XAXIS//' '//YAXIS 

C+++ 

C+++ Prepare the labels 

C+++ Label 1 displays the title of the plot 

C+++ Label 2 to 4 display the initial conditions 

C+++ 

UNITS$ = ' M ' 

RXC$ = ' (X)= ' 

TYC$ = ' (Y)= ' 
pzc$ = ' ( z ) = ' 

IF ( COND ( 3 ) .EQ. ' R ' ) GOTO 10 
UNITS$ = ' DEG' 



// 
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RXC$ = 1 ()R = 1 
TYC$ = ' ?Q) = 1 
PZC$ = 1 ?F) = 1 
10 CONTINUE 

IF (COND(l) .EQ. ' F ' ) FIELD$ = 1 (F)AR' 

IF (COND(l) .EQ. ' N ' > FIELD$ = 1 (N)EAR' 

LABEL2=FIELD$//' FIELD (GTD) PLOT WITH STARTING COND. AS FOLLOWS$' 
LABEL3= ' FREQ = '//FS$//' (GH)Z 1 //RXC$//RX£// ' M$ 1 

LABEL4= 1 '//TYC$//TY$//UNITS$//' ' //PZC$//PZ^//UNITS$// ' $ 1 

c++*+ 

C+++ Writes the control statement and labels onto the control file 



C+++ 



WRITE 1 


;3, 1 (ago; 


1 1 ) 


WRITE < 


3, 1 (A60 


1 1 J 


WRITE 1 


3, 1 (A60 


1 1 J 


WRITE 1 


3, 1 (A60 


| 1 \ 


WRITE ( 


3, 1 (a6o; 


) ' ) 



RETURN 

END 



CNTRL 

LABEL1 

LABEL2 

LABEL3 

LABEL4 



C- 

C 

c 

c 

c 



* SUBROUTINE WRITED * 
************************* 



C+++ 

C+++ This subroutine writes the appropriate data to GTDPDATA file 
C+++ 

SUBROUTINE WRITED (DATAIN, COND ,NPTS , IYPARM) 

REAL DATAIN (1000, 5) 

CHARACTERS C0ND(4) 

CHARACTERS YDIM,XAXIS,YAXIS 
CHARACTERS XQ , YQ , XP ARM , YPARM- 

COMMON / DPARMS / XQ , YQ , XP ARM , YP ARM , YD I M , XAXI S , YAXI S 
DUMMY =0. 

IY = IYPARM + 1 
DO 10 I = 1 ,NPTS 

IF (YDIM .EQ. 'DB ') THEN 

IF (DATAIN (I , IY) .LE. l.E-20) DATAIN (I , IY)=1 .E-20 
IF ( COND ( 4 ) .EQ. ' E ' ) DATA=20 .*AL0G10 (DATAIN (I , IY) ) 
IF ( COND ( 4 ) .EQ. ' H 1 ) DATA=20 .*ALOG10( DATAIN (I , IY) ) 
IF ( COND ( 4 ) .EQ. ' P ' ) DATA=10 .*ALOG10( DATAIN (I , IY) ) 

ELSE 

DATA=DATAIN ( I , I Y) 

ENDIF 

WRITE (2,100) DATAIN (1,1), DATA, DUMMY 
100 FORMAT (5X,3F10. 3) 

10 CONTINUE 
RETURN 
END 
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OUTPUT DATA PLOTTING CODE (GTDPLOT) 
a. Flow Chart 



Mafn Program 
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Flow Chart For GTDPLOT (continue) 
Subroutine PARMS 




return 



b. Program Listing 



c 

c 

c 

c 

c 

c 

c 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

1010 CONTINUE 
C+++ • 

C+++ Read control statements and labels to be plotted. 

C+++ 

CALL PARMS (END) 

IF (END) GO TO 1020 

C+++ 

C+++ Plot the graph 
C+++ 

CALL PLOT 

C+++ 

C+++ Indicate that we have completed at least one plot. 

C+++ 

STRTED = .TRUE. 

GO TO 1010 

C+++ 

C+++ End of plot. Otherwise, program abort through an error handling 
C+++ routine. 

C+++ 

1020 CONTINUE 

CALL DISEND 

STOP 

END 



kkkkkkkkkkkkkkkkkkkkkkkkkkkk 

* GTDPLOT * 

* MODIFIED FROM NFPLOT * 

* BY KOH WEE JIN * 

* ON MAY 29 1986 * 

* AT NPS * 

**************************** 

This program reads in the Control and Data files created by 
GTDRWSEL . 

It then select the appropriate labels and scales and plot types 
and plot them using Disspla subroutines 

LOGICAL*4 STRTED 
L0GICAL*4 END /.FALSE./ 

COMMON /DISPLA/ STRTED 

Turn off the start indicator 

STRTED = .FALSE. 

Initialize Disspla 

CALL DISSET 

Plot until 'END' control statement is detected 



c 

C 

C 

C 

C 

C 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



*********************** 

* SUBROUTINE ALGPLT * 

k k k k k kk -k k k * A kk k k k -k k k k k k 

This subroutine replaces the Disspla subroutine ALGPLT. 

It provides a rounded axis origin and log cycle length in 
inches per cycle. 

Modification is provided to eliminate truncated axes symptomatic 
in logarithmic plots, only full log cycles are returned 

SUBROUTINE ALGPLT (AMIN, AMAX, AXLEN, A0, ACY) 

REAL 5 ^ AMIN 
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REAL*4 


AMAX 


REAL*4 


AXLEN 


REAL*4 


A0 


REAL*4 


ACY 


REAL*8 


DMIN 


REAL*8 


DMAX 


REAL*8 


DXLEN 


REAL*8 


DAO 


REAL*8 


OMAGO 


REAL*8 


OMAG1 


INTEGER*4 


MAGDIF 



C+++ 

C+++ Do it all in double precision 
C+++ 

DMIN = AMIN 
DMAX = AMAX 
DXLEN = AXLEN 

C+++ 

C+++ Get the order of magnitude of the minimum value 
C+++ 

OMAGO = DLOGIO(DMIN) 

C+++ 

C+++ If the minimum is an even power of 10, then use it as the lower 
C+++ limit. Otherwise, go to next lower order of magnitude. 

C+++ 

IF ( IDINT (OMAGO ) .NE. OMAGO) OMAGO = IDINT (OMAGO - 1.) 

C+++ 

C+++ Do the same for the maximum. 

C+++ 

OMAG1 = DLOGIO (DMAX) 

C+++ 

C+++ If the minimum is an even power of 10, then use it as the upper 
C+++ limit. Otherwise, go to next higher order of magnitude. 

C+++ 

IF ( IDINT (0MAG1) .NE. .OMAG1 ) 0MAG1 = IDINT(0MAG1 +1.) 

C+++ 

C+++ Compute the origin 
C+++ 

1010 CONTINUE 

DAO = 10 . DO** IDINT (OMAGO) 

IF (DAO .GT. DMIN) THEN 

OMAGO = IDINT (OMAGO) - 1. 

GO TO 1010 
ENDIF 

C+++ 

C+++ Compute the difference in order of magnitude and the number of 
C+++ user units (usually in inches) per cycle of log. 

C+++ 

1020 CONTINUE 

MAGDIF = (OMAG1 - OMAGO) - 1. 

IF ((DAO * 10 . D0**MAGDIF ) .LT. DMAX) THEN 
OMAG1 = OMAG1 + 1. 

GO TO 1020 
ENDIF 
A0 = DAO 

ACY = DXLEN / MAGDIF 

RETURN 

END 



C 
C 
C 
C 
C 

C+++ 

C+++ 

C+++ 

C+++ 



*********************** 

* SUBROUTINE BLSRCH * 
*********************** 



This subroutine scans for the first non-blank character and 
returns with the pointer (I) set to point to the character. 

SUBROUTINE BLSRCH(LINE , I) 

CHARACTER*! LINE (80) 
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CHARACTERS BLANK /■ '/ 

C+++ 

C+++ If we have gone beyond the end of the line, complain and abort. 
C+++ 

IX = I 

IF (I .GE. 81) GO TO 8010 

C+++ 

C+++ Scan for a non-blank 
C+++ 

DO 1010 I = IX, 80 
IF (LINE(I) .NE. BLANK) RETURN 
1010 CONTINUE 



C+++ 

C+++ Error message. 

C+++ 

8010 CONTINUE 

CALL ERRMSG(IX, 'Unexpected end of line detected. $') 
CALL DISEND 
STOP 400 
END 



C 
C 
C 
C 
C 

C+++ 

C+++ 

C+++ 



9010 



*********************** 

* SUBROUTINE DISEND * 
*********************** 



End of Disspla 

SUBROUTINE DISEND 
LOGICAL*4 STRTED 
COMMON /DISPLA/ STRTED 
IF (STRTED) CALL DONEPL ' 

IF (.NOT. STRTED) WRITE(6,9010) 

RETURN 

FORMAT ('0N0 PLOTS GENERATED... PROGRAM ABORT.' 
END 



) 



C 

C 

c 

c 

c 



C+++ 

C+++ 

C+++ 



cxxx 

C+++ 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



*********************** 



* SUBROUTINE DISSET * 
*********************** 



This subroutine establishes the parameters for Disspla. 

SUBROUTINE DISSET 
CALL COMPRS 

CALL SETCPR(0 , 1, 0, 0) 

CALL PAGE (8 . 5 , 11.0) 

Lettering is FUTURA with upper case and lower case. 

Upper case letters are delimited by '(' escape character. 

CALL BLOWUP (.5) 

CALL FUTURA 

CALL BASALF ( ' L/CSTANDARD ' ) 

CALL MIXALF ( 1 STANDARD ' ) 

CALL MX3ALF ( 1 INSTRUCTION' , ' | ' ) 

CALL MX4ALF( 'L/CGREEK' , '?') 

CALL NOBRDR 

Integer (or rounded) axes with Y axis labelling at 0 degrees. 



C 

C 

C 



CALL INTAXS 
CALL YAXANG(0. ) 
RETURN 
END 



*********************** 
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* SUBROUTINE ERRMSG * 

it -k'k'k'k'k'k-k'k 



Cl 

c| 

C+++ 

C+++ This subroutine writes an error message to the logical unit 6. 
C+++ 

SUBROUTINE ERRMSG (N , MSG) 

CHARACTER*! DOLLAR /'$'/ 

CHARACTER* 1 MSG(l) 

CHARACTER* 1 BAR / 1 '/ 

CHARACTER*! BLANK 7' 7 

C+++ 

C+++ Indicate the error by writing the bar below it. 

C+++ 

WRITE(6,9010) (BLANK, 1=2, N) , BAR 

C+++ 

C+++ Find the end of the message. 

C+++ 

DO 1010 I = 1, 80 

IF (MSG(I) . EQ. DOLLAR) GO TO 1020 
1010 CONTINUE 
I = 81 



C+++ 

C+++ Write the message. 

C+++ 

1020 CONTINUE 

1 = 1-1 

WRITE (6 , 9010) (MSG(J), J = 1, I) 
RETURN 

9010 FORMAT ( IX , 80A1 ) 

END 



C 
C 
C 
C 
C 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



1010 

1020 



*********************** 

* SUBROUTINE EXPONE * 
*********************** 



This subroutine concatenates the exponent (if required). 

SUBROUTINE EXPONE ( INVEC , * EXPON) 

CHARACTER*! INVEC (64) 

INTEGER*4 EXPON 

INTEGER*4 . FILL /Z000000F0/ 

CHARACTER* 1 HEXCHR(4) 

INTEGER*4 DUMMY 

EQUIVALENCE (HEXCHR(l) , DUMMY) 

IF (EXPON .EQ. 0) RETURN 

Find the end of the input vector and generate the exponent. 

DO 1010 I = 1. 64 

IF (INVEC (I ) .EQ. '$') GO TO 1020 
CONTINUE 

CALL ABEND ( 'EXPONE' , 'Unable to concatenate. Vector full. 
CONTINUE 
INVEC(I) = ' ' 

1 = 1 + 1 

CALL TIMS10(INVEC(I) , I) 

IXP = EXPON 
IF (IXP .LT. 0) THEN 
IXP = IABS(IXP) 

INVEC(I) = '-' 

1 = 1 + 1 
ENDIF 

IXX = IXP / 10 
IF (IXX .NE. 0) THEN 
DUMMY = IXX + FILL 
INVEC ( I ) = HEXCHR(4) 

1 = 1 + 1 

IXP = IXP - (IXX * 10) 



$' 



) 
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c 
c 
c 
c 
c 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



ENDIF 

DUMMY = IXP + FILL 
INVEC(I) = HEXCHR(4) 
1 = 1 + 1 
INVEC(I) = '$' 

RETURN 

END 



*********************** 

* SUBROUTINE GETTOK * 
*********************** 



This subroutine extracts an 8 character token. 

SUBROUTINE GETTOK(LINE , I, TOKEN) 

CHARACTERS LINE (80) 

CHARACTERS TOKEN 
CHARACTERS TOK 
CHARACTER* 1 CHAR(8) 

EQUIVALENCE (CHAR(l), TOK) 

CHARACTERS BLANKS /' '/ 

CHARACTERS BLANK /' '/ . 

CHARACTERS COMMA / ' , ' / 

Clear the token. 

TOK = BLANKS 

Find the first non-blank character. 

CALL BLSRCH(LINE , I) 

Extract the token (up to 8 characters). 



N = 1 
IX = I 

DO 1010 I = IX, 80 
IF (LINE(I) .EQ. BLANK) GO TO 1030 
IF (LINE ( I ) .EQ. COMMA) GO TO 1020 
IF (N .EQ. 9) GO TO 8010 

C+++ 

C+++ Character Ok. So put it into the token buffer. 
C+++ 

CHAR(N) = LINE (I) 

N = N + 1 
1010 CONTINUE 
I = 81 
RETURN 

C+++ 

C+++ Exit if comma found after bumping the pointer. 
C+++ 

1020 CONTINUE 

1 = 1 + 1 
1030 CONTINUE 

TOKEN = TOK 
RETURN 



C+++ 

C+++ Error, token length is > 8. 

C+++ 

8010 CONTINUE 

CALL ERRMSG( I , 

'Command word or token length exceeds 8 characters .$ 1 ) 
CALL DISEND 
STOP 500 
END 



C 

C 

C 



*&**** jtjs'kjt'k'k'k'k'k'k'k'k'kJt'k'k'k 
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* SUBROUTINE LINLIN * 
*********************** 



£111 

c+++ 

C+++ 

C+++ 



c 
c 
c 
c 
c 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



c 

c 

c 

c 

c 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



c 

c 

c 

c 

c 



This subroutine sets up DISSPLA for a Linear-Linear axis plot. 

SUBROUTINE LINLIN 
REAL *4 XMIN 

REAL*4 XMAX 

REAL*4 YMIN 

REAL*4 YMAX 

COMMON /RANGE/ XMIN, XMAX, YMIN, YMAX 

CALL GRAF (XMIN, 'SCALE', XMAX, YMIN, 'SCALE', YMAX) 

RETURN 

END 



^^ii^^iz'k’k'k'k'k'k'k'k'k'k'k'k'k'k'k'k’k 

* SUBROUTINE LINLOG * 
*********************** 



This subroutine sets up DISSPLA for a Linear-Log axis plot. 

SUBROUTINE LINLOG 
REAL*4 XMIN 

REAL*4 XMAX 

REAL*4 YMIN 

REAL*4 YMAX 

COMMON /RANGE/ XMIN, XMAX, YMIN, YMAX 

Get the proper scaling for both axes. 

CALL AXSPLT(XMIN, XMAX, 6.0, XORIG, XSTEP, DUMMY) 

CALL ALGPLT(YMIN, YMAX, 6.0, YORIG, YCYCLE) 

Set up DISSPLA for a Log scale on the Y axis. 

CALL YLOG( XORIG , XSTEP, YORIG, YCYCLE) 

RETURN 

END 



************* ******** 

* SUBROUTINE LGLG * 
********************* 



This subroutine sets up DISSPLA for a Log-Log axis plot. 

SUBROUTINE LGLG 
REAL*4 XMIN 

REAL*4 XMAX 

REAL*4 YMIN 

REAL*4 YMAX 

COMMON /RANGE/ XMIN, XMAX, YMIN, YMAX 

Get the proper scaling for both axes. 

CALL ALGPLT(XMIN, XMAX, 6.0, XORIG, XCYCLE) 

CALL ALGPLT(YMIN, YMAX, 6.0, YORIG, YCYCLE) 



Set up DISSPLA for a Log-Log scale plot. 



CALL LOGLOG (XORIG, XCYCLE, 

RETURN 

END 



YORIG, YCYCLE) 



*********************** 

* SUBROUTINE LOGLIN * 
*********************** 
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C+++ 

C+++ This subroutine sets up DISSPLA for a Log-Linear axis plot. 
C+++ 

SUBROUTINE LOGLIN 
REAL*4 XMIN 

REAL*4 XMAX 

REAL*4 YMIN 

REAL*4 YMAX 

COMMON /RANGE/ XMIN, XMAX, YMIN, YMAX 

C+++ 

C+++ Get the proper scaling for both axes. 

C+++ 

CALL ALGPLT(XMIN, XMAX, 6.0, XORIG, XCYCLE) 

CALL AXSPLT(YMIN, YMAX, 6.0, YORIG, YSTEP , DUMMY) 

C+++ 

C+++ Set up DISSPLA for a Log scale on the X axis. 

C+++ 

CALL XLOG( XORIG, XCYCLE, YORIG, YSTEP) 

RETURN 

END 



*********************** 

* SUBROUTINE NUMBER * 

*********************** 

This subroutine extracts an integer from the data stream. 

SUBROUTINE NUMBER (LINE , I, N) 

CHARACTER* 1 LINE (80) 

CHARACTER* 1 T0KEN(8) 

CHARACTER* 1 ZERO /'O'/ 

CHARACTER* 1 NINE / 1 9 ' / 

CHARACTER* 1 BLANK /' ■/ 

INTEGER*4 ZIP 

CHARACTER* 1 NONE (4) /3*Z00,'0'/ 

EQUIVALENCE (NONE(l), ZIP) 

INTEGER*4 NUM 

CHARACTER*! BYTES (4) /4*Z00/ 

EQUIVALENCE (NUM, BYTES) 

Extract the token. 

CALL GETTOK(LINE , I, TOKEN) 

Get the number. 

N = 0 

DO 1010 J = 1, 8 
IF (TOKEN(J) .EQ. BLANK) RETURN 

IF ((TOKEN(J) .LT. ZERO) .OR. (TOKEN (J) .GT. NINE)) GO TO 8010 

Extract the digit and get it into the number. 

BYTES (4) = TOKEN(J) 

N = (N * 10) + (NUM - ZIP) 

1010 CONTINUE 
RETURN 

C+++ 

C+++ Error. Not a digit. 

C+++ 

8010 CONTINUE 

CALL ERRMSG(I, 'Not a digit.$') 

CALL DISEND 

STOP 

END 

C 

C ********************** 



c 

c 

c 

c 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 
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* SUBROUTINE PARMS * 
********************** 



C 

c 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



This subroutine reads control parameters and data 
The control statement syntax is: 

line 1 — 



<STATUSXNPTSXX QXX PARMXY QXY PARMXY UNITXX AXISXY AXIS> 



where <STATUS> 
<NPTS> 

<X Q> 

<X PARM> 
<Y Q> 



<Y PARM> 

<Y UNIT> 
<X AXIS> 
<Y AXIS> 



is 'OLD' , 'NEW 1 , or 'END' , 

is the number of points to be plotted (omit if 
'OLD 1 of 'END' is selected), 

is 'X', 'Y', 'Z 1 , 'R', 'THETA', 'PHI' or 'FREQ' 
is 'DIST', 'MAG', 'ANGLE' or 'FREQ' 
si 'ER', 'EX', 'EPHI', 'EY', 'ETHETA', ' EZ ' , 
'EPEAK', 'HR 1 , 'HX', 'HPHI', 'HY', ' HTHETA" , 
'HZ', 'HPEAK', 'PR', 'PX', 'PPHI', 'PY', 

' PTHETA ' , 'PZ' or ' PPEAK ' 

is 'DIST', 'MAG', or 'PHASE' (<X PARM> and 
<Y PARM> are mutually exclusive), 

is 'LIN' or ' DB ' 

is 'LIN' or 'LOG', and 

is 'LIN' or 'LOG'. 



line 2 to 5 are <LABEL> ' s 
where <LABEL> 



is a label up to 60 characters (letters to 
be capitalized must be enclosed in parentheses) 
terminated by the '$' character. 



SUBROUTINE PARMS (END) 

REAL *4 DATA\1000,3) 



INTEGERS 
INTEGERS 
INTEGERS 
CHARACTER*1 
CHARACTER* 1 
CHARACTERS 
CHARACTER*! 
INTEGER*4 
INTEGER*4 
INTEGER*4 
INTEGER*4 
INTEGER*4 



NPTS 
PI 
P2 

LABEL! 

LABEL2 
LABEL3 
LABEL4 
XTYPE 
YTYPE 
XQUAL 
YgUAL 

COMMON /DPARMS/ DATA, NPTS, PI 
LABEL4 , XTYPE, 
END 

CTL(80) 

SELECT(5) 



L0GICAL*4 

CHARACTER*! 

CHARACTER*8 



, P2 , LABEL1, LABEL2, LABEL3 , 
YTYPE, XQUAL, YQUAL, YU 



/'DIST' , 
'MAG' , 

' PHASE ' 

' ANGLE ' 

' FREQ ' / 

CHARACTER*8 XQUALS(7) /'X', 

' Y' , 

' Z* , 

'R' , 

' PHI ' , 

' THETA ' 

' F ' / 

YUNIT(2) /'LIN', 
'DB 1 / 

YQUALS(21) /'ER', 

' EPEAK ' 
'HR', 1 
'HPEAK' 
'PR 1 , ' 

'PPEAK' 



CHARACTER*8 
CHARACTER* 8 



'EX' 

HX' 

PX' 

/ 



ETHETA ' , 


'EY' , 


'EPHI' , 


'EZ 


HTHETA ' , 


'HY' , 


' HPHI ' , 


'HZ 


PTHETA ' , 


'PY' , 


'PPHI' , 


'PZ 
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CHARACTER*8 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 

1010 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 

1020 



C+++ 

C+++ 

C+++ 

1030 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



COMAND ( 3 ) /'NEW', 

' OLD ' , 

'END '/ 

CHARACTERS TYPES (2) /'LIN' 

| ' LOG ' / 

LOGICAL*4 DATAIN /.FALSE./ 

Read the control information and echo it. 

READ (5 ,9010 ,END=8010) CTL 
WRITE (6, 9020) CTL 

Parse the control information 

1 = 1 



Data source: 

OLD — use data already read in. 

NEW — read a new block of <NPTS> points 
END — end of job 

CALL XTRACT(CTL , I, COMAND, 3, JUMP) 

GO TO (1010, 1030, 1020), JUMP 

'NEW' selected. 



CONTINUE 
INUMER = I 

CALL NUMBER (CTL, I, NPTS) 

only 1000 or less points allowed. 

IF (NPTS .GT. 10C0) GO TO 8020 
CALL READIN 
GO TO 1030 



'END' selected. 

CONTINUE 
END = .TRUE. 

RETURN 

Process the X and Y PARM specifications 
CONTINUE 

CALL XTRACT (CTL , I, XQUALS, 7, XQUAL) 
CALL XTRACT (CTL, I, SELECT, 5, PI) 

CALL XTRACT (CTL, I, YQUALS , 21, YQUAL) 
CALL XTRACT (CTL, I, SELECT, 5, P2) 

PRINT * , ' XQUALS ( XQUAL ) = ' , XQUALS ( XQUAL ) 
PRINT *, ' YQUALS ( YQUAL )=', YQUALS (YQUAL) 
PRINT *, 'XSELECT(P1)=' ,SELECT(P1) 

PRINT *, ' YSELECT(P2)= ' ,SELECT(P2) 

Select the Y axis unit. 



CALL XTRACT (CTL , I , YUNIT , 2 ,YU) 
PRINT*, 'YUNIT(YU)=' ,YUNIT(YU) 



Select the axis type 



CALL XTRACT (CTL 
CALL XTRACT ( CTL 
PRINT *, 'TYPES 
PRINT *, 'TYPES 



, I, TYPES, 2, XTYPE ) 
, I, TYPES, 2, YTYPE ) 
(XTYPE )=' , TYPES (XTYPE 
(YTYPE )=' , TYPES (YTYPE 



! 



Read in the labers 
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C+++ 

C+++ 

C+++ 

8010 



8020 



8030 

9010 

9020 



C 

c 

c 

c 

c 



C+++ 

C+++ 

C+++ 



READ ( 5 , 9010 ,END=8030) LABEL1 

WRITE (6 ,9020) LABEL 1 

READ ( 5 , 9010 , END=8030) LABEL2 

WRITE (6 ,9020) LABEL2 

READ ( 5 , 9010 ,END=8030) LABEL3 

WRITE(6 , 9020) LABEL3 

READ (5 , 9010 ,END=8030) LABEL4 

WRITE (6 ,9020) LABEL4 

RETURN 

Error messages 
CONTINUE 

CALL ERRMSG(1, 'No END statement found.$') 

CALL DISEND 
STOP 100 
CONTINUE 

CALL ERRMSG(INUMER, 'Number of points exceeds 1000.$') 
CALL DISEND 
STOP 200 
CONTINUE 

CALL ABEOF('PARMS ', 5) 

FORMAT (80A1) 

FORMAT ( IX, 80A1) 

END 



********************* 
* SUBROUTINE PLOT * 

-k k k k k k k * * * * * k * * * * * * * * 



This subroutine plots the data. 



SUBROUTINE 
REAL*4 
INTEGERS 
INTEGER*4 
INTEGERS 
CHARACTER*1 
CHARACTER*1 
CHARACTER*! 
CHARACTER* 1 
INTEGER*4 
INTEGER*4 
INTEGER*4 
INTEGER*4 
INTEGER*4 



PLOT 

DATA( 1000 , 3 ) 

NPTS 

PI 

P2 

LABEL1 (80' 
LABEL2(80 
LABEL3(80 
LABEL4(80 j 
XTYPE 
YTYPE 
XQUAL 
YQUAL 
YU 



COMMON /DP ARMS/ DATA, NPTS, PI 
LABEL4, XTYPE, 
EXPONX 



INTEGER*4 
INTEGER*4 
REAL*4 
REAL*4 
REAL*4 
REAL*4 
COMMON 
CHARACTER*64 
CHARACTER*12 



P2 , LABEL 1 , LABEL2, 
YTYPE, XQUAL, YQUAL, 



EXPONY 
XMIN 
XMAX 
YMIN 
YMAX 

/RANGE/ XMIN, 
NAMVEC 
NAME (5) 



CHARACTER*16 XUNITS(7) 



XMAX, YMIN, YMAX 



I STANCE 
AGNITUDE 
HASE 1 , 
NGLE ' , 
REQUENCY 

x 





LABEL3 , 
YU 



62 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



CHARACTER*40 YUNITS(21,2 
DATA ( (YUNITS ( I , J) , 1=1 , 7 
iX LXHX) ((V)/M 
IY LXHX 
iZ LXHX 
IR LXHX 
iX LXHX 



E 
E 
E 
E 
E 
E 
E 
E 
E 
E 

DATA 



LH.75 

LH.75 

LH.75 

LH.75 

LH.75 

LH.75 

LH.75 

LH.75 



V)/M 

V)/M 

D(BV 




LH. 75)R | LXHX) { 



•?Q) ILXHXF8 
'YfLXHX) (( 



/M( )$ 



LXHXF8) ( 
, LXHXF8) ( 
(PEAK) | LXHX) 



V)/M{ 

ywiu 



?F) | LXHXF8 

LH.75 )Z| LXHX) (()D 
LH.75') (PEAK) I LXHX) 
( (YUNITS ( I, J) ,1=8 
H LH . 7 5 ) X LXHX) ((A 
H LH . 7 5 ) Y LXHX 
H LH . 7 5 ) Z LXHX 
H LH . 7 5 ) R LXHX 
H LH . 7 5 ) X LXHX 
H LH.75)?Q) | LXHXF8 
H LH.75)YfLXHX) (( 

H LH . 7 5 > ? F ) | LXHXF8 
H LH . 7 5 ! Z | LXHX ) (( 

H LH.75) (PEAK) II"" 
DATA ‘ 



BV 
()D(BV 
BV)/M( 



5 )$ 1 




/iN M $', 

. H( ) ) $ 1 , 

D(BV)/M( ) )$ 1 / 

, J=1 ,2) / 1 (Hi LH.75 )R | LXHX) 

LH . 7 5 ) ?Q ) I LXHXF8 ) (i 

LH . 7 5 ) ? F ) I LXHXF 8 ) (i 
LH.75) (PEAK) |LXHX) 

7m())$' , 

BA)/M( ) )$ 



.(A)/M( 

(A)/M( 




BA)/M( 
.. X)(()D- 

( (YUNITS ( I , J) , 1=15 , 21 ) 




LH.75 

LH.75' 

LH.75' 

LH.75 

LH.75 

LH.75 

LH.75’ 

LH.75 

LH.75 

LH.75 

LH.75 

LH.75 

LH.75 

LH.75 



(PEAK) 1 LXHX) ( ( )D(BA) /M( ) ) $ 1 
- -- -* jsl.2)/ 






R I LXHX) 

X | LXHX) 

?Q) | LXHXF8) 

YfLXHX) ((W)/SQ ... 
? F ) | LXHXF8) ((W)/S 
Z | LXHX) ( (W)/SQ M ’ 
(PEAK) LXHX) 

R I LXHX) ( ( )D 
X| LXHX) (( 

?Q) | LXHXF8 
YfLXHX) (( 

?F) | LXHXF8 
Z | LXHX) (( 



. f 

(tW)/SQ M( 
BW)/SQ Mj|, 



)D(BW)/S 



(.).P(??)/S^ ) M0)$ I 



(PEAK) | LXHX) 



BW) /SO w V//v . 

()D(BV?)/SQ M(5 >$‘ , 

BW)/SQ MU )$• , 
(()D(BW)/SQ M ())$' 



Find minimum and maximum 



IP1 = 1 

CALL MINMAX(DATA( 1 , IP1 ) , NPTS , XMIN, XMAX) 

IP2 = 2 

CALL MINMAX(DATA( 1 , IP2 ) , NPTS, YMIN , YMAX) 

PRINT*, 1 XMIN= 1 , XMIN, ' XMAX=',XMAX 

PRINT *, 1 YMIN= 1 ,YMIN, 1 YMAX= 1 , YMAX 

PRINT*, 'PI = 1 ,P1 , 1 P2 = 1 ,P2 

PRINT*, 'NAME(P1)=' ,NAME(P1) 

PRINT*, , NAME(P2)=' ,NAME(P2) 

CALL AREA2D(6 . , 6. ) 

Write the headers 

CALL HEADIN(LABEL1 , 100, 1.0, 4) 

CALL HEADIN ( LABEL2 , 100, 1.0, 4) 

CALL HEADIN (LABEL3, 100, 1.0, 4) 

CALL HEADIN (LABEL4, 100, 1.0, 4) 

Label the axes 

NAMVEC=NAME ( PI ) / / 1 1 / /XUNITS ( XQUAL ) 

PRINT*, 1 NAMVEC=' ,NAMVEC 

CALL SCLDAT(DATA(l , IP1 ) , NPTS, EXPONX, XMIN, XMAX) 
CALL EXPONE (NAMVEC , EXPONX) 

PRINT * , ' XNAME= ' , NAMVEC 
CALL XNAME (NAMVEC, 100) 

NAMVEC=NAME (P2 ) // ' 1 //YUNITS ( YQUAL , YU) 

PRINT*, 'NAMVEC=' , NAMVEC 
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C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



c 
c 
c 
c 
c 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



CALL SCLDAT(DATA(1,IP2), NPTS, EXPONY, YMIN, YMAX) 
CALL EXPONE (NAMVEC , EXPONY) 

PRINT * , 'YNAME=' , NAMVEC 
CALL YNAME (NAMVEC, 100) 

PRINT * , ' XTYPE= 1 , XTYPE , ' YTYPE= ' , YTYPE 

Dispatch to the correct routine. 

IF ((XTYPE .EQ. 1) .AND. (YTYPE .EQ. 1)) CALL LINLIN 

IF ((XTYPE .EQ. 1) .AND. (YTYPE .EQ. 2)) CALL LINLOG 

IF ((XTYPE .EQ. 2) .AND. (YTYPE .EQ. 1)) CALL LOGLIN 

IF ((XTYPE .EQ. 2) .AND. (YTYPE .EQ. 2)) CALL LGLG 

Draw the curve 

CALL CURVE (DATA(1 , IP1) . DATA(1,IP2), NPTS, 0) 

CALL RSCDAT(DATA(1 , IP1 ) , NPTS, EXPONX, XMIN, XMAX) 
CALL RSCDAT(DATA( 1 , IP2) , NPTS, EXPONY, YMIN, YMAX) 
CALL ENDPL(O) 

RETURN 

END 



*********************** 

* SUBROUTINE READIN * 

**+:**+:********■& -k-k-kic-k-k-k-k 

This subroutine reads in the data. 



SUBROUTINE READIN 
REAL*4 DATA (100 0,3) 

INTEGER*4 NPTS 
INTEGER*4 PI 
INTEGERS P2 
CHARACTER’ 1 ' 1 LABEL 1 
CHARACTER*! LABEL2 
CHARACTER* 1 LABEL3 
CHARACTER* 1 LABEL4 
INTEGER*4 XTYPE 
INTEGER*4 YTYPE 
INTEGER*4 XQUAL 
INTEGER*4 YQUAL 
COMMON /DP ARMS/ DATA, NPTS, PI 
LABEL4 , XTYPE, 



80 

80 

80 

80 



, P2 , LABEL1, 
YTYPE, XQUAL 



LABEL2, 

YQUAL 



LABEL3 



Read the data 



DO 1010 1=1, NPTS 

READ (1,9010, END=8010) (DATA(I,J), J = 1, 3) 
1010 CONTINUE 
RETURN 



C+++ 

C+++ Abends 
C+++ 

8010 CALL ABEOF( 'READIN 1 , 1) 
9010 FORMAT (5X,3F10. 3) 

END 



C 

c 

c 

c 

c 

C+++ 

C+++ 



C+++ 



*********************** 



* SUBROUTINE RSCDAT * 
*********************** 



This subroutine re-scales the data 

SUBROUTINE RSCDAT (DATA, NPTS, EXPON, 
REAL*4 DATA (NPTS) 

INTEGERS NPTS 

INTEGER*4 EXPON 



MIN, MAX) 
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C+++ 

C+++ Generate scale factor and compute the data 
C+++ 

SF = 10**EXPON 
DO 1010 1=1, NPTS 

DATA(I) = DATA(I) * SF 
1010 CONTINUE 

MIN = MIN * SF 
MAX = MAX * SF 
RETURN 
END 



C 
C 

c 
c 
c 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



1010 



c 



cm 



kkkkkkkkkkkkkkkkkkkkkkk 

* SUBROUTINE SCLDAT * 

k k k k k k * k k k k k k k * * * * * * * k k 



This subroutine scales the data (if necessary) 

SUBROUTINE SCLDAT (DATA. NPTS, EXPON, MIN, MAX) 

REAL *4 DATA (NPTS) 

INTEGERS NPTS 

INTEGERS EXPON 

REAL*4 MIN 

REAL*4 MAX 

Get order of magnitude of minimum and maximum 

IF (MIN .NE. 0) THEN 

IOMMIN = LOG10 ( ABS (MIN ) ) 

ELSE 

IOMMIN = 0 
ENDIF 

IF (MAX .NE. 0) THEN 

IOMMAX = LOG10 (ABS (MAX) ) 

ELSE 

IOMMAX = 0 
ENDIF 

If the order of magnitude is within 10**(+/-3), then just set 
exponent to 0 and don't scale 

IF (((IOMMIN .GE. -3) .AND. (IOMMIN .LE. 3)) .AND. 

| ((IOMMAX .GE. -3) .AND. (IOMMAX .LE. 3))) THEN 
EXPON = 0 
ELSE 



Compute a suitable scaling factor 

IF (ABS (IOMMAX) .GT. 3) THEN 

IOMMAX = (IOMMAX / 3) * 3 
SF = 10 .** IOMMAX 
EXPON = IOMMAX 

ENDIF 

IF (ABS (IOMMIN) .GT. 3) THEN 

IOMMIN = (IOMMIN / 3) * 3 
SF = 10 . **IOMMIN 
EXPON = IOMMIN 

ENDIF 

DO 1010 1=1, NPTS 

DATA ( I ) = DATA ( I ) / SF 
CONTINUE 

MIN = MIN / SF 

MAX = MAX / SF 
ENDIF 
RETURN 
END 
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c 

c 

c 

C+++ 



C+++ 

C+++ 



* SUBROUTINE TIMS10 * 

’k'k'k'k'k’kjs'k'k'k'k'k'k'k'A’k'k’k'kjs'kjcjz 



This subroutine puts the '*10**' string in 

SUBROUTINE TIHS10(INVEC, I) 

CHARACTER* 11 INVEC 

INTEGERS I 

INVEC = ‘X 10 | EH. 75) 1 

I = I + 11 

RETURN 

END 



C 
C 
C 
C 
C 

C+++ 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



* -k * * * * * * * * * * * * * ’k * * * -k * * * 



* SUBROUTINE XTRACT * 

^i^^'k'k'k'kiK'k'k'k'k'k'k'k'k'k'k'k'k'k'k’k 



This subroutine scans for a token, checks it against the token 
table (word) and returns the index of the token in the table 

SUBROUTINE XTRACT (LINE, I, WORD, NWORDS , INDEX) 

CHARACTERS LINE (80) 

CHARACTERS WORD (NWORDS) 

CHARACTERS TOKEN 

Get the token 



DO 10 J=l, NWORDS 
PRINT *, WORD(J) 

10 CONTINUE 
IERR = I 

CALL GETTOK(LINE , I, TOKEN) 

C+++ 

C+++ Scan the table for the token. If it is not there, write an error . 
C+++ message and terminate. 

C+++ 



DO 1010 INDEX = 1, NWORDS 
IF (WORD (INDEX) .EQ. TOKEN) RETURN 
1010 CONTINUE 



C+++ 

C+++ Error message 
C+++ 

CALL ERRMSG(IERR, 'Unidentifiable command or token. $') 
CALL DISEND 
STOP 300 
END 
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3. GEOMETRIES PLOTTING CODE (GTDGEOMP) 
a. Flow Chart 



Main Program 



OKlirn 

wftai** 



Mti uo co—on 
block* 



Mtl UO 3 1 OCX 

au 




HOC : mx. No. of cylinder* 

M>X : MX. No. of Ditto* 

NS AX t MX. No. of *oux* *rr*y* 
NSX t mx. No. of xxrcn 




StflJL 



Pff*- 



cyl 1no*r at* 4 
coaoutM *nd 
ctP** loci 



5 U 3 R. pfflCPL 


J 


•na of p»oC 

» 


1 



stop 



I r«oi Tn 
I ^^eouorcy tn 
CHZ 




rnoi to dtu 



© 
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Flow Cnart For GTDGEOMP (continue) 
Mafn Program 




6 $ 



Flow Chart For GTDGEOMP (continue) 
Main Program 
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Flow Chart For GTDGEOMP (continue) 
Subroutine PLOT 



start 



suaa- 






cxtraaa 


«noa of 


• 11 00 


I«ct* « 


Q«t. af 


fil-a la 



SUlfl. 



vAwaj 
rtioiTn 
#vtng poln 
loci it on 

I 

StM. P 



1SS£I_ 



MU UC 15 

plot 



CflSIflY 

no I oc. 



SUIH. 
cat. trw loc. 
Va llna* Join 
tr» cyltnoar 
end paw 
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b. Program Listing 



C+++ ****************************** 

C+++ * GTDGEOMP * 

C+++ * * 

C+++ * WRITTEN BY KOH WEE JIN * 

C+++ * ON 9 SEP 1986 * 

C+++ * AT NPS * 

£+- 1 — {- ******* **************** ******* 

C+++ THIS PROGRAM READS IN THE DATA FROM GTD INPUT CARDS 
C+++ IT THEN PLOTS THE FOLLOWING GEOMETRIES 
C+++ CYLINDER(S) 

C+++ INFINITE GROUND PLANE 
C+++ PLATE (S) 

C+++ SOURCE (S) 

C+++ SOURCE ARRAY ( S ) 

C+++ THE PLOTS INCLUDES THE EFFECTS DUE TO 

C+++ AXES ROTATION 

C+++ ORIGIN TRANSLATION 

C+++ SCALING 

C+++ 

C+++ 



C+++ MAXIMUM NUMBER OF PLATES = NPX 

C+++ MAXIMUM NUMBER OF EDGES PER PLATES = NEX 

C+++ MAXIMUM NUMBER OF CYLINDERS = NCX 

C+++ MAXIMUM NUMBER OF POINTS FOR CYLINDER END CAPS = NCECP 

C+++ MAXIMUM NUMBER OF SINGLE SOURCES = NSX 

C+++ MAXIMUM NUMBER OF SOURCE ARRAYS = NSAX 

C+++ MAXIMUM NUMBER OF ELEMENTS PER SOURCE ARRAY = NESAX 



C+++ 

PARAMETER 
PARAMETER 
PARAMETER 
PARAMETER 
PARAMETER 
PARAMETER 
PARAMETER 

C+++ 

C+++ DECLARE VARIABLES 
C+++ 

C+++ PLATE 

C+++ CORNER LOCATIONS : XPC,YPC,ZPC 

C+++ NO. OF PLATES : IPG 

C+++ NO. OF CORNERS PER PLATE : MEP 

C+++ CYLINDER 

C+++ CENTER : XCCYL , YCCYL , ZCCYL 

C+++ END CAP CENTERS : ZCN.ZCP 

C+++ END CAP LOCI : XCYLP,YCYLP,ZCYLP,XCYLN,YCYLN,ZCYLN 

C+++ END CAP CUT ANGLES : THTN , THTP 

C+++ AXES ROTATION MATRIX : RCYL 

C+++ AXES LENGTH FROM CENTER : AC,BC 

C+++ NO. OF CYLINDER : ICG 

C+++ SINGLE SOURCE 

C+++ CORNER LOACTIONS : XSC,YSC,ZSC 

C+++ NO. OF SOURCES : ISG 

C+++ SOURCE ARRAY 

C+++ CORNER LOACTIONS : XSAC , YSAC , ZSAC 

C+++ NO. OF SOURCES s ISA 

C+++ NO. OF ELEMENTS PER SOURCE ARRAY s MSAX 

C+++ TOTAL NO. OF ELEMENTS : MSAXT 

C+++ INFINITE GROUND PLANE 

C+++ WITH GROUND PLANE : LGP = .TRUE. 

C+++ WITHOUT GROUND PLANE •. LGP = .FALSE. 

C+++ MIN. & MAX. X, Y & Z VALUES COMPARISON 

C+++ CYLINDER : COMPCG 

C+++ PLATE : COMPPG 

C+++ SOURCE : COMPSG 

C+++ SOURCE ARRAY : COMPSA 



i NPX=50) 
NEX=12) 
NCX=100) 
NCECP=201) 
NSX=1000 ) 
NSAX=5) 
NESAX=500) 
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C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



OBSERVATION POINT IN GLOBAL PLANE 
ANGLE FROM X-AXIS : PHI 
ANGLE FROM Z-AXIS s THETA 
RADIUS FROM ORIGIN : RADIUS 
GLOBAL AXES ROTATION MATRIX : ROT 
GLOBAL AXES ORIGIN TRANSLATION s XCR,YCR,ZCR 
FREQUENCY IN GHZ : FREQ 
SCALING FACTOR : SCALEF 
INTEGER DEFINED BY UN CARD : IUN 
INTEGER DEFINED BY US CARD : IUS 

REAL COMPPG(3 , 2 ) ,COMPCG(3 ,2) ,COMPSG(3 ,2) ,COMPSA(3,2) 

REAL ROT (3,3) , ROBS ( 3 , 3 ) 

REAL XPC(NPX,NEX) , YPC(NPX,NEX) ,ZPC(NPX,NEX) 

REAL XCYLP (NCX , NCECP ) , YCYLP (NCX , NCECP ) , ZCYLP (NCX , NCECP ) 

REAL XCYLN( NCX, NCECP) ,YCYLN(NCX, NCECP) ,ZCYLN (NCX, NCECP) 

REAL XCCAPP (NCX) ,YCCAPP(NCX) ,ZCCAPP(NCX) 

REAL XCCAPN(NCX) ,YCCAPN(NCX) ,ZCCAPN(NCX) 

REAL XSC(NSX,4) ,YSC(NSX,4) ,ZSC(NSX,4) 

REAL XSAC (NSAX , NESAX , 4 ) ,YSAC(NSAX,NESAX,4) ,ZSAC(NSAX,NESAX,4) 
INTEGER MEP(NPX) ,MSAX(NSAX) 

CHARACTER*60 TITLE, LABEL 
CHARACTERS C2 
LOGICAL LGP,LGEOM 

SET UP COMMON BLOCKS 



^ COMMON / CGDATA/ 

COMMON /OBDATA/ 
COMMON / OSDATA/ 
COMMON /PGDATA/ 
COMMON /PIDATA/ 
COMMON /RTDATA/ 
COMMON /SADATA/ 
COMMON /SCALER/ 
COMMON / SGDATA/ 
COMMON /VUDATA/ 
TITLE= 1 $ 1 
LGP= . FALSE . 
LGEOM=. FALSE. 
IST=0 



I CG , COMPCG , XCCAPP , YCCAPP , ZCCAPP , 
XCCAPN , YCCAPN , ZCCAPN 
XOBS, YOBS, ZOBS, ROBS 
POSMAX , TOSMAX , PHIOS , THETOS 
IPG , MEP , COMPPG 
PI , RTD 

XCR , YCR , ZCR , ROT 

ISA , MSAX , COMPS A , MSAXT 

FREQ, IUN, IUS, SCALEF 

ISG , COMPSG 

PHI, THETA, RAD IUS 



* 



INITIALIZE DISSPLA 



CALL DISSIN 



READ IN KEY LETTERS AND SELECT THE BRANCHES 



C+++ 

10 CONTINUE 



READ (1 


eg 


, IX, A60) 1 ) 


C2, LABEL 


IF < 


C2 


.EQ. 


'CT' 


i GOTO 


20 


IF < 


C2 


.EQ. 


'CG' 


» GOTO 


30 


IF ( 


C2 


.EQ. 


■EN' 


i GOTO 


40 


IF < 


C2 


.EQ. 


'FR' 


i GOTO 


50 


IF ( 


C2 


.EQ. 


'GP' 


i GOTO 


60 


IF < 


C2 


.EQ. 


'NC 


» GOTO 


70 


IF < 


C2 


.EQ. 


'NG' 


i GOTO 


80 


IF < 


C2 


.EQ. 


'NP' 


1 GOTO 


90 


IF < 


C2 


.EQ. 


'NS' 


i GOTO 


100 


IF < 


C2 


.EQ. 


'NX' 


> GOTO 


110 


IF < 


C2 


.EQ. 


'PG' 


» GOTO 


120 


IF 1 


C2 


.EQ. 


' RT ' J 


> GOTO 


130 


IF 1 


C2 


.EQ. 


' SA 1 ; 


I GOTO 


140 


IF I 


C2 


.EQ. 


'SG' 


1 GOTO 


150 


IF i 


C2 


.EQ. 


'UF' 


1 GOTO 


160 


IF i 


C2 


.EQ. 


'UN' 


> GOTO 


170 


IF 1 


[ C2 


.EQ. 


'US' 


> GOTO 


180 


IF 1 


I C2 


.EQ. 


• xq • ; 


) GOTO 


190 
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GOTO 10 

C+++ 

C+++ CT : COMMENT TITLE CARD 
C+++ 

20 CONTINUE 

WRITE (2,*) 

WRITE (2,*) 'CT' 

TITLE=LABEL 

WRITE (2,*) ' TITLE =' , TITLE 

GOTO 10 

C+++ 

C+++ CG : CYLINDER GEOMETRY CARD 
C+++ 

30 CONTINUE 

WRITE (2,*) 

WRITE (2,*) 'CG' 

LGEOM=.TRUE. 

ICG=ICG+1 

WRITE (2,*) ' ICG =' , ICG 

IF (ICG .GT. NCX) GOTO 35 

CALL CYLCAP ( XCYLP , YCYLP , ZCYLP , XCYLN , YCYLN , ZCYLN ) 
GOTO 10 

35 PRINT *, 'NO. OF CYLINDERS EXCEEDS NCX' 

CALL DONE PL 
STOP 

C+++ 

C+++ EN : END CARD 
C+++ 

40 CONTINUE 

WRITE (2 ,*) 

WRITE (2,*) 'EN' 

CALL DONEPL 
STOP 

C+++ 

C+++ FR : FREQUENCY CARD 
C+++ 

50 CONTINUE 

WRITE (2,*) 

WRITE (2,*) 'FR' 

LGEOM= . TRUE . 

READ (1,*) FREQ 

WRITE (2,*) ' FREQ = ' , FREQ 

GOTO 10 

C+++ 

C+++ GP s GROUND PLANE CARD 
C+++ 

60 CONTINUE 

WRITE (2,*) 

WRITE (2,*) 'GP' 

LGEOM= . TRUE . 

LGP= . TRUE . 

WRITE (2,*) ' LGP =' , LGP 

READ (1,*) LSLAB 
IF (LSLAB .EQ. 0) GOTO 10 
READ (1,*) DUMMY, DUMMY, DUMMY, DUMMY 
GOTO 10 

C+++ 

C+++ NC : NO CYLINDER CARD 
C+++ 

70 CONTINUE 

WRITE (2 , *) 

WRITE (2,*) 'NC' 

LGEOM= .TRUE . 

ICG=0 

CALL NOCYL ( XCYLP , YCYLP , ZCYLP , XCYLN , YCYLN , ZCYLN ) 
GOTO 10 

C+++ 

C+++ NG ; NO GROUND PLANE CARD 
C+++ 
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80 CONTINUE 

WRITE (2,*) 

WRITE (2 , * ) 'NG' 

LGEOM=.TRUE. 

LGP= . FALSE . 

GOTO 10 

C+++ 

C+++ NP : NO PLATE CARD 
C+++ 

90 CONTINUE 

WRITE (2,*) 

WRITE (2,*) 'NP' 

LGEOM=.TRUE. 

IPG=0 

CALL RESET 1 ( 1 , NPX , NEX , XPC , YPC , ZPC , COMPPG ) 

CALL RESET2(NPX,MEP , 1 , DUMMY, 1 , DUMMY) 

GOTO 10 

C+++ 

C+++ NS s NO SOURCE CARD 
C+++ 

100 CONTINUE 

WRITE (2,*) 

WRITE (2,*) 'NS' 

LGEOM= . TRUE . 

ISA=0 

ISG=0 

CALL RESET1 { 1 ,NSX , 4 , XSC , YSC , ZSC , COMPSG) 

CALL RESET1 (NSAX , NESAX , 4 , XSAC , YSAC , ZSAC , COMPSA) 

CALL RESET2 (NSAX, MSAX, 1, DUMMY, 1, DUMMY) 

GOTO 10 

C+++ 

C+++ NX : NEXT SET OF NEW DATA CARD 
C+++ 

110 CONTINUE 

WRITE (2 , *) 

WRITE (2,*) 'NX' 

LGEOM= . FALSE . 

CALL NXT ( XPC , YPC , ZPC , XCYLP , YCYLP , ZCYLP , XCYLN , YCYLN , ZCYLN , 
*XSC , YSC , ZSC , XSAC , YSAC , ZSAC , LGP , TITLE ) 

GOTO 10 

C+++ 

C+++ PG PLATE CARD 
C+++ 

120 CONTINUE 

WRITE (2,*) 

WRITE (2,*) 'PG' 

LGEOM=.TRUE. 

IPG=IPG+1 

WRITE (2,*) ' IPG =' , IPG 

IF (IPG .GT. NPX) GOTO 125 
CALL PLATEG ( XPC , YPC , ZPC ) 

GOTO 10 

125 PRINT *, 'NO. OF PLATES EXCEEDS NPX' 

CALL DONEPL 
STOP 

C+++ 

C+++ RT : ROTATE TRANSLATE AXES CARD 
C+++ 

130 CONTINUE 

WRITE (2,*) 

WRITE (2,*) 'RT' 

LGEOM=.TRUE. 

READ(1 , *) XCR , YCR , ZCR 

READ ( 1 , *) THZP , PHZP , THXP , PHXP 

CALL ROT ATE ( THZP , PHZP , THXP , PHXP , ROT ) 

GOTO 10 

C+++ 

C+++ SA : SOURCE ARRAY CARD 
C+++ 
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140 CONTINUE 

WRITE (2,*) 

WRITE (2,*) • SA ' 

LGEOM= . TRUE . 

ISA=ISA+1 

WRITE (2,*) 1 ISA =' , ISA 

IF (ISA .GT. NSAX) GOTO 145 
IST=ISA+ISG+MSAXT 
IF (1ST .GT. NSX) GOTO 147 
CALL SORCEA ( XSAC , YSAC , ZSAC ) 

GOTO 10 
145 CONTINUE 

PRINT *, 'THE TOTAL NO. OF SOURCE ARRAYS = 1 ,ISA 
PRINT *, 'NO. OF SOURCE ARRAYS EXCEEDS NSAX' 

CALL DONEPL 
STOP 

147 CONTINUE 

PRINT *, 'THE TOTAL NO. OF SOURCES = ',IST 
PRINT *, 'NO. OF SOURCES EXCEEDS NSX' 

CALL DONEPL 
STOP 

C+++ 

C+++ SG : SINGLE SOURCE CARD 
C+++ 

150 CONTINUE 

WRITE (2 , *) 

WRITE (2,*) 'SG' 

LGEOM= . TRUE . 

ISG=ISG+1 

WRITE (2,*) ' ISG =' , ISG 

IST=ISA+ISG+MSAXT 
IF (1ST .GT. NSX) GOTO 155 
CALL SORCEG ( XSC , YSC , ZSC ) 

GOTO 10 
155 CONTINUE 

PRINT *, 'THE TOTAL NO. OF SOURCES = ',IST 
PRINT *, 'NO. OF SOURCES EXCEEDS NSX' 

CALL DONEPL 
STOP 

C+++ 

C+++ UF : GEOMETRY SCALE FACTOR CARD 
C+++ 

160 CONTINUE 

WRITE (2,*) 

WRITE (2,*) 'UF' 

LGEOM= .TRUE . 

READ ( 1 , *) SCALEF 

WRITE (2,*) ' SCALEF =' , SCALEF 

GOTO 10 

C+++ 

C+++ UN : PLATES AND CYCLINDERS GEOMETRY SCALING CARD 
C+++ 

170 CONTINUE 

WRITE (2,*) 

WRITE (2,*) 'UN' 

LGEOM= . TRUE . 

READ ( 1 , * ) I UN 

WRITE (2,*) ' I UN =' , I UN 

GOTO 10 

C+++ 

C+++ US s SOURCES GEOMETRY SCALING CARD 
C+++ 

180 CONTINUE 

WRITE (2,*) 

WRITE (2,*) 'US' 

LGEOM= .TRUE . 

READ(1 ,*) IUS 

WRITE (2 ,*) ' IUS =' ,IUS 

GOTO 10 
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C+++ 

C+++ XQ s EXECUTION CARD 
C+++ 

190 CONTINUE 

WRITE (2,*) 

WRITE (2,*) 'XQ' 

IF (.NOT. LGEOM) GOTO 10 

CALL PLOT ( XPC , YPC , ZPC , XCYLP , YCYLP , ZCYLP , XCYLN , YCYLN , ZCYLN , 

*XSC , YSC , ZSC , XSAC , YSAC , ZSAC , LGP , TITLE ) 

LGEOM=. FALSE. 

GOTO 10 
END 

C 

C 

c*** 

C*** SUBROUTINE BLANK C 
C*** 

C+++ 

C+++ THIS SUBROUTINE DETERMINES THE PORTION OF ELLIPSE OF THE 

C+++ CYLINDER TO BE BLANKED OFF 

C+++ 

SUBROUTINE BLANKC ( ICYL , XCYL , YCYL , ZCYL , ICBDY , PN$ , LBLNK1 , LBLNK2 ) 
PARAMETER (NCX=100) 

PARAMETER (NCECP=201) 

REAL ROBS (3,3) , COMPCG (3,2) 

REAL XCYL(NCX, NCECP) ,YCYL(NCX, NCECP) ,ZCYL(NCX,NCECP) 

REAL X(NCECP) ,Y(NCECP) ,Z(NCECP) 

REAL XCCAPP(NCX) , YCCAPP(NCX) ,ZCCAPP(NCX) 

REAL XCCAPN(NCX) ,YCCAPN(NCX) ,ZCCAPN(NCX) 

CHARACTER PN$ 

LOGICAL LBLNK1 , LBLNK2 

COMMON / CGDATA/ ICG, COMPCG, XCCAPP , YCCAPP ,ZCCAPP , 

* XCCAPN , YCCAPN ,-ZCCAPN 

COMMON /OBDATA/ XOBS , YOBS ,ZOBS , ROBS 
LBLNK1= .TRUE . 

LBLNK2=.TRUE. 

WRITE (2,*) 

WRITE (2,*) 'SUBROUTINE BLANKC' 

IF (PN$ .EQ. ' P ' ) WRITE (2,*) ' POSITIVE END CAP' 

IF (PN$ .EQ. ' N ' > WRITE (2,*) ' NEGATIVE END CAP' 

DO 10 1=1 , NCECP 
X( I )=XCYL (ICYL , I ) 

Y ( I ) =YCYL ( I CYL , I ) 
z(l)=ZCYL(lCYL,I) 

10 CONTINUE 
C+++ 

C+++ TRANSFER THE LOCUS OF THE END CAP TO OBSERVATION CO-ORDINATE 
C+++ 

CALL ROTRAN (NCECP , X , Y , Z , XOBS , YOBS , ZOBS , ROBS ) 

C+++ 

C+++ TRANSFER THE CENTER POINTS OF THE TWO END CAPS TO OBSERVATION 

C+++ CO-ORDINATE 

C+++ 

XCP=XCC APP ( I CYL ) 

YCP=YCCAPP (ICYL ) 

ZCP=ZCCAPP (ICYL ) 

XCN=XCCAPN (ICYL ) 

YCN=YCCAPN ( ICYL ) 

zcn=zccapn(icyl) 

CALL ROTRAN ( 1 , XCP , YCP , ZCP , XOBS , YOBS , ZOBS , ROBS ) 

CALL ROTRAN ( 1 , XCN , YCN , ZCN , XOBS , YOBS , ZOBS , ROBS ) 

C+++ 

C+++ CONSIDER ONLY THE X-Y PLANE 

C+++ DETERMINE THE INTERCEPTION POINTS BETWEEN THE ELLIPSE AND 
C+++ THE CENTER LINE FORMED BY JOINING THE CENTERS OF THE TWO 
C+++ END CAPS 
C+++ 

NHALF= (NCECP- 1 )/2 
ISLOPE=l 
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IF (ABS (XCP-XCN) .LT. 0.00001) GOTO 15 
SLOPE= ( YCP- YCN ) / (XCP-XCN ) 

DO 20 I=ICBDY, ICBDY+NHALF 
11=1 

IF (II .GE. NCECP) I1=I1-NCECP+1 
12 = 11+1 

YI l=YCP+SLOPE* (X(Il)-XCP) 

YI 2=YCP+SLOPE* ( X ( I 2 ) - XCP ) 

IF (Y(I1) .GE. YI1 .AND. Y(I2) .LE. YI2) GOTO 30 
IF (Y(I1) .LE. YI1 .AND. Y(I2) .GE. YI2) GOTO 30 
20 CONTINUE 
GOTO 80 
15 CONTINUE 
ISLOPE=0 

DO 25 I=ICBDY, ICBDY+NHALF 
11=1 

IF (II .GE. NCECP) I1=I1-NCECP+1 
12 = 11+1 



IF (X(I1) 


.GE. XCP .AND. X< 


; i2 ) 


.LE. 


XCP) 


GOTO 


30 


IF (X(I1) 
25 CONTINUE 
GOTO 80 


.LE. XCP .AND. X< 


\12) 


.GE. 


XCP) 


GOTO 


30 



30 CONTINUE 
IBLNK1=I1 
IBLNK2=I1+NHALF 

IF (IBLNK2 .GT. NCECP) IBLNK2=IBLNK2 -NCECP 

C+++ 

C+++ DETERMINE THE DISTANCES FROM THE INTERCEPTION POINTS TO THE CENTER 

C+++ POINT OF THE OPPOSITE END CAP 

C+++ 

XC=XCP 

YC=YCP 

ZC=ZCP 

IF (PN$ .EQ. ' N 1 ) GOTO 40 
XC=XCN 
YC=YCN 
ZC=ZCN 
40 CONTINUE 

D1=SQRT((X(IBLNK1)-XC)*(X(IBLNK1)-XC)+ 
*(Y(IBLNK1)-YC)*(Y(IBLNK1)-YC) ) 

D2=SQRT ( ( X ( IBLNK2 ) -XC ) * ( X ( IBLNK2 ) -XC ) + 

* ( Y ( IBLNK2 ) - YC ) * ( Y ( IBLNK2 ) - YC ) ) 

C+++ 

C+++ COMPARE THE TWO DISTANCES AND INCULDING IN THE DEPTH (Z CO-ORD) 

C+++ TO DETERMINE IF THE PORTION IS BEING BLOCKED OR NOT 

C+++ 

IF (D1 .EQ. D2) GOTO 50 
IF (ISLOPE .EQ. 0) GOTO 55 
IF (D1 .GT. D2) GOTO 60 

ZCI=ZCN+(ZCP-ZCN)*(X(IBLNK1)-XCN)/ (XCP-XCN) 

IF (ABS(Z(IBLNK1) ) .GT. ABS(ZCI)) LBLNK1= .FALSE . 

GOTO 50 
60 CONTINUE 

ZCI=ZCN+ (ZCP-ZCN ) * (X ( IBLNK2 ) -XCN ) / (XCP-XCN ) 

IF (ABS (Z ( IBLNK2 ) ) .GT. ABS(ZCI)) LBLNK2= . FALSE . 

GOTO 50 
55 CONTINUE 

IF (D1 .GT. D2) GOTO 70 

ZCI=ZCN+ (ZCP-ZCN )*(Y(IBLNK1) -YCN)/ (YCP- YCN) 

IF (ABS (Z (IBLNK1 ) ) .GT. ABS(ZCI)) LBLNK1= . FALSE . 

GOTO 50 
70 CONTINUE 

ZCI=ZCN+ (ZCP-ZCN )*(Y( IBLNK2 ) -YCN ) / (YCP -YCN ) 

IF (ABS (Z ( IBLNK2 ) ) .GT. ABS(ZCI)) LBLNK2= . FALSE . 

50 CONTINUE 

WRITE (2,*) 'SUBROUTINE BLANKC END' 

RETURN 

80 PRINT *, 'THE INTERCEPTION POINT WAS NOT FOUND ' 

PRINT *, 'PROBABLY THE ALGORITHM IS NOT GOOD ENOUGH' 
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c — 

C*** 

C*** 

C*** 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



c— 

c*** 

c*** 

C*** 

C+++ 

C+++ 

C+++ 

C+++ 



PRINT *, 'TRY OTHER METHOD' 

CALL DONEPL 

STOP 

END 



BLOCK DATA 

BLOCK DATA 

SET UP PARAMETERS 

PARAMETER (NPX=50) 

PARAMETER (NCX=100) 

PARAMETER (NSX=1000) 

PARAMETER (NSAX=5) 

DECLARE VARIABLES AND ARRAYS 

REAL COMPPG (3,2), COMPCG (3,2), COMPSG (3,2), COMPSA (3,2) 
REAL ROT (3,3), ROBS (3,3) 

REAL XCCAPP(NCX) ,YCCAPP(NCX) ,ZCCAPP(NCX) 

REAL XCCAPN (NCX) ,YCCAPN(NCX) ,ZCCAPN(NCX) 

INTEGER MEP(NPX) ,MSAX(NSAX) 

SET UP COMMON BLOCKS 

COMMON / CGDATA/ ICG , COMPCG ,XCCAPP ,YCCAPP ,ZCCAPP , 

* XCCAPN, YCCAPN,ZCCAPN 

COMMON /OBDATA/ XOBS , YOBS , ZOBS , ROBS 
COMMON / OSDATA/ POSMAX, TOSMAX, PHIOS , THETOS 
COMMON /PGDATA/ IPG, MEP, COMPPG 
COMMON /PIDATA/ PI,RTD 
COMMON /RTDATA/ XCR,YCR,ZCR,ROT 
COMMON /SADATA/ ISA,MSAX, COMPSA, MSAXT 
COMMON /SCALER/ FREQ , IUN, IUS , SCALEF • 

COMMON /SGDATA/ I SG, COMPSG 
COMMON /VUDATA/ PHI , THETA, RADIUS 

INITIALIZE DATA SETS 

DATA COMPCG /3*10 . OEIO , 3*-10 . OEIO/ 

DATA COMPPG /3*10 . OEIO , 3*-10 . OEIO/ 

DATA COMPSA /3*10 . OEIO , 3*-10 .OEIO/ 

DATA COMPSG /3*10 . OEIO , 3*-10 . OEIO/ 

DATA XCR,YCR,ZCR / 3*0.0/ 

DATA ROT /l. 0,3*0. 0,1. 0,3*0. 0,1.0/ 

DATA FREQ, IUN, IUS, SCALEF /0 . 2997925 , 1 , 0 , 1 .0/ 

DATA PHI , THETA, RADIUS /30 . 0, 60 . 0 , 1 . 0/ 

DATA ICG, IPG, ISA, ISG, MSAXT /5*0/ 

DATA PI ,RTD /3 . 14159 , 57 . 29578/ 

DATA POSMAX, TOSMAX/4. 0,9.0/ 

END 



SUBROUTINE COMP 



THIS SUBROUTINE COMPARES THE OVERALL MIN. & MAX. X, Y & Z VALUES 
IT ALSO COMPUTES THE PLOTTING SIZE NEEDED BY DISSPLA 

SUBROUTINE COMP (C , Cl , C2 , C3 , C4 ,L) 

REAL C(3,2),C1(3,2),C2(3,2),C3(3,2),C4(3,2),L(3) 

WRITE (2,*) 

WRITE (2,*) 'SUBROUTINE COMP' 

DO 10 1=1.3 

IF (C(I,1) .GT. C1(I,1)) C(I,1)=C1(I,1) 

IF (C(I,1) .GT. C2 ( I , 1 ) ) C(I,1)=C2(I,1) 

IF (C(I,1) .GT. C3 ( I , 1 ) ) C( I , 1 )=C3 (1,1) 
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IF (C( I , 1 ) .GT. C4 ( I , 1 ) ) C(I , 1)=C4(I , 1) 
10 CONTINUE 

DO 20 1=1,3 

IF (C(l,2) .LT. Cl (I , 2) ) C(I , 2)=C1 (I ,2) 
IF (c(l,2) .LT. C2 (i , 2) ) c(l , 2)=C2 (1 , 2) 
IF (c(l,2) .LT. C3 (i , 2) ) c(l , 2)=C3 ( I , 2) 
IF (c(l,2) .LT. C4( I , 2) ) C(I , 2)=C4(I , 2) 
L(I)=C(I,2)-C(I,1) 

IF (L(I) .LT. 0.01) THEN 

CENTER= (C(I,2)+C(I,l))/2.0 
C(I , 2)=CENTER+0 . 005 
c( 1 , 1 )=CENTER-0 . 005 
L(I)=0.01 
END IF 

20 CONTINUE 

WRITE (2,*) 'SUBROUTINE COMP END' 

RETURN 

END 

C 

c*** 

C*** SUBROUTINE COMPAR 
r*** 

C+++ 



C+++ 

C+++ 



THIS SUBROUTINE COMPARES THE MIN. & MAX. 

SUBROUTINE COMPAR(N,C,X,Y,Z) 

REAL C(3,2),X(N),Y(N),Z(N) 



IF 
IF 
IF 
IF 
IF 
IF 

10 CONTINUE 
RETURN 
END 



(l,lj .GT. X(I) 




1 C 1,1 


>=X(I) 


(1,2) .LT. X(I 




1 c(l , 2 


l=X I 


(2,1) .GT. Y(l 




1 C(2,l 


l=Y I 


(2,2) .LT. Y(l 




1 C(2,2 


l=Y('I 


(3,1) .GT. Z(l 




» C 3,1 


l=Z I 


(3,2) .LT. Z(I 




1 C(3,2l 


l=z(l) 



C*** 

c*** 

C*** 



SUBROUTINE CYLBDY 



C+++ 



X, Y & Z VALUES 



C+++ THIS SUBROUTINE COMPUTES THE POSITION OF THE POINT ICBDY 
C+++ NEEDED TO DRAW THE TWO LINES JOINING THE ELLIPSES 
C+++ 

SUBROUTINE CYLBDY ( I CYL , XCYLP , YCYLP , ZCYLP , ICBDY) 
PARAMETER (NCX=100) 

PARAMETER (NCECP=201) 

REAL XCYLP ( NCX, NCE CP) ,YCYLP(NCX,NCECP) , ZCYLP (NCX,NCE CP) 
REAL ROBS (3,3), COMPCG (3,2) 

REAL XCCAPP (NCX) ,YCCAPP(NCX) ,ZCCAPP(NCX) 

REAL XCCAPN(NCX) ,YCCAPN(NCX) ,ZCCAPN(NCX) 

REAL X(NCECP) ,Y(NCECP) ,Z(NCECP) 

INTEGER ICBDY (NCX) 

COMMON / CGDATA/ ICG, COMPCG, XCCAPP, YCCAPP,ZCCAPP, 

* XCCAPN , Y CCAPN , ZCCAPN 

COMMON /OBDATA/ XOBS , YOBS ,ZOBS , ROBS 
COMMON /PI DATA/ PI,RTD 
COMMON /VUDATA/ PHI , THETA, RADIUS 
WRITE (2,*) 

WRITE (2,*) 'SUBROUTINE CYLBDY' 

WRITE (2,*) ' ICYL =',ICYL 

C+++ 



C+++ TRANSFER THE LOCUS OF THE POSITIVE END CAP TO 

C+++ OBSERVATION CO-ORDINATE 

C+++ 



DO 10 1=1 ,NCECP 
X ( I ) =XCYLP ( I CYL , I ) 
Y ( I )=YCYLP ( ICYL, I ) 
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C+++ 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



Z( I )=ZCYLP (ICYL , I ) 

10 CONTINUE 

CALL ROTRAN (NCECP , X , Y , Z , XOBS , YOBS , ZOBS , ROBS ) 

TRANSFER THE CENTERS OF THE TWO END CAPS TO 
OBSERVATION CO-ORDINATE 

XCP=XCCAPP (ICYL ] 

YCP=YCCAPP(ICYL 
ZCP=ZCCAPP ( ICYL 
XCN=XCCAPN(ICYL 
YCN=YCCAPN(ICYL 
ZCN=ZCCAPN ( I CYL ) 

CALL ROTRAN ( 1 , XCP , YCP , ZCP , XOBS , YOBS , ZOBS , ROBS ) 

CALL ROTRAN ( 1 , XCN , YCN , ZCN , XOBS , YOBS , ZOBS , ROBS ) 

THE PLANE FORM BY X(OBS)-AXIS AND Y(OBS)-AXIS IS THE PLANE 
THAT WILL BE DISPLAYED ON THE SCREEN. 

HENCE, WE CAN REMOVE THE Z(OBS) COMPONENTS AND WORK ONLY WITH 
X(OBS) AND Y(OBS) COMPONENTS. 

THE POINT THAT THE LINES (BODY OF THE CYLINDER) MEET THE ELLIPSES 
OF THE END CAPS WILL HAVE MAXIMUM PERPENDICULAR DISTANCE H 
FROM THE LINE JOINING THE CENTERS OF THE TWO END CAPS. 

ONLY ONE POINT (ICBDY) NEED TO BE COMPUTED BECAUSE THE OTHER 
POINTS ARE SYMMETRICAL TO IT AND CAN BE COMPUTED EASILY 

H1=0 
H2=0 

D=SQRT( (XCP-XCN)*(XCP-XCN)+(YCP-YCN)*(YCP-YCN) ) 

IF (D .LT. 0.00001) GOTO 50 
DO 40 1=1 , NCECP 

D1=SQRT((X(I)-XCP)*(X(I)-XCP)+(Y(I)-YCP)*(Y(I)-YCP)) 
D2=SQRT((X(I)-XCN)*(X(I)-XCN)+(Y(I)-YCN)*(Y(i)-YCN)) 

IF (D1 .LT. 0.000001) GOTO 20 
DUMMY=(Dl*Dl+D*D-D2*D2)/2 .0/D/D1 
IF (ABS (DUMMY ) .GT. 1.0001) GOTO 70 

IF (ABS (DUMMY) .GT. 1.0 .AND. ABS (DUMMY) .LE. 1.0001) THEN 
DUMMY=S IGN ( 1 . 0 , DUMMY ) 

END IF 

ALPHA=ACOS (DUMMY) 

H3=D1*SIN (ALPHA) 

IF (I .LE. 2) GOTO 30 

IF (H2 .GE. HI .AND. H2 .GE. H3) GOTO 60 
GOTO 30 
CONTINUE 
H3=0. 000001 
CONTINUE 
H1=H2 
H2=H3 
CONTINUE 
GOTO 60 
CONTINUE 
ICBDY ( ICYL)=1 
CONTINUE 
ICBDY (I CYL) =1-1 

WRITE (2 , *) 'SUBROUTINE CYLBDY END' 

RETURN 



20 

30 

40 

50 

60 

70 



PRINT * 


, 1 ARGUMENT 


FOR ACR 


COS EXCEEDS 


LIMIT OF +/- 1.0 


WRITE < 


f 2 ,*) ' 


ARGUMENT 


FOR ACR COS 


EXCEEDS +/- 1.0' 


WRITE ( 


2,* 1 


D 


= ' ,D 


WRITE ( 


2 ,*) 1 


D1 


= ' ,D1 




WRITE 1 


\ 2 ' l ' 


D2 


= ' , D2 




WRITE 1 


2,*) 1 


ARG 


=' , DUMMY 





CALL 

STOP 

END 



DONEPL 
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SUBROUTINE CYLCAP 



C*** 

r*** 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



THIS SUBROUTINE COMPUTES THE POINTS NEEDED TO DRAW THE ELLIPSES 



ANGLES NEEDED TO DEFINE THE CYLINDER CO-ORDINATE 



THETA TO ZT AXIS 
PHI TO ZT AXIS 
THETA TO XT AXIS 
PHI TO XT AXIS 



TCLZ 

PCLZ 

TCLX 

PCLX 



SUBROUTINE CYLCAP (XCYLP , YCYLP , ZCYLP , XCYLN , YCYLN , ZCYLN) 
PARAMETER (NCX=100) 

PARAMETER (NCECP=201) 

REAL COMPCG(3 ,2) ,ROT(3 , 3 ) ,R(3 , 3) 

REAL XCYLP (NCX,NCE CP) , YCYLP (NCX,NCECP) , ZCYLP (NCX,NCECP) 
REAL XCYLN (NCX,NCECP) , YCYLN(NCX,NCECP) , ZCYLN (NCX,NCECP) 
REAL XCCAPP(NCX) , YCCAPP(NCX) ,ZCCAPP(NCX) 

REAL XCCAPN(NCX) ,YCCAPN(NCX) ,ZCCAPN(NCX) 

REAL XP(NCECP) ,YP(NCECP) ,ZP(NCECP) 

REAL XN(NCECP) ,YN(NCECP) ,ZN(NCECP) 

COMMON /CGDATA/ ICG,COMPCG,XCCAPP,YCCAPP,ZCCAPP, 

* XCCAPN , YCCAPN , ZCCAPN 

COMMON /PIDATA/ PI,RTD 
COMMON /RTDATA/ XCR, YCR,ZCR,ROT 
COMMON /SCALER/ FREQ , I UN , IUS , SCALEF 
WRITE (2 , *) 

WRITE (2 
WRITE (2 , : 

READ 
READ 
READ 
READ 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 



'SUBROUTINE CYLCAP 1 
' ICG =' , ICG 



XCCYL , YCCYL , ZCCYL 
TCLZ, PCLZ, TCLX, PCLX 
AC , BC 

ZCN , THTN , ZCP , THTP 
' XCCYL = ' , XCCYL 

' YCCYL =', YCCYL 

' ZCCYL =', ZCCYL 

1 TCLZ =', TCLZ, ' 

' TCLX =', TCLX, ' 

' AC =' , AC , ' BC =' , BC 

1 ZCN =' , ZCN, 1 ZCP =' ,ZCP 

' THTN = 1 , THTN , ' THTP =',THTP 



PCLZ =' , PCLZ 
PCLX = l *,PCLX 



COMPUTE THE LOCI OF THE ELLIPSES IN CYLINDER CO-ORDINATE 



NHALF= (NCECP+1 ) /2 
RTHP=THTP/RTD 
RTHN=THTN / RTD 
XSTEP=2 . 0*AC/ (NHALF- 1 ) 

XSTART=-1 . 0*AC-XSTEP 
DO 10 1=1, NHALF 
XP ( I ) =XSTART+REAL ( I ) *XSTEP 
YP(I )=BC*SQRT( 1 . 0- (XP( I )/AC)**2 ) 
IF (THTP .EQ. 90.0) THEN 
ZP(I)=ZCP 
ELSE 

ZP ( I ) =ZCP+XP ( I ) /TAN ( RTHP ) 

END IF 
XN ( I ) =XP ( I ) 

YN ( I ) =YP ( I ) 

IF (THTN .EQ. 90.0) THEN 
ZN ( I ) =ZCN 
ELSE 

ZN ( I ) =ZCN+XN ( I ) /TAN ( RTHN ) 

END IF 

XP(NCECP-I+1)=XP(I) 

YP (NCECP-I+1 )=-l . 0^YP( I ) 
ZP(NCECP-I+1)=ZP(I) 

XN (NCECP-I+1 )=XN( I) 
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YN(NCECP-I+1)=-1.0*YN(I) 

ZN (NCECP- 1+1 )=ZN(I) 

10 CONTINUE 

YP(NCECP)=YP(1) 

YN(NCECP)=YN(1) 

YP (NHALF ) =- 1 . 0*YP (NHALF ) 

YN (NHALF ) =- 1 . 0*YN (NHALF ) 

C+++ 

C+++ COMPUTE THE AXES ROTATION MATRIX FOR CYLINDER CO-ORDINATE 
C+++ 

CALL ROTATE (TCLZ , PCLZ , TCLX , PCLX , R ) 

C+++ 

C+++ TRANSFER THE POINTS TO GLOBAL CO-ORDINATE 



C+++ 

CALL ROTRAN (NCECP , XP , YP , ZP , XCCYL , YCCYL , ZCCYL , R) 

CALL ROTRAN (NCECP , XN , YN , ZN , XCCYL , YCCYL , ZCCYL , R) 

CALL ROTRAN (NCECP , XP , YP , ZP , XCR , YCR , Z CR , ROT ) 

CALL ROTRAN (NCECP , XN , YN , ZN , XCR , YCR , ZCR , ROT ) 

XCCAPP? I CG)=0 
YCCAPP? I CG)=0 
ZCCAPP? I CG)=ZCP 

xccapn(icg)=o 

yccapn?icg)=o 

zccapn?icg)=zcn 

CALL ROTRAN ( 1 , XCCAPP (ICG), YCCAPP (ICG), ZCCAPP (ICG), 

A XCCYL YCCYL ZCCYL R) 

CALL ROTRAN ( 1 , XCCAPN (ICG) ,YCCAPN(ICG) ,ZCCAPN( ICG) , 

* XCCYL, YCCYL, ZCCYL, R) 

CALL ROTRAN (1, XCCAPP (ICG) ,YCCAPP(ICG) ,ZCCAPP(ICG) , XCR, YCR, ZCR, ROT 
CALL ROTRAN (l, XCCAPN ( ICG) ,YCCAPN?ICG) ,ZCCAPN(lCG) , XCR, YCR, ZCR, ROT 

C+++ 

C+++ SCALE THE POINTS 



C+++ 



CALL 

CALL 

CALL 

CALL 



SCALES (NCECP , 1 G 1 ,XP,YP,ZP) 
SCALES(NCECP, ' G ’ , XN , YN , ZN ) 
SCALES ( 1 , ' G ' , XCCAPP (ICG) , YCCAPP 
SCALES (l, ' G * , XCCAPN? ICG) ,YCCAPN 



1=1, NCECP 
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XCYLP(ICG,I] 


l=XP< 


;n 


YCYLP(lCG,I 


l=YP< 


;i) 


ZCYLP( ICG , I 


l=ZP< 


i) 


xcyln(icg,i 


l=XN( 


i) 


ycyln(icg,i 


!=YN< 


i) 


zcyln(icg,i; 


l=ZNl 


;i) 


CONTINUE 








ZCCAPP (ICG)) 

zccapn(icg)) 



C+++ 



C+++ COMPUTE THE MIN. & MAX. X, Y & Z VALUES FOR CYLINDER GEOMETRY 
C+++ 

CALL COMP AR (NCECP , COMPCG , XP , YP , ZP ) 

CALL COMP AR( NCECP, COMPCG, XN,YN,ZN) 

WRITE (2 , *) 'SUBROUTINE CYLCAP END 1 

RETURN 

END 



C 

c*** 

c *** SUBROUTINE DISSET 
C*** 



C+++ 



C+++ THIS SUBROUTINE SETS UP THE GRAPHIC DIMENSIONS, AXES AND TITLE 
C+++ 

SUBROUTINE D I S S ET ( COMPAL , LENGTH , TITLE) 

REAL COMPAL (3,2), LENGTH ( 3 ) 

CHARACTER*? PHI $ ,THETA$ , RADI $ 

CHARACTER*60 TITLE , LABEL1 , LABEL2 , LABEL3 
COMMON /VUDATA/ PHI , THETA, RAD I US 
WRITE (2,*) 

WRITE (2,*) 'SUBROUTINE DISSET' 

OPEN (3, STATUS= ' NEW ' ) 

WRITE (3,100) PHI, THETA, RADIUS 



82 



100 FORMAT (IK, F7 . 2 , IX, F7 . 2 , IX, F7 . 2 ) 

CLOSE (3) 

OPEN (3, STATUS= 1 OLD 1 ) 

READ (3,200) PHI$ , THETAS ,RADI$ 

200 F0RMAT(1X,A7 , IX, A7 , IX, A7) 

LABEL1=' VIEWING POINT : ?F) = '//PHI$//' DEG$ 1 
LABEL2= 1 ?Q) = 1 / /THETA$// 1 DEG$ 1 

LABEL3= 1 ( )R = '//RADI$//' M $' 

CLOSE (3) 

CALL AREA2D(8. 0,8.0) 

CALL VOLM3D ( LENGTH ( 1 ) , LENGTH ( 2 ) , LENGTH ( 3 ) ) 

CALL HEAD IN (TITLE ,100,1.0,4) 

CALL HEADINi LABEL1 , 100 ,1.0,4) 

CALL HEAD IN (LABEL2 , 100 ,1.0,4) 

CALL HEAD IN (LABEL3 ,100,1.0,4) 

CALL X3NAME ( 1 ( X ) IN M$',100) 

CALL Y3NAME( ' (Y) INM$',100) 

CALL Z3NAME( 1 (Z) INM$',100) 

THETA V=90. 0-THETA 

CALL VUANGL( PHI, THETAV. RADIUS) 

CALL GRAF3D ( COMPAL (1,1),' SCALE 1 , COMPAL (1,2), 

*COMPAL (2,1),' SCALE 1 , COMPAL(2 , 2 ) , COMPAL (3 , 1 ) , 'SCALE' , C0MPAL(3 , 2 ) ) 
WRITE (2,*) 'SUBROUTINE DISSET END 1 
RETURN 
END 

C 

c*** * 

C*** SUBROUTINE DISSIN 

C*** 

C+++ 

C+++ THIS SUBROUTINE INITIALIZE DISSPLA 
C+++ 

SUBROUTINE DISSIN 

CALL COMPRS 

CALL PAGE (8.5,11.0) 

CALL FUTURA 

CALL BASALF ( ' L/ CSTD 1 ) 

CALL MIXALF ( 1 STANDARD 1 ) 

CALL MX3ALF ( 1 L/CGREEK 1 , 1 ? 1 ) 

CALL NOBRDR 
CALL YAXANG(O.O) 

RETURN 

END 

C 

C*** SUBROUTINE INVERS 
c*** 

C+++ 

C+++ THIS SUBROUTINE INVERSES THE 3 X 3 MATRIX RIN 
C+++ 

SUBROUTINE INVERS (RIN, ROUT) 

REAL RIN (3,3), ROUT (3,3) 

DETERM=RIN (1,1) *RIN ( 2 , 2) *RIN ( 3 , 3) +RIN { 1 , 2) *RIN ( 2 , 3 ) *RIN (3,1)+ 

* RIN(1,3)*RIN(2,1)*RIN(3,2)-RIN(1,1)*RIN(2,3)*RIN(3,2)- 

* RIN(1,2)*RIN(2,1)*RIN(3,3)-RIN(1,3)*RIN(2,2)*RIN(3,1) 

IF (DETERM .EQ. 0.0) GOTO 10 

ROUT (1 , 1 )=(RIN(2 , 2)*RIN(3 , 3) -RIN(2 , 3 )*RIN(3 , 2) ) /DETERM 
ROUT (1 , 2 )=(RIN( 1 ,3 )*RIN(3 , 2 ) -RIN(1 , 2)*RIN(3 , 3 ) ) /DETERM 
ROUT(l ,3)=(RIN(1, 2 )*RIN( 2 , 3 ) -RIN( 1 , 3 )*RIN(2 , 2 ) ) /DETERM 
ROUT (2 , 1 )=(RIN(2 , 3)*RIN(3 , 1 ) -RIN(2 , 1 )*RIN(3 , 3 ) ) /DETERM 
ROUT (2 , 2 )=(RIN( 1 ,1)*RIN(3,3)-RIN(1, 3 )*RIN(3 , 1 ) ) /DETERM 
ROUT(2 , 3)=(RIN(1 ,3)*RIN(2,1)-RIN(1, 1 )*RIN(2 , 3) ) /DETERM 
ROUT (3 ,1)=(RIN(2, 1 )*RIN(3 , 2 ) -RIN( 2 , 2 )*RIN(3 , 1 ) ) /DETERM 
ROUT (3 , 2 )=(RIN ( 1 , 2)*RIN(3 , 1) -RIN ( 1 , 1 )*RIN(3 , 2 ) ) /DETERM 
ROUT (3 ,3) = (RIN(1, 1 ) 5<: RIN(2 , 2) -RIN ( 1 , 2 )*RIN( 2 , 1 ) )/ DETERM 
RETURN 
10 CONTINUE 

PRINT *, 'THE DETERMINANT OF RIN IS EQUAL TO ZERO' 

WRITE (2,*) 'SUBROUTINE INVERS 1 
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WRITE (2,*) 

WRITE (2,*) 

WRITE (2,*) 

WRITE (2,*) 

WRITE (2,*) 

CALL DONEPL 
STOP 
END 

C 

c*** 

C*** SUBROUTINE NOCYL 
r*** 

C+++ 

C+++ THIS SUBROUTINE RESETS ALL DATA BELONG TO CYLINDER GEOMETRY 
C+++ 

SUBROUTINE NOCYL ( XCYLP , YCYLP , ZCYLP , XCYLN , YCYLN , ZCYLN ) 

PARAMETER (NCX=100) 

PARAMETER (NCECP=20l) 

REAL COMPCG (3,2) , R ( 3 , 3 ) 

REAL XCYLP (NCX , NCECP ) , YCYLP ( NCX , NCECP ) , ZCYLP (NCX,NCECP) 

REAL XCYLN(NCX, NCECP) ,YCYLN(NCX , NCECP) , ZCYLN (NCX, NCECP) 

REAL XCCAPP (NCX) , YCCAPP (NCX) , ZCCAPP (NCX) 

REAL XCCAPN (NCX) ,YCCAPN(NCX) ,ZCCAPN(NCX) 

COMMON / CGDATA/ ICG, COMPCG, XCCAPP, YCCAPP, ZCCAPP, 

* XCCAPN, YCCAPN,ZCCAPN 

CALL RESET1 ( 1 , NCX , NCECP , XCYLP , Y.CYLP , ZCYLP , COMPCG ) 

CALL RESET1 ( 1 , NCX , NCECP , XCYLN , YCYLN , ZCYLN , COMPCG) 

CALL RESET2 (NCX , XCCAPP , NCX , YCCAPP , NCX , ZCCAPP ) 

CALL RESET2 (NCX , XCCAPN , NCX , YCCAPN , NCX , ZCCAPN ) 

RETURN 

END 

C - - 

c *** 

C*** SUBROUTINE NXT 

C*** * 

C+++ 

C+++ THIS SUBROUTINE RESET ALL GEOMETRY DATA 
C+++ 

SUBROUTINE NXT ( XPC ,YPC ,ZPC , XCYLP , YCYLP , ZCYLP , XCYLN, YCYLN, ZCYLN, 
*XSC , YSC , ZSC , XSAC , YSAC , ZSAC , LGP , TITLE ) 

PARAMETER (NPX=50) 

PARAMETER (NEX=12) 

PARAMETER (NCX=100) 

PARAMETER (NCECP=201) 

PARAMETER (NSX=1000) 

PARAMETER (NSAX=5) 

PARAMETER (NESAX=500) 

REAL COMPPG (3,2) , COMPCG(3 , 2) , COMPSG(3 , 2) , COMPSA(3 , 2) 

REAL ROT (3 , 3) ,ROBS(3.3) 

REAL XCYLP (NCX, NCECP) , YCYLP (NCX, NCECP) , ZCYLP (NCX, NCECP) 

REAL XCYLN (NCX, NCECP) , YCYLN (NCX, NCECP) ,ZCYLN(NCX, NCECP) 

REAL XCCAPP (NCX) ,YCCAPP(NCX) , ZCCAPP (NCX) 

REAL XCCAPN(NCX) .YCCAPN(NCX) , ZCCAPN(NCX) 

REAL XPC(NPX,NEX) , YPC(NPX.NEX) ,ZPC(NPX,NEX) 

REAL XSC(NSX, 4) ,YSC(NSX,4) ,ZSC(NSX,4) 

REAL XSAC (NSAX , NESAX , 4 ) ,YSAC(NSAX,NESAX,4) ,ZSAC(NSAX,NESAX,4) 
INTEGER MEP(NPX) ,MSAX(NSAX) 

CHARACTERS TITLE 
LOGICAL LGP 

COMMON /CGDATA/ ICG, COMPCG, XCCAPP, YCCAPP, ZCCAPP, 

* XCCAPN, YCCAPN, ZCCAPN 
COMMON /OBDATA/ XOBS , YOBS , ZOBS , ROBS 
COMMON /PGDATA/ IPG ,MEP , COMPPG 
COMMON /RTDATA/ XCR,YCR,ZCR,ROT 
COMMON /SADATA/ ISA,MSAX,COMPSA,MSAXT 
COMMON /SCALER/ FREQ , IUN , IUS , SCALEF 
COMMON /SGDATA/ ISG,COMPSG 

COMMON /VUDATA/ PHI , THETA , RADIUS 

CALL RESET1 ( 1 , NPX , NEX , XPC , YPC , ZPC , COMPPG ) 



THE DETERMINANT OF RIN IS EQUAL TO ZERO' 
THE MATRIX RIN HAS THE FOLLOWING VALUES' 
' ,RIN(1 , 1 ) , ' ' ,RIN(1 ,2) , ' ' ,RIN(l ,3) 

' , RIN(2 , 1 ) , ' ' , RIN(2 , 2 ) , ' ',RIN(2,3) 

' ,RIN(3 , 1) , ' ' , RIN(3 , 2) , ' ',RIN(3,3) 
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CALL RESET1 ( 1 ,NCX ,NCECP , XCYLP , YCYLP , ZCYLP , COMPCG) 

CALL RESET 1 ( 1 , NCX , NCE CP , XCYLN , YCYLN , ZCYLN , COMPCG ) 

CALL RESET1 ( 1 , NSX , 4 , XSC , YSC , ZSC , COMPSG) 

CALL RESET1 (NSAX , NESAX , 4 , XSAC , YSAC , ZSAC , COMPSA ) 

CALL RESET2 ( NPX , MEP , NSAX , MS AX . 1 , FREQ ) 

CALL RESET2 ( 1 ,XCR, 1 , YCR, 1 ,ZCR) 

CALL RESET2 ( 1 ,XOBS , 1 , YOBS , 1 ,ZOBS ) 

CALL RESET2 (NCX , XCCAPP , NCX , YCCAPP , NCX , ZCCAPP ) 

CALL RESET2 (NCX , XCCAPN , NCX , YCCAPN , NCX , ZCCAPN) 

CALL ROTATE (0.0, 0.0 ,90. 0,0.0, ROT) 

CALL ROTATE (0.0, 0.0, 90. 0,0.0, ROBS) 

LGP= . FALSE . 

TITLE= 1 $ 1 

SCALEF=1 . 0 

PHI=60 . 0 

THETA=30 . 0 

RADIUS=1 . 0 

ICG=0 

IPG=0 

ISA=0 

ISG=0 

IUN=1 

IUS=1 

RETURN 

END 

C 

C*** 

C*** SUBROUTINE OBSERV 
C*** 

C+++ 

C+++ THIS SUBROUTINE COMPUTES THE AXES ROTATION MATRIX AND 
C+++ ORIGIN TRANSLATION CO-ORDINATES FO.R THE CONVERSION FROM 
C+++ GLOBAL CO-ORDINATE TO OBSERVATION CO-ORDINATE 
C+++ 

C+++ IT ASSUMES THAT THE CENTER OF OBSERVATION IS AT THE CENTER 

C+++ OF THE WORK BOX DEFINED BY DISSPLA 

C+++ 

SUBROUTINE OBSERV (COMPAL) 

• REAL COMPAL (3,2), ROBS (3,3),R(3,3) 

CHARACTER ANSS 
CHARACTERS THETA$ ,PHI$ 

CHARACTERSO MSSG1,MSSG2 

COMMON /OBDATA/ XOBS , YOBS ,ZOBS , ROBS 

COMMON / OSDATA/ POSMAX , TOSMAX , PHIOS , THETOS 

COMMON /PI DATA/ PI,RTD 

COMMON /VUDATA/ PHI , THETA, RADIUS 

WRITE (2,*) 

WRITE (2,*) 'SUBROUTINE OBSERV' 

C+++ 

C+++ THE OFFSET ANGLES ARE NEEDED TO COMPENSATE FOR THE DIFFERENCES 
C+++ BETWEEN THE VIEWING ANGLES INTERPRETED BY THE PROGRAM AND THAT 
C+++ SEEN BY DISSPLA 
C+++ 

THETOS=REAL(INT(-l . 0*TOSMAX*SIN(2 .0*THETA/RTD) ) ) 

PHIOS=REAL ( INT ( - 1 . 0*POSMAX*SIN ( THETA/RTD ) ) ) 

C+++ 

C+++ DETERMINE THE CENTRE OF THE WORK BOX 
C+++ 

XCWB= ( COMPAL (1,2) +COMPAL ( 1 , 1 ) ) / 2 . 0 
YCWB= ( COMPAL (2,2) +COMPAL ( 2 , 1 ) ) / 2 . 0 
ZCWB= ( COMPAL (3,2) +COMPAL ( 3 , 1 ) ) / 2 . 0 

C+++ 

C+++ COMPUTE THE OBSERVATION POINT WITH RESPECT TO THE CENTER 

C+++ OF THE WORK BOX 

C+++ 

RPHI= (PHI+PHIOS ) /RTD 
RTHETA= ( THETA+THETOS ) /RTD 
XCOBS=RADIUS*SIN(RTHETA) *COS (RPHI ) 

YCOBS=RADIUS*S IN ( RTHETA ) *S IN ( RPHI ) 
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ZCOBS=RADIUS*COS (RTHETA) 

C+++ 

C+++ DETERMINE THE OBSERVATION POINT IN GLOBAL CO-ORDINATE 
C+++ 

XOBS=XCWB+XCOBS 

YOBS=YCWB+YCOBS 

ZOBS=ZCWB+ZCOBS 

C+++ 

C+++ DETERMINE THE NEW AXES ANGLES CENTERED AT THE OBSERVATION 
C+++ POINT WITH RESPECT TO THE GLOBAL CO-ORDINATE WITH THE 
C+++ Z(OBS) AXIS POINTING TOWARDS THE CENTER OF THE WORK BOX 
C+++ 

THETAZ=180 . 0- (THETA+THETOS ) 

PHIZ=180 .0+(PHI+PHI0S ) 

IF (PHIZ .GT. 360.0) PHIZ=PHIZ-360 .0 
THETAX=ABS ( 9 0 . 0 -THETAZ ) 

IF (THETAZ .LT. 90.0) PHIX=180 . 0+PHIZ 
IF (THETAZ .GE. 90.0) PHIX=PHIZ 
IF (PHIX .GT. 360.0) PHIX=PHIX-360 .0 

C+++ 

C+++ FORM AXES ROTATION MATRIX 
C+++ 

CALL ROTATE (THETAZ . PHIZ , THETAX , PHIX , R) 

CALL INVERS(R,ROBS) 

WRITE (2,*) 1 THE MATRIX ROBS HAS THE FOLLOWING VALUES' 

WRITE (2,*) ' ' ,ROBS(l,l) , ' ' , ROBS ( 1 , 2 ) , ' ',ROBS(l,3). 

WRITE (2,*) ' 1 , ROBS (2 , 1 ) , ' ' , ROBS (2 , 2 ) , ' ',ROBS(2, 3) 

WRITE (2,*) 1 ' ,ROBS(3,l) , ' ' ,ROBS(3,2) , ' ',ROBS(3,3) 

C+++ 

C+++ COMPUTE THE ORIGIN TRANSLATION 
C+++ 

XOBS=-1.0*XOBS 
YOBS=-1.0*YOBS 
ZOBS=-l . 0*ZOBS 

CALL ROTRAN ( 1 , XOBS , YOBS , ZOBS ,0.0,0. 0,0.0, ROBS ) 

WRITE (2,*) 1 THE ORIGIN TRANSLATION HAS THE FOLLOWING VALUES' 

WRITE (2,*) ' XOBS =',XOBS 

WRITE (2,*) ' YOBS =',YOBS 

WRITE (2,*) ' ZOBS =',ZOBS 

WRITE (2,*) 'SUBROUTINE OBSERV END' 

RETURN 

END 

C - 

c*** 

C*** SUBROUTINE PLATEG 
c*** 

C+++ 

C+++ THIS SUBROUTINE READS IN THE PLATE CORNERS AND TRANSFER THEM TO 

C+++ GLOBAL CO-ORDINATE 

C+++ 

SUBROUTINE PLATEG (XPC,YPC,ZPC) 

PARAMETER (NPX=50) 

PARAMETER (NEX=12) 

REAL XPC(NPX,NEX) ,YPC(NPX,NEX) ,ZPC(NPX,NEX) ,X(NEX) ,Y(NEX) ,Z(NEX) 
REAL COMPPG(3,2) ,R(3,3) ,R0T(3,3) 

INTEGER MEP(NPX) 

COMMON /PGDATA/ IPG ,MEP , COMPPG 
COMMON /RTDATA/ XCR,YCR,ZCR,ROT 
COMMON /SCALER/ FREQ , IUN, IUS , SCALEF 
WRITE (2,*) 

WRITE (2 ,*) 'SUBROUTINE PLATEG' 

C+++ 

C+++ READ IN DATA 
C+++ 

READ(1,*) MEP(IPG) ,LSLAB 

WRITE (2,*) ' MEP(IPG) =',MEP(IPG) 

IF (MEP(IPG) .GT. NEX) GOTO 50 
IF (LSLAB .EQ. 0) GOTO 20 
READ ( 1 , *) NS LAB 
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DO 10 I=1,NSLAB 

READ(1 , *) DUMMY, DUMMY, DUMMY, DUMMY, DUMMY 
10 CONTINUE 
20 CONTINUE 

DO 30 1=1 ,MEP(IPG) 

READ(1,*)X(I),Y(I),2(I) 

30 CONTINUE 
C+++ 

C+++ TRANSFER THE POINTS TO GLOBAL CO-ORDINATE 
C+++ 

CALL ROTRAN(MEP( IPG) ,X, Y,2,XCR,YCR,2CR,R0T) 

C+++ 

C+++ SCALE THE POINTS 
C+++ 

CALL SCALES (MEP( IPG) , 1 G 1 ,X,Y,2) 

CALL COMPAR(MEP( IPG) ,C0MPPG,X,Y,2) 

WRITE (2,*) 1 IPG =' ,IPG 

DO 40 1=1 .MEP(IPG) 

XPC (IPG , I )=X( I ) 

YPC(IPG,I)=Y(I) 

2PC(IPG , I )=2( I ) 

40 CONTINUE 

WRITE (2,*) 'SUBROUTINE PLATEG END' 

RETURN 

50 PRINT *, 'NO. OF CORNERS EXCEEDS NEX' 

CALL DONEPL 
'STOP 
END 



C 

C*** 

Qkkk 

rkkk 

C+++ 
C+++ 
* C+++ 



SUBROUTINE PLOT 



THIS SUBROUTINE PLOTS THE GEOMETRY OF THE INPUT DATA 

SUBROUTINE PLOT ( XPC , YPC , 2PC , XCYLP , YCYLP , 2CYLP , XCYLN , YCYLN , 2CYLN , 
*XSC , YSC , 2SC , XSAC , YSAC , 2SAC , LGP , TITLE ) 

PARAMETER (NPX=50) 

PARAMETER (NEX=12) 

PARAMETER (NCX=100) 

PARAMETER (NCECP=201) 

PARAMETER (NSX=1000) 

PARAMETER (NSAX=5) 

PARAMETER (NESAX=500) 

REAL XPC(NPX,NEX) , YPC(NPX,NEX) ,2PC(NPX,NEX) 

REAL XCYLP (NCX , NCECP ) , YCYLP (NCX,NCECP) ,2CYLP(NCX,NCECP) 

REAL XCYLN(NCX, NCECP) ,YCYLN(NCX, NCECP) ,2CYLN(NCX, NCECP) 

REAL XCCAPP (NCX) ,YCCAPP(NCX) ,2CCAPP(NCX) 

REAL XCCAPN(NCX) ,YCCAPN(NCX) ,2CCAPN(NCX) 

REAL XSC(NSX,4) ,YSC(NSX,4) ,ZSC(NSX,4) 

REAL XSAC (NS AX , NESAX , 4 ) , YSAC(NSAX,NESAX.4) ,2SAC(NSAX,NESAX,4) 
REAL COMPAL(3,2) , COMPCG(2 , 3) , COMPPG(3 , 2 ) , COMPSA(3 , 2 ) , COMPSG(3 , 2) 
REAL ROT (3,3). ROBS (3,3) 

REAL LENGTH(3) 

INTEGER MEP(NPX) ,MSAX(NSAX) 

CHARACTER*60 TITLE 
CHARACTER ANS$ 

LOGICAL LGP 

COMMON /CGDATA/ ICG , COMPCG , XCCAPP , YCCAPP , 2CCAPP , 

* XCCAPN , YCCAPN , 2CCAPN 

COMMON /OBDATA/ XOBS , YOBS , 20BS , ROBS 
COMMON /OSDATA/ POSMAX , TOSMAX , PHIOS , THETOS 
COMMON /PGDATA/ IPG,MEP,COMPPG 
COMMON /PI DATA/ PI,RTD 
COMMON /RTDATA/ XCR,YCR,2CR,R0T 
COMMON /SCALER/ FREQ , I UN , IUS , SCALEF 
COMMON /SGDATA/ ISG,COMPSG 
COMMON /SADATA/ ISA, MS AX, COMPSA,MSAXT 
COMMON /VUDATA/ PHI , THETA, RAD IUS 
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DATA COMPAL/3*10 . 0E10 , 3*-10 . OEIO/ 

WRITE (2,*) 

WRITE (2,*) 'SUBROUTINE PLOT' 

C+++ 

C+++ DETERMINE OVERALL MIN. & MAX. X, Y & Z VALUES 
C+++ 

CALL COMP ( COMPAL , COMPCG , COMPPG , COMPSA , COMPSG , LENGTH ) 

10 CONTINUE 
C+++ 

C+++ DETERMINE THE VIEWING ANGLES 
C+++ 

CALL VANGLE (TITLE, LENGTH ) 

C+++ 

C+++ SET UP DISSPLA FOR 3-D PLOTTING 
C+++ 

CALL DISSET( COMPAL, LENGTH, TITLE) 

C+++ 

C+++ COMPUTE ROTATION MATRIX AND TRANSLATION POSITION FOR 

C+++ OBSERVATION CO-ORDINATE 

C+++ 

CALL OBSERV( COMPAL) 

C+++ 

C+++ PLOT INFINITE GROUND PLANE ? 

C+++ 

IF (.NOT. LGP) GOTO 20 
CALL PLOTGP( COMPAL) 

20 CONTINUE 
C+++ 

C+++ PLOT PLATE (S) ? 

C+++ 

IF (IPG .EQ. 0) GOTO 30 
CALL PLOTPG ( XPC , YPC , ZPC ) 

30 CONTINUE 
C+++ 

C+++ PLOT CYLINDER(S) ? 

C+++ 

IF (ICG .EQ. 0) GOTO 40 

CALL PLOTCG ( XCYLP , YCYLP , ZCYLP , XCYLN , YCYLN , ZCYLN) 

40 CONTINUE 
C+++ 

C+++ PLOT SOURCE (S) ? 

C+++ 

IF ( ISG .EQ. 0) GOTO 50 
CALL PLOTSG(XSC,YSC,ZSC) 

50 CONTINUE 
C+++ 

C+++ PLOT SOURCE ARRAY (S) ? 

C+++ 

IF (ISA .EQ. 0) GOTO 60 
CALL PLOTSA(XSAC,YSAC ,ZSAC) 

60 CONTINUE 

CALL ENDPL(O) 

70 CONTINUE 

PRINT *, 'DO YOU WANT ANOTHER VIEW OF THE PLOT (Y OR N) ' 
READ ( 5 , ' (Al) ' ) ANS$ 

IF (ANS$ .NE. ' Y ' .AND. ANS$ .NE. ' N ' ) GOTO 70 
IF (ANS$ .EQ. ' N ' ) GOTO 80 
GOTO 10 
80 CONTINUE 

WRITE (2,*) 'SUBROUTINE PLOT END' 

RETURN 

END 

C 

C*** 

C *** SUBROUTINE PLOTCG 

i r-*** 

C+++ 

C+++ THIS SUBROUTINE PLOTS THE CYLINDER(S) 

C+++ 
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SUBROUTINE PLOTCG ( XCYLP , YCYLP , ZCYLP , XCYLN , YCYLN , ZCYLN ) 
PARAMETER (NCX=100) 

PARAMETER (NCECP=201) 

REAL ROT (3 ,3) ,ROBS(3,3) ,COMPCG(3,2) 

REAL XCYLP (NCX , NCECP ) , YCYLP (NCX , NCECP ) , ZCYLP (NCX,NCECP) 

REAL XCYLN ( NCX , NCECP ) . YCYLN (NCX , NCECP ) , ZCYLN (NCX , NCECP ) 

REAL X2(2) ,Y2(2) ,Z2(2) 

REAL XCCAPP(NCX) ,YCCAPP(NCX) ,ZCCAPP(NCX) 

REAL XCCAPN(nCX) ,YCCAPN(NCX) ,ZCCAPN(NCX) 

INTEGER ICBDY (NCX) 

LOGICAL LBLNK1 ,LBLNK2 

COMMON /CGDATA/ ICG,COMPCG,XCCAPP,YCCAPP,ZCCAPP, 

* XCCAPN , YCCAPN , ZCCAPN 

COMMON /OBDATA/ XOBS , YOBS ,ZOBS , ROBS 
COMMON /PIDATA/ PI,RTD 
COMMON /RTDATA/ XCR,YCR,ZCR,ROT 
COMMON /SCALER/ FREQ , IUN, IUS , SCALEF 
COMMON /VUDATA/ PHI , THETA, RADIUS 
WRITE (2,*) 

WRITE (2,*) 'SUBROUTINE PLOTCG' 

DO 10 1=1, ICG 

C+++ 

C+++ PLOT THE LINES JOINING THE ELLIPSES 
C+++ 

CALL CYLBDY ( I , XCYLP , YCYLP , ZCYLP , I CBDY ) 

X2(1)=XCYLP(I , ICBDY (I ) ) 

Y2(1)=YCYLP(I,ICBDY(I)) 

Z2 (1 )=ZCYLP( I , ICBDY( I ) ) 

X2 ( 2 ) =XCYLN (1,1 CBDY ( I ) ) 

Y2(2 }=YCYLN ( I , ICBDYiI ) ) 

Z2(2)=ZCYLN(I,ICBDY(I)) 

CALL CURV3D ( X2 . Y2 , Z2 , 2 , 0 ) 

NHALF= (NCECP- l)/2 

X2 ( 1 ) =XCYLP ( I , ICBDY ( I ) +NHALF ) 

Y2 ( 1 ) =YCYLP ( I , ICBDY ( I ) +NHALF ) 

Z2 ( 1 ) =ZCYLP ( I , ICBDY ( I ) +NHALF ) 

X2 ( 2 ) =XCYLN ( I , ICBDY ( I ) +NHALF ) 

Y2 ( 2 ) = YCYLN ( I , ICBDY ( I ) +NHALF ) 

Z2 ( 2 ) =ZCYLN ( I , ICBDY ( I ) +NHALF ) 

CALL CURV3D(X2 , Y2 , Z2 , 2 , 0) 

C+++ 

C+++ PLOT THE POSITIVE END CAP ELLIPSE 
C+++ BLANK OFF THE REGION THAT IS BLOCKED 
C+++ 

NHALF= (NCECP- l)/2 

CALL BLANKC(I , XCYLP, YCYLP, ZCYLP, ICBDY( I) , ' P ' , LBLNK1 , LBLNK2 ) 
CALL PLOTEP ( I , XCYLP , YCYLP , ZCYLP , ICBDY ( I ) , LBLNK1 , LBLNK2 ) 

C+++ 

C+++ PLOT THE NEGATIVE END CAP ELLIPSE 
C+++ BLANK OFF THE REGION THAT IS BLOCKED 
C+++ 

CALL BLANKC(I , XCYLN, YCYLN, ZCYLN, ICBDY ( I ) , ' N ' , LBLNK1 , LBLNK2 ) 
CALL PLOTEP ( I , XCYLN , YCYLN , ZCYLN , I CBDY ( I ) , LBLNK1 , LBLNK2 ) 

10 CONTINUE 

WRITE (2,*) 'SUBROUTINE PLOTCG END' 

RETURN 

END 

C 

C*** 

C*** SUBROUTINE PLOTEP 

c*** 

C+++ 

C+++ THIS SUBROUTINE PLOTS THE END CAP ELLIPSE OF CYLINDER 
C+++ 

SUBROUTINE PLOTEP ( ICYL , XCYL , YCYL , ZCYL , ICBDY , LBLNK1 , LBLNK2 ) 
PARAMETER (NCX=100) 

PARAMETER (NCECP=201) 

REAL XCYL (NCX, NCECP) , YCYL (NCX, NCECP) , ZCYL (NCX, NCECP) 

REAL X(NCECP) ,Y(NCECP) ,Z(NCECP) 



89 



LOGICAL LBLNK1 ,LBLNK2 
WRITE (2,*) 

WRITE (2,*) 'SUBROUTINE PLOTEP' 

WRITE (2,*) 1 LBLNK1=' ,LBLNK1 

WRITE (2,*) ' LBLNK2=' ,LBLNK2 

NHALF= (NCECP-1 ) /2 
IF (LBLNK1 ) GOTO 20 

C+++ 

C+++ PLOT THE ONLY FIRST HALF OF THE ELLIPSE 
C+++ 

DO 10 J=1 .NHALF+l 

WRITE (2,*) ' DO LOOP 1' 

K=ICBDY+NHALF+J-1 

IF (K .GT. NCECP) K=K-NCECP+1 

X ( J ) =XCYL ( I CYL , K } 

Y ( J ) =Y CYL ( I CYL , K ) 

2( J)=2CYL( ICYL ,K) 

WRITE (2,*) J,' ',X(J),' 1 , Y( J) , ' 1 ,2(J) 

10 CONTINUE 

CALL CURV3D (X,Y ,2, NHALF+l ,0) 

GOTO 60 
20 CONTINUE 

IF (LBLNK2) GOTO 40 

C+++ 

C+++ PLOT ONLY THE SECOND HALF OF THE ELLIPSE 
C+++ 

DO 30 J=l, NHALF+l 

WRITE (2,*) 1 DO LOOP 2' 

K=ICBDY+J-1 

IF (K .GT. NCECP) K=K-NCECP+1 
X ( J ) =XCYL ( I CYL , K ) 

Y ( J ) =YC YL ( I CYL , K ) 

2 ( J ) =2CYL ( I CYL , K ) 

WRITE (2 ,*) J,' ',X(J),' 1 , Y( J) , 1 ' ,2(J) 

30 CONTINUE 

CALL CURV3D(X ; Y,2, NHALF+l ,0) 

GOTO 60 
40 CONTINUE 
C+++ 

C+++ PLOT THE FULL ELLIPSE 
C+++ 

DO 50 J=1 .NCECP 

WRITE (2,*) 1 DO LOOP 3' 

X ( J ) =XCYL ( I CYL , J ) 

Y ( J ) =YCYL ( I CYL , J ) 

2 ( J ) =2CYL ( I CYL , J ) 

WRITE (2,*) J,' 1 ,X( J) , 1 1 / Y( J) , 1 1 ,2(J) 

50 CONTINUE 

CALL CURV3D ( X , Y , 2 , NCECP , 0 ) 

60 CONTINUE 

WRITE ( 2 ,*) 'SUBROUTINE PLOTEP END' 

RETURN 

END 

C 

C*** SUBROUTINE PLOTGP 
C*** 

C+++ 

C+++ THIS SUBROUTINE PLOTS THE INFINITE GROUND PLANE 
C+++ 

SUBROUTINE PLOTGP ( COMP AL) 

REAL COMPAL(3,2) ,X(5) ,Y(5) ,2(5) 

WRITE (2,*) 

WRITE (2 ,*) 'SUBROUTINE PLOTGP' 



x(i; 


l=COMPAL< 


l,i 


1*1.2 


X(2 


l=COMPAL ( 


l,i 


1*1.2 


X(3 


l=COMPAL< 


1,2 


1*1.2 


X(4 


t=COMPAL ( 


1,2 


1*1.2 


X(5 


»=X(1) 
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Y ( 1 )=C0MPAL(2 , 2)*1 . 2 
Y(2)=COMPAL(2,l)*1.2 
Y(3 )=C0MPAL(2 , 1 )*1 . 2 
Y(4)=C0MPAL(2,2)*1.2 
Y(5)=Y(1) 

Z(1)=0.0 
Z(2)=0.0 
Z(3)=0.0 
Z(4 =0.0 
Z(5}=0.0 

WRITE (2,*) 1 CORNERS OF INFINITE GROUND PLANE (GLOBAL)' 

WRITE (2,*) 1 I',' X(I)',' Y(I ) 1 , 1 Z(I) 

DO 10 1=1,5 

WRITE (2,*) I,X(I),Y(I),Z(I) 

10 CONTINUE 

CALL CURV3D(X, Y,Z , 5,0) 

WRITE (2,*) 'SUBROUTINE PLOTGP END' 

RETURN 

END 



C 

c*** 

c*** SUBROUTINE PLOTPG 

c*** 



C+++ 

C+++ 

C+++ 



THIS SUBROUTINE PLOTS THE PLATE (S) 

SUBROUTINE PLOTPG(XPC,YPC,ZPC) 

PARAMETER (NPX=50) 

PARAMETER (NEX=12) 

REAL XPC(NPX.NEX) ,YPC(NPX,NEX) ,ZPC(NPX,NEX) 
REAL X(NEX+1) ,Y(NEX+1) ,Z(NEX+1) 

INTEGER MEP(NPX) 

COMMON /PGDATA/ IPG,MEP , COMPPG 
WRITE (2,*) 

(2,*) 'SUBROUTINE PLOTPG' 



WRITE 
DO 30 1=1 .IPG 



WRITE 

WRITE 

WRITE 



PLATE NO . = ' , I 
CORNERS OF PLATE (GLOBAL) ' 
IV X(DV 



1 

2 ,*) 1 
2 *) 1 

DO"iO J=1 ,MEP(I ) 

X(J)=XPC(I/J) 

Y(J)=YPC(I,J) 

z(j)=zpc(i,j) 

10 CONTINUE 
M=MEP(I)+1 
X(M)=X( 1) 

y(m)=y(i) 

Z(M)=Z( 1 ) 

DO 20 K=1.M 

WRITE (2,*) K,X(K) , Y(K) ,Z(K) 

20 CONTINUE 

CALL CURV3D ( X , Y , Z , M , 0 ) 

30 CONTINUE 

WRITE (2,*) 'SUBROUTINE PLOTPG END' 

RETURN 

END 



Y(I)', 



Z(I) 



C 

C*** SUBROUTINE PLOTSA 
C+++ 

C+++ THIS SUBROUTINE PLOTS THE SOURCE ARRAY ( S ) 

C+++ 

SUBROUTINE PLOTSA ( XSAC ,YS AC , ZSAC) 

PARAMETER (NSAX=5) 

PARAMETER (NESAX=500) 

REAL XSAC (NS AX , NES AX , 4 ) ,YSAC(NSAX,NESAX,4) ,ZSAC(NSAX,NESAX,4) 
REAL X(5) ,Y(5) ,Z(5) 

INTEGER MSAX(NESAX) 



I 



I 
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COMMON /SADATA/ ISA,MSAX,COMPSA,MSAXT 
WRITE (2,*) 

WRITE (2,*) 'SUBROUTINE PLOTSA' 

DO 40 1=1, ISA 

WRITE (2,*) ' SOURCE ARRAY NO. =',I 

DO 30 J=1 ,MSAX(I ) 

WRITE (2,*) ' SOURCE ELEMENT NO. =',J 

WRITE (2,*) 1 CORNERS OF SOURCE ELEMENT (GLOBAL)' 

DO 10 K=1 ,4 
X(K)=XSAC(I, J,K) 

Y(K)=YSAC(I , J,K) 

Z(K)=ZSAC(l, J,K) 

10 CONTINUE 
X(5)=X(1) 

Y(5)=Y(1) 

Z(5)=Z(l) 

WRITE (2,*) ' I',' X(I) ' , ' Y(I) ' , 1 Z(I) 

DO 20 K=1 , 5 

WRITE (2,*) K,X(K) ,Y(K) ,Z(K) 

20 CONTINUE 

CALL CURV3D ( X , Y , Z , 5 , 0 ) 

30 CONTINUE 
40 CONTINUE 

WRITE (2,*) 'SUBROUTINE PLOTSA END' 

RETURN 

END 

C 

C*** SUBROUTINE PLOTSG 
C*** 

C+++ 

C+++ THIS SUBROUTINE PLOTS THE SINGLE SOURCE(S) 

C+++ 

SUBROUTINE PLOTSG(XSC,YSC,ZSC) 

PARAMETER (NSX=1000) 

REAL XSC(NSX,4) ,YSC(NSX,4) ,ZSC(NSX,4) 

REAL X(5) ,Y( 5) ,Z( 5) 

COMMON /SGDATA/ ISG,COMPSG 
WRITE (2,*) 

WRITE (2,*) 'SUBROUTINE PLOTSG' 

DO 30 1=1. ISG 

WRITE (2,*) ' SOURCE NO. =',I 

WRITE (2,*) ' CORNERS OF SOURCE (GLOBAL)' 

WRITE (2,*) ' I',' X( I ) ' , ' Y(I) ' , 1 Z(I) 

DO 10 J=1 ,4 

X(J)=XSC(I,J) 

Y( J )=YSC ( I , J) 

Z(j)=ZSC(l , J) 

10 CONTINUE 
X(5)=X(1) 

Y(5)=Y(1) 

Z(5)=Z(1) 

DO 20 K=1 , 5 

WRITE (2,*) K,X(K) ,Y(K) ,Z(K) 

20 CONTINUE 

CALL CURV3D(X,Y, 2,5,0) 

30 CONTINUE 

WRITE (2,*) 'SUBROUTINE PLOTSG END 1 

RETURN 

END 

C 

c*** 

C*** SUBROUTINE RESET1 

Q'k'k'k 

C+++ 

C+++ THIS SUBROUTINE INITIALIZES THE GEOMETRICAL CO-ORDINATES 

C+++ AS WELL AS THE COMPARISON VARIABLE 

C+++ 

SUBROUTINE RESET1 (II , 12 , 13 ,X, Y,Z,C) 
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10 

20 

30 



40 

50 



C 

C*** 

Qkkk 

rkkk 

C+++ 

C+++ 

C+++ 



10 

20 

30 



C 

Qkkk 

Qkkk 

Qkkk 

C++ + 
C+++ 
C+++ 



C 

C 

c 

c 

c 

c 

c 



REAL X(I1,I2,I3),Y(I1,I2,I3),Z(I1,I2,I3),C(3,2) 
DO 30 J1=1,I1 
DO 20 J2=l , 12 
DO 10 J3=l . 13 
X(J1,J2,J3)=0.0 
Y(J1, J2,J3)=0.0 
Z(J1,J2,J3)=0.0 
CONTINUE 
CONTINUE 
CONTINUE 
A=10 . 0E10 
DO 50 1=1,2 
DO 40 J=1 , 3 
C ( J , I ) =A 
CONTINUE 
A=-10 . 0E10 
CONTINUE 
RETURN 
END 



SUBROUTINE RESET2 



THIS SUBROUTINE RESETS 3 SINGLE DIMENSIONED ARRAYS TO ZERO 

SUBROUTINE RESET2 (II ,X1 , 12 ,X2 , 13 ,X3) 

REAL X1(I1) ,X2(I2) ,X3(I3) 

DO 10 J=1,I1 
X1(J)=0.0 
CONTINUE 
DO 20 J=1 , 12 
X2 ( J)=0 .0 
CONTINUE 
DO 30 J=1 , 13 
X3 ( J)=0 . 0 
CONTINUE 
RETURN 
END 



SUBROUTINE ROTATE 



THIS SUBROUTINE COMPUTES THE AXES ROTATION MATRIX 

SUBROUTINE ROTATE (THETAZ , PHIZ , THETAX , PHIX , R) 

REAL R ( 3 3 ) 

COMMON /PIDATA/ PI,RTD 

THZ=THETAZ/RTD 

PHZ=PHIZ/RTD 

THX=THETAX/RTD 

PHX=PHIX/RTD 

R(l,l)=COS(PHX)*SIN(THX) 

R(2,1)=SIN(PHX)*SIN(THX) 

R(3,l)=COS(THX) 

R(l,3)=COS(PHZ)*SIN(THZ) 

R(2,3)=SIN(PHZ)*SIN(THZ) 

R(3,3)=COS(THZ) 

D1=R( 1 ,3)*R(1,1)+R(2.3)*R(2,1)+R(3,3)*R(3,1) 

IF (ABS(Dl) .GT. 0.1) GOTO 10 
R(1,1)=R(1,1)-R(1,3)*D1 
R(2,1)=R(2,1)-R(2,3)*D1 
R(3 1 ) =R ( 3 i)-R(3 3)*D1 

D2=SQRT(R(i,l)*R(i,l)+R(2,l)*R(2 / l)+R(3 > l)*R(3,l)) 
R(1 , 1)=R( 1 , 1 )/D2 
R(2,1)=R(2,1)/D2 
R(3 , 1 )=R(3 , 1 )/D2 

R(1/2)=R(2,3)*R(3,1)-R(3,3)*R(2,1) 
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10 



R(2,2j=R(l,l)*R(3,3 



R(3;2)=R(1,3)*R(2 
RETURN 






R(l,3>*Rf3,l) 

R(2,3)*R(l,l) 



PRINT 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 



2 * 

2 '? 
2 * 

2 * 
2 ,* 
2 * 
2 ,*< 



■AXES SPECIFIED ARE NOT ORTHOGONAL 1 



CALL DONEPL 

STOP 

END 



SUBROUTINE ROTATE' 

AXES SPECIFIED ARE NOT ORTHOGONAL' 
THE INPUT PARAMETERS ARE AS FOLLOWS' 
THETA-Z =' ,THETAZ, ' DEG' 

PHI-Z = 1 ,PHIZ , ' DEG' 

THETA-X =' , THETAX, ' DEG' 

PHI-X = 1 ,PHIX, ' DEG' 



C 

c*** 

c *** SUBROUTINE ROTRAN 
C*** 



C+++ 

C+++ THIS SUBROUTINE ROTATES AND TRANSLATE CO-ORDINATES 



C+++ 

SUBROUTINE ROTRAN(N,X,Y,Z.XC,YC,ZC,R) 
REAL X(N) ,Y(N) ,Z(N) , R(3 , 3 ) 

DO 10 1=1, N 
X1=X(I) 

Y1=Y(I) 

Z1=Z(I) 

X(I)=R(1,1)*X1+R(1,2)*Y1+R(1,3)*Z1+XC 
Y(I)=R(2,1)*X1+R(2,2)*Y1+R(2,3)*Z1+YC 
Z(I)=R(3,1)*X1+R(3,2)*Y1+R(3 ; 3)*Z1+ZC 
10 CONTINUE 
RETURN 
END 



C 

c*** 

C*** SUBROUTINE SCALES 

C*** 



C+++ 

C+++ THIS SUBROUTINE SCALES THE GEOMETRIES 

C+++ ALPHABET 'S' FOR SOURCES 

C+++ AND ' G ' FOR ALL OTHER GEOMETRIES 



C+++ 



SUBROUTINE SCALES(N,A$,X,Y,Z) 
REAL X(N) ,Y(N) ,Z(N) 

CHARACTER A$ 



COMMON 
IF (A$ 
IF 
IF 
IF 

END IF 
IF (A$ 



/SCALER/FREQ , I UN , IUS , SCALEF 



.EQ. 
’I UN 
IUN 
IUN 
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)=SCALE*X( 

>=SCALE*Y( 

)=SCALE*Z( 



CONTINUE 

RETURN 

END 



THEN 

1) SCALE=SCALEF 

2) SCALE=0 . 3048*SCALEF 

3) SCALE=0.0254*SCALEF 

THEN 

0) SCALE=0.2997925/FREQ*SCALEF 
SCALE=SCALEF 
SCALE=0 . 3048*SCALEF 
SCALE=0 • 0254*SCALEF 



C 

c*** 

C*** SUBROUTINE SORCEA 

C*** 



C+++ 
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C+++ THIS SUBROUTINE READS IN THE DATA FOR SOURCE ARRAY AND COMPUTES 
C+++ THE CORNERS OF EACH ELEMENTS. 

C+++ IT ALSO COMPARE THE MIN. & MAX. X, Y & Z VALUES 
C+++ 



C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 



CENTER OF EACH ELEMENT s XCSA,YCSA,ZCSA 
ANGLES NEEDED TO DEFINE THE SOURCE ARRAY CO-ORDINATE 
THETA TO ZT AXIS : THSAZ 
PHI TO ZT AXIS : PHSAZ 
THETA TO XT AXIS : THSAX 
PHI TO XT AXIS : PHSAX 

AXES ROTATION MATRIX FOR SOURCE ARRAY CO-ORDINATE : RSA 



SUBROUTINE SORCEA ( XSAC , YSAC , ZSAC ) 

PARAMETER (NSX=1000) 

PARAMETER (NSAX=5) 

PARAMETER (NESAX=500) 

REAL XSAC(NSAX,NESAX,4) ,YSAC(NSAX,NESAX.4) ,ZSAC(NSAX,NESAX,4) 
REAL XCSA(NESAX) ,YCSA (NESAX) , ZCSA(NESAX) 

REAL X(4) ,Y(4) ,Z(4) 

REAL COMPSA(3 , 2) ,R(3 , 3 ) ,RSA(3 , 3) , COMPSG(3 , 2) 

INTEGER MSAX(NSAX) 

COMMON /RTDATA/ XCR,YCR,ZCR,R 
COMMON /SCALER/ FREQ , IUN, IUS , SCALEF 
COMMON /SADATA/ ISA,MSAX,COMPSA,MSAXT 
COMMON /SGDATA/ ISG,COMPSG 
WRITE " 



WRITE 



f!:*S 



'SUBROUTINE SORCEA' 



READ IN DATA 



C+++ 

WRITE (2 ,*) ' ISA =' , ISA 

READ ( 1 , * ) MSAX(ISA) 

WRITE (2,*) ' MSAX(ISA) =',MSAX(ISA) 

IF (MSAX(ISA) .GT. NESAX) GOTO 40 
MSAXT=MSAXT+MSAX( ISA) 

IST=ISA+ISG+MSAXT 

WRITE (2, *) 1 ISG = 1 , ISG 

WRITE (2,*) ' MSAXT =' < MSAXT 

WRITE ( 2 ,*) ' 1ST = 1 , 1ST 

IF (1ST .GT. NSX) GOTO 50 
DO 10 1=1 ,MSAX(ISA) 

READ(1 , *) XCSA(I) ,YCSA(I) ,ZCSA(I) 

10 CONTINUE 

CALL SCALES (MSAX( ISA) , ' G • ,XCSA,YCSA,ZCSA) 

READ ( 1 , * ) THSAZ , PHSAZ , THSAX , PHSAX 
WRITE (2,*) 1 THSAZ= 1 , THSAZ , ' PHSAZ= 1 , PHSAZ 

WRITE (2,*) ' THSAX=' , THSAX , 1 PHSAX=' , PHSAX 

READ ( 1 , * ) DUMMY, HS , HAWS 

WRITE (2 *) 1 HS = 1 ,HS , 1 HAWS =',HAWS 

READ ( 1 , * ) DUMMY, DUMMY 

C+++ 



C+++ COMPUTE THE AXES ROTATION MATRIX FOR SOURCE ARRAY 
C+++ 

CALL ROTATE (THSAZ , PHSAZ , THSAX , PHSAX , RSA ) 

C+++ 

C+++ DETERMINE THE CORNERS OF EACH SOURCE ELEMENTS IN SOURCE ARRAY 

C+++ CO-ORDINATE 

C+++ 

DO 30 1=1 ,MSAX(ISA) 

X(l)=HAWS/2.0 
X ( 2 ) =HAWS / 2 . 0 
X(3 )=-l . 0*HAWS/2 . 0 
X(4)=-l . 0*HAWS/2 . 0 
Y(1)=0.0 
Y 2=0.0 
Y(3 )=0 . 0 
Y(4)=0 . 0 
Z(l)=HS/2.0 
Z(2)=-l . 0*HS/2 . 0 
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Z(3)=-1.0*HS/2.0 

Z(4)=HS/2.0 

CALL SCALES (4, 'S' , X,Y,Z) 

C+++ 

C+++ TRANSFER THE CORNER LOCATIONS TO GLOBAL CO-ORDINATE 
C+++ 

CALL R0TRAN(4,K,Y / Z,KCSA(I) ,YCSA(I) ,ZCSA(I) ,RSA) 

CALL ROTRAN ( 4 , X , Y , Z , XCR , YCR , ZCR , R) 

DO 20 J=1 , 4 

XS AC (ISA, I , J)=X( J) 

YSAC(ISA,I,J)=Y(J) 

ZSAC(ISA, I, J)=Z( J) 

20 CONTINUE 
C+++ 

C+++ COMPARE THE MIN. & MAX. X, Y & Z VALUES 
C+++ 

CALL C0MPAR(4, COMPSA,X,Y,Z) 

30 CONTINUE 

WRITE (2,*) 'SUBROUTINE SORCEA END' 

RETURN 
40 CONTINUE 

PRINT *, 'THE NO. OF SOURCE ELEMENTS = ' ,MSAX(ISA) 

PRINT *, 'NO. OF SOURCE ELEMENTS PER ARRAY EXCEEDS NESAX' 

CALL DONEPL 

STOP 

50 CONTINUE 

PRINT *, 'THE TOTAL NO. OF SOURCES = ',IST 
PRINT *, 'TOTAL NO. OF SOURCES EXCEEDS NSX' 

CALL DONEPL 

STOP 

END 



C 

c*** 

C*** 

r*** 

C+++ 

C+++ 

C+++ 

c+++ 

C+++ 



C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



SUBROUTINE SORCEG 



THIS SUBROUTINE READS IN THE SOURCE DATA AND COMPUTES THE CORNERS 
LOCATIONS 

IT ALSO COMPARES THE MIN. & MAX. X, Y & Z VALUES 



SUBROUTINE SORCEG(XSC,YSC,ZSC) 

PARAMETER (NSX=1000) 

REAL XSC(NSX,4).YSC(NSX,4) ,ZSC(NSX,4) ,X(4) ,Y(4) ,Z(4) 
REAL COMPSG(3,2) ,R(3,3),RS(3,3) 

COMMON /RTDATA/ XCR, YCR, ZCR, R 
COMMON /SCALER/ FREQ , IUN , IUS , SCALEF 
COMMON / SGDATA/ ISG , COMPSG 
WRITE (2,*) 

WRITE (2,*) 'SUBROUTINE SORCEG' 



READ IN DATA 



WRITE (2 *) ' ISG =' , ISG 

READ ( 1 , *) XCS , YCS ,ZCS 

WRITE (2,*) ' XCS =',XCS,' YCS =',YCS,' ZCS =',ZCS 

CALL SCALES (1, ' G ' , XCS , YCS , ZCS) 

READ ( 1 , *) THSZ , PHSZ , THSX , PHSX 
WRITE (2,*) ' THSZ= ' , THSZ , ' PHSZ=',PHSZ 

WRITE (2 *) ' THSX=' ,THSX, ' PHSX=',PHSX 

READ ( 1 , *) DUMMY, HS, HAWS 

WRITE (2 *) ' HS = ' ,HS , ' HAWS =',HAWS 

READ(1 , *) DUMMY, DUMMY 

COMPUTE THE AXES ROTATION MATRIX FOR SOURCE CO-ORDINATE 



CALL ROTATE (THSZ , PHSZ , THSX , PHSX , RS ) 
COMPUTE THE CORNERS OF THE SOURCE 



I 
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X(1 )=HAWS/ 2 . 0 

X(2)=HAWS/2.0 

X(3)=-1.0*HAWS/2.0 

X(4)=-1.0*HAWS/2.0 

Y(1)=0.0 

Y( 2 )=0.0 

Y(3 5=0 . 0 

Y(45=0.0 

Z(l5=HS/2.0 

Z(2)=-1.0*HS/2.0 

Z(3)=-1.0*HS/2.0 

Z(4)=HS/2.0 

CALL SCALES(4, 'S' ,X,Y,Z) 

C+++ 

C+++ TRANSFER THE CORNER LOCATIONS TO GLOBAL CO-ORDINATE 
C+++ 

CALL ROTRAN (4 , X , Y , Z , XCS , YCS , ZCS , RS ) 

CALL ROTRAN ( 4 , X , Y , Z , XCR , YCR , ZCR , R ) 

C+++ 

C+++ COMPARE THE MIN. & MAX. X, Y & Z VALUES 
C+++ 

DO 10 1=1,4 
XSC(ISG, l)=X(I) 

YSC(ISG,i5=Y(i5 

zsc(isg,i5=z(i5 

10 CONTINUE 

CALL COMP AR ( 4 , COMPSG , X , Y , Z ) 

WRITE (2,*) 'SUBROUTINE SORCEG END' 

RETURN 

END 

C 

C*** 

C*** SUBROUTINE VANGLE 
r*** 

C+++ 

C+++ THIS SUBROUTINE DETERMINES THE VIEWING ANGLES 
C+++ 

SUBROUTINE VANGLE (TITLE , LENGTH) 

REAL LENGTH ( 3 ) 

CHARACTER ANS$ 

CHARACTER*60 TITLE 

COMMON /VUDATA/ PHI , THETA, RADIUS 

WRITE (2,*) 

WRITE (2, *) 'SUBROUTINE VANGLE' 

C+++ 

C+++ COMPUTE THE MINIMUM DISTANCE RADMIN AWAY FROM THE CENTER OF 
C+++ THE WORK BOX 

C+++ ALSO ENSURE THAT THE RADIUS CHOSEN IS GREATER THAN RADMIN 
C+++ 

RADMIN=1 . 0 
DO 10 1=1,3 

IF (LENGTH(I)*5.0 ,GT. RADMIN) RADMIN=REAL( INT (LENGTH (I )*5 . 0) ) 
10 CONTINUE 

IF (RADIUS .LT. RADMIN) RADIUS=RADMIN 
20 CONTINUE 

CALL FRTCMS ( ' CLRSCRN ') 

PRINT *, TITLE 
PRINT * 

PRINT *, 'DO YOU WANT TO CHANGE THE VIEWING ANGLE?' 

PRINT *, 'THE PRESENT SETTINGS (FROM CENTER OF WORK BOX) ARE:' 

PRINT *, ' PHI IN DEG =',PHI 

PRINT *, ' THETA IN DEG =' , THETA 

PRINT *, ' RADIUS IN M =' , RADIUS 

PRINT *, 'ANSWER Y OR N PLEASE' 

READ(5, ' (Al) ' ) ANS$ 

IF (ANS$ .NE. ' Y ' .AND. ANS$ .NE. ' N ' ) GOTO 20 
IF ( ANS$ .EQ. 'N ' ) GOTO 40 
PRINT *, 'ENTER PHI IN DEG' 

READ(5,*) PHI 
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PRINT * 'ENTER THETA IN DEG 1 
READ (5,*) THETA 

30 PRINT *, 'ENTER RADIUS IN M 1 
READ ( 5 , *) RADIUS 
IF (RADIUS .GE. RADMIN) GOTO 40 
PRINT *, 'RADIUS MUST BE GREATER OR EQUAL TO 
GOTO 30 
40 CONTINUE 

CALL FRTCMS ( 1 CLRSCRN ' ) 



WRITE (2,*) 


1 PHI =' 


, PHI 


WRITE (2,*) 


1 THETA 


=' , THETA 


WRITE (2,*) 


1 RADIUS 


=' , RADIUS 


WRITE (2,*) 


1 SUBROUTINE 


VANGLE END 


RETURN 






END 







, RADMIN, ' M' 
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4. SA CARD GENERATING CODE (GTDSA) 
a. Flow Chart 




SUiR. tfCOC 
rNn In 
fCC-ISC cownai 



m 

SUBR. 

raao* In geo*. 

of array 4 
c**put*« carter 
location* 



1 > 



SUBR TARgfL 

••Tecta 
tapering typ# 



-■SUBfL- I Attft. . 

CO«OUte* th* 

E-ftold •«p, 
tapering ret16 



zn 

SUBR. EflELD 
compute* 
E-Ma Id 
Mgnftuo* for 
-Mg? J.' ffP tt. 



?L'9B. VRITEC 
■rite* output 
cata file 
according to 
NEC-3SC fonata 



1 1 
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b. Program Listing 



C 

c 

c 

c 

c 

c 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



********************** 

* GTDSA * 

* BY KOH WEE JIN * 

* ON OCT 10 1986 * 

* AT NPS * 

********************** 



THIS PROGRAM IS TO CREATE THE SA CARD NEEDED BY NEC-BSC CODE 
FOR PLANAR ARRAY 

It is capable of producing the following Taperings for both the 
Xt and Zt directions : 

Uniform 

Cosine 

Cosine square 
Triangular and 
binomial 



with either TE 10 mode or Uniform distribution for all elements 

The code uses the following assumptions: 

Simple rectangular array grid 

All elements are confined to Xt-Zt plane 

All elements have the same length, width and orientation 

E field is orientated along Xt direction 

The variabel names used and their meanings are listed below: 

Maximum Dimensions 



NESAX : Max. No. of elements per array 
NELSAX : Max. No. of elements per side 

Note that NELSAX is needed to limit the 
dimensions of the array variables 

Units and Wavelength 

IUN : Integer used to define geometry unit 
IUS : Integer used to define source unit 

UNITS : Dimension of unit specified 
UNITS : Unit conversion ratio to m 



FREQ : Frequency in GHz 
LAMBDA : Wavelength in m 

Rotation Angles and Matrix 



THSZ 

PHSZ 

THSX 

PHSX 

ROT 

LORTHO 



Theta angle of Zt axis w. r. t. global co-ord. 
Phi angle of Zt axis w. r. t. global co-ord. 
Theta angle of Xt axis w. r. t. global co-ord. 
Phi angle of Xt axis w. r. t. global co-ord. 
Axes rotation matrix 

Logic to check if Xt and Zt are orthogonal 



Array Area 

XALENG : Width of array area including edge spacing 
ZALENG : Length of array area including edge spacing 



XCAREA : X co-ord. of the center of array area 
YCAREA : Y co-ord. of the center of array area 
ZCAREA : Z co-ord. of the center of array area 
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C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 

C+++ 



C+++ 

C+++ 

C+++ 



Element parameters 

NX 
NZ 

A 
B 

XCE 
YCE 
ZCE 

XELENM 
ZELENM 

XELENG 

ZELENG 



XESPAC 

ZESPAC 

Tapering 

IXTAP 

IZTAP 



AMPPED 
XAMP 
ZAMP 

Power and E 
PTOTAL 
E 

ERATIO 
ESQSUM 
EMAX 

Constants 

C 

DTR 
ZO 

PARAMETER (NESAX=45000) 

PARAMETER (NELSAX = 270) 

REAL XCE (NELSAX,NELSAX) ,YCE (NELSAX,NELSAX) ,ZCE(NELSAX,NELSAX) 
REAL XAMP (NELSAX) , ZAMP (NELSAX) 

REAL E (NELSAX, NELSAX) , LAMBDA, ROT (3, 3) 

CHARACTER* 1 ANS$ 

LOGICAL LORTHO 

COMMON /E ANGLE/ THSZ,PHSZ,THSX,PHSX 

COMMON /GENDAT/ C,PI,DTR,Z0 

COMMON /NECDAT/ IUN, IUS , FREQ, LAMBDA 

COMMON /RTDATA/ ROT , XCARE A , Y CAREA , ZCAREA 

COMMON /XZDATA/ NX,NZ,XALENG,ZALENG,A,B, IDIST 

C=0. 2997925 

PI=3. 1415926 

DTR=PI/180 .0 

Z0=376 . 7 

Get NEC-BSC commands on units and frequency 
CALL NECCMD 



: No. of elements along Xt direction 
s No. of elements along Zt direction 

: Length of element in Zt direction 
: Width of element in Xt direction 

: X co-ord. of the center of each element 

. Y co-ord. of the center of each element 

: Z co-ord. of the center of each element 

: Element width in Xt direction in m 

: Element length in Zt direction in m 

s Element width in Xt direction in the same unit 
specified for geometry (IUN) 
s Element length in Zt direction in the same unit 
specified for geometry (IUN) 

: Element spacing in Xt direction 
: Element spacing in Zt direction 



: Integer used to define the type of element 
tapering along Xt direction 
: Integer used to define the type of element 
tapering along Zt direction 

s Amplitude ratio of pedestal 
: Amplitude tapering along Xt direction 
: Amplitude tapering along Zt direction 

Field 

: Total input power in watts 

s E field for each element in V/m 
s E field amplitude ratio for each element 
: Sum of the square of E field ratio 
s Max. E field in V/m 



s Speed of light divided by 1E9 
: Degree to radian conversion ratio 
s INtrinsic wave impedance in ohm 
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C+++ 

C+++ Compute the center locations of each elements 
C+++ 

CALL ELMCEN (XCE ,YCE ,ZCE) 

C+++ 

C+++ Select the tapering types and pedestal height 
C+++ 

CALL TAPSEL(IXTAP , IZTAP , AMPPED) 

C+++ 

C+++ Compute the amplitude tape rings 
C+++ 

CALL TAPER ( I XTAP , 1 X 1 , XAMP , XCE ) 

CALL TAPER ( I ZTAP , ' Z ' , ZAMP , ZCE ) 

C+++ 

C+++ Computer the E-field magnitude 
C+++ 

CALL EFIELD (XAMP, ZAMP, E, AMPPED) 

C+++ 

C+++ Write out the SA card data 
C+++ 

CALL WRITED(XCE , YCE , ZCE , E) 

STOP 

END 

C 



^^^^•k-k-k-k-k'k'k’k'k'k’k'k’k’k'k'k’k'k’k 

* SUBROUTINE EFIELD * 

****** -kit *********** 

This subroutine computes the E field in V/m for each element 

SUBROUTINE EFIELD (XAMP , ZAMP , E , AMPPED ) 

PARAMETER (NESAX=45000) 

PARAMETER (NELSAX=270) 

REAL XAMP(NELSAX) ,ZAMP(NELSAX) , LAMBDA 
REAL ERATIO(NELSAX,NELSAX) ,E (NELSAX,NELSAX) 

REAL UNITS(0:3) 

CHARACTER* 1 ANS$ 

COMMON /GENDAT/ C,PI,DTR,Z0 
COMMON /NECDAT/ I UN , I US , FREQ , LAMBDA 
COMMON /XZDATA/ NX ,NZ . XALENG , ZALENG , A,B , IDIST 
DATA UNITS ( 1 ) , UNITS ( 2 ) , UNITS (3 ) /I .0,0. 3048 , 0 . 0254/ 

UNITS (0)=LAMBDA 
WRITE (2,*) 

WRITE (2,*) 'Subroutine EFIELD' 

WRITE (2,*) ' I Xamp ' 

DO 10 1=1, NX 

WRITE (2,*) I , XAMP ( I ) 

10 CONTINUE 

WRITE (2,*) ' I Zamp' 

DO 20 1=1, NZ 

WRITE (2,*) I , ZAMP ( I ) 

20 CONTINUE 
C+++ 

C+++ Convert element lengths to m 
C+++ 

XELENM=B*UNITS ( IUS ) 

ZELENM=A*UNITS ( IUS ) 

C+++ 

C+++ Compute the relative amplitude 
C+++ 

ESQSUM=0 . 0 
DO 40 1=1, NZ 

DO 30 J=1,NX 

ERATIO (I , J) =XAMP ( J ) *ZAMP ( I ) * ( 1 . 0- AMPPED ) + AMPPED 
ESQSUM=ESQSUM+ERATIO ( I , J ) *^2 
30 CONTINUE 
40 CONTINUE 



C 
C 
C 
C 

C+++ 

C+++ 

C+++ 
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WRITE (2,*) ■ ESQSUM = 1 ,ESQSUM 

C+++ 

C+++ Enter total input power 
C+++ 

PRINT * 

PRINT *, 'Enter total input power in Watt' 

PRINT * 

READ (5 ,*) PTOTAL 
PRINT * 

WRITE (2 , *) ' PTOTAL = PTOTAL , 1 W' 

C+++ 

C+++ Enter element distribution type 

C+++ TE 10 of Uniform 

C+++ 

50 CONTINUE 
PRINT * 

PRINT *, 'Enter element distribution type' 

PRINT *, '0) Uniform distribution' 

PRINT *, '1) TE 10 mode' 

PRINT * 

READ (5 ,*) IDIST 

IF (IDIST .NE. 0 .AND. IDIST .NE. 1) GOTO 50 
IF (IDIST .EQ. 0) GOTO 55 

C+++ 

C+++ Compute maximum E field 
C+++ 

EMAX=SQRT(PT0TAL*4.0*Z0/ZELENM/XELENM/ESQSUM) 

GOTO 58 
55 CONTINUE 

EMAX=SQRT(PTOTAL*Z0/ZELENM/XELENM/ESQSUM) 

58 CONTINUE 

WRITE (2,*) ' EMAX = ' , EMAX, ' V/m 1 

C+++ 

C+++ Check if Emax exceeds 2.9 MV/m 
C+++ 

IF (EMAX .GT. 1.0631E6) THEN 
60 CONTINUE 
PRINT * 

PRINT *, 'Emax = ' , EMAX , ' V/m' 

PRINT *, 'Emax is greater than the max. power flux constrain' 
PRINT *, 'of 300 KW/sq-cm or equivalent to 10.631 KV/cm' 

PRINT *, 'Do you want to continue?' 

PRINT * 

READ (5, ' (Al) ' ) ANS$ 

IF (ANS$ .NE. ' N ' .AND. ANS$ .NE. ' Y ' ) GOTO 60 
IF (ANS$ .EQ. 'N') GOTO 90 
END IF 

C+++ 

C+++ Compute individual excitations 
C+++ 

DO 80 1=1, NZ 

DO 70 J=1 ,NX 

E (I , J)=EMAX*ERATIO(I , J) 

70 CONTINUE 
80 CONTINUE 

WRITE (2 ,*) 'Subroutine EFIELD end' 

RETURN 
90 CONTINUE 
STOP 
END 



C 

c 

c 

c 

c 



C+++ 

C+++ 

C+++ 

C+++ 



*********************** 

* SUBROUTINE ELMCEN * 

This subroutine reads in the geometry data of the array elements 
and computes the center locations of each element 
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SUBROUTINE ELMCEN ( XCE , YCE , ZCE ) 

PARAMETER (NESAX=45000) 

PARAMETER \NELSAX = 270) 

REAL XCE (NELS AX , NELS AX ) , YCE (NELSAX,NELSAX) ,ZCE(NELSAX / NELSAX) 
REAL ROT (3,3), LAMBDA , UNITS ( 0 : 3 ) 

CHARACTER* 10 UNITS (0:3) 

CHARACTER* 1 ANS$ 

LOGICAL LORTHO 

COMMON /E ANGLE/ THSZ, PHSZ, THSX, PHSX 

COMMON /GENDAT/ C,PI,DTR,Z0 

COMMON /NECDAT/ IUN , I US , FREQ , LAMBDA 

COMMON /RTDATA/ ROT , XCAREA , Y CAREA , ZCARE A 

COMMON /XZDATA/NX , NZ , XALENG , ZALENG , A , B , IDIST 

DATA UNITS ( 1 ) , UNITS ( 2 ) , UNITS ( 3 ) / 1 . 0 , 0 . 3048 ,0.0254/ 



C+++ 

C+++ 

C+++ 



DATA UNITS/ ' wavelength m 
UNITS (0)=LAMBDA 
WRITE (2,*) 

WRITE (2,*) 'Subroutine ELMCEN 



ft 



in 



10 



Enter array element orientation and compute axes rotation matrix 

CALL FRTCMS ( ' CLRSCRN ') 

CONTINUE 

PRINT *, 'Enter Array Element orientation (in deg) that 1 
PRINT *, 'defines Zt and Xt directions as used in NEC-BSC code' 
PRINT *, 'i.e. THSZ , PHSZ , THSX , PHSX ' 

PRINT * 

READ (5,*) THSZ, PHSZ, THSX, PHSX 
CALL ROTATE (LORTHO) 

IF (LORTHO) GOTO 20 
PRINT * 

THSZ 
PHSZ 
THSX 
PHSX 



20 



PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
GOTO 10 
CONTINUE 



k 

*' 

*' 

*' 

*' 

*' 

*' 



1 ,THSZ , 1 
1 , PHSZ, 1 



1 ,THSX, 1 
1 ,PHSX, 1 

Axes specified are not orthogonal 1 
Re-enter the angles' 



deg' 

deg 1 

deg' 

deg' 



C+++ 

C+++ 

C+++ 



WRITE 


(2,*; 


1 1 THSZ 


= 1 ,THSZ, 


WRITE 


(2 , * 


1 1 PHSZ 


= ' , PHSZ , 


WRITE 


2 * 


1 1 THSX 


= ' ,THSX, 


WRITE 


( 2 ,*j 


1 1 PHSX 


= ' ,PHSX, 



C+++ 

C+++ 



deg' 
deg' 
deg' 
deg' 

Enter center location of the Array Area (in Global co-ord) 

CALL FRTCMS ( ' CLRSCRN ' ) 

PRINT * 

PRINT *, 'Enter center co-ordinate of Array Area (X,Y,Z) ' 

PRINT * 1 * ' ' 

PRINT * 

READ (5 
WRITE ( 

WRITE ( 

WRITE ( 

Enter number of elements 



'The unit is in ', UNITS (IUN) 



,*), 
2 ,* 


XCAREA , Y CARE A , Z CAREA 
) ' XCAREA = ', XCAREA, UNITS 1 


:iun 


2 * 


1 ' Y CAREA = 


' , YCAREA, UNITS i 


IUN 


2 , *! 


) ' ZCAREA = 


' , ZCAREA, UNITS 1 


IIUN 



C+++ 

CALL FRTCMS ( ' CLRSCRN ' ) 

30 CONTINUE 

PRINT *, 'Enter No. of elements along Xt axis' 

PRINT * 

READ (5,*) NX 
IF (NX .GT. NELS AX) THEN 
PRINT * 

PRINT *, 'No. of elements along Xt axis = ' ,NX 

PRINT *, 'No. of elements along Xt axis exceeds NELSAX' 

PRINT *, 're-enter the No. of elements' 
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35 



PRINT * 

GOTO 30 
END IF 
CONTINUE 
PRINT * 

PRINT *, 'Enter No. of elements along Zt axis' 
READ (5,*) NZ 
IF (NZ .GT. NELSAX) THEN 
PRINT * 



'No. of elements along Zt axis = ' ,ZX 

'No. of elements along Zt axis exceeds NELSAX' 

're-enter the No. of elements' 



PRINT * 

PRINT * 

PRINT * 

PRINT * 

GOTO 35 
END IF 

NTOTAL=NX*NZ 

IF (NTOTAL .GT. NESAX) THEN 
PRINT * 

'Total No. of elements = ', NTOTAL 
'Total No. of elements exceeds NESAX 
're-enter the No. of elements' 



C+++ 

C+++ 

C+++ 



PRINT * 
PRINT * 
PRINT * 
PRINT * 
GOTO 30 
END IF 
WRITE (2,*) 
WRITE (2,*) 



NX 

NZ 



= ' , NX 
= ' ,NZ 



Enter the size of element 



40 



CONTINUE 

CALL FRTCMS ( 1 CLRSCRN ' ) 

PRINT * 

Enter the size of the element 1 
Note that all elements. are assumed to have' 



PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 



'the same size' 

'Enter the element length along Xt axis direction' 
'The unit is in ',UNIT$(IUS) 



C 

C 

C 

C 

c 



READ (5,*) B 
PRINT * 

PRINT *, 'Enter the element length along Zt axis direction' 
PRINT *, 'The unit is in ', UNITS (IUS) 

PRINT * 

READ (5 ,*) A 
IF (B .GT. A) THEN 
PRINT * 

PRINT *, 'THE ELEMENT LENGTH = ' , A,UNIT$(IUS) 

PRINT *, 'THE ELEMENT WIDTH = ' ,B,UNIT$(lUS) 

PRINT * 'THE WIDTH MUST BE SHORTER THAN THE LENGTH' 



c 


PRINT *, 


'RE-ENTER 


ELEMENT SIZE' 




c 


PRINT * 








c 


GOTO 40 








c 


END IF 










WRITE (2,*) 


' A 


= ' , A, UNITS 


(IUS) 




WRITE (2,*) 


' B 


= ' ,B, UNITS 


(IUS) 



C+++ 

C+++ 

C+++ 



Check if TE 10 mode exist 

ZELENG=A*UNITS ( IUS ) 

IF (ZELENG .LT. LAMBDA/ 2.0) THEN 



PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 



* 
* 
*' 
*' 
*' 
*' 
* ' 
*' 



m‘ 
1 1 



'The operating freq = 'FREQ,' GHz' 

'The corresponding wavelength = ', LAMBDA, 

'The length of array element = ', ZELENG, ' m 1 
'For TE 10 mode to exist, the length of the element' 
'must be greater than or equal to half wave length 1 
'Re-enter element size' 
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GOTO 40 
END IF 

C+++ 

C+++ Convert the element lengths to the geometry units 
C+++ 

XELENG=B*UNITS ( IUS ) /UNITS ( IUN) 

ZELENG=A*UNITS ( IUS ) /UNITS ( IUN) 

WRITE (2,*) 1 XELENG = ' ,XELENG, UNITS (IUN) 

WRITE (2,*) ' ZELENG = ' ,ZELENG, UNITS (IUN) 

C+++ 

C+++ Enter gap size between elements 
C+++ 

CALL FRTCMS ( 1 CLRSCRN ') 

PRINT *, 'Enter the gap size between elements in Xt direction 1 
PRINT *, 'The unit is in ' , UNITS (IUS) 

PRINT * 

READ (5,*) XGAP 
PRINT * 

PRINT *, 'Enter the gap size between elements in Zt direction' 
PRINT *, 'The unit is in ' ,UNIT$(IUS) 

PRINT * 

READ (5,*) ZGAP 

XGAP1=XGAP*UNITS ( IUS ) /UNITS ( IUN) 

ZGAP1=ZGAP*UNITS ( IUS ) /UNITS ( IUN) 

XESPAC=XGAP1+XELENG 
ZESPAC=ZGAP1+ZELENG 
XALENG=XESPAC*REAL(NX+1 ) 

IF (NX .EQ. 1) XALENG=XELENG 
ZALENG=ZESPAC*REAL(NZ+1 ) 

IF (NZ .EQ. 1) ZALENG=ZELENG 



WRITE 


2,* 


1 1 XGAP 


= 


' , XGAP, UNITS (IUS) 


WRITE 


2,* 


1 1 ZGAP 


= 


' , ZGAP, UNITS ( IUS ) 


WRITE 


2,* 


1 ' XESPAC 


= 


' , XESPAC, UNITS 1 


IUN) 


WRITE 


(2 , * 


1 ' ZESPAC 


= 


' , ZESPAC, UNITS! 


IUN 


WRITE 


2,* 


1 ' XALENG" 


=r 


' ,XALENG, UNITS! 


IUN 


WRITE 


2,* 


1 ' ZALENG 


= 


' , ZALENG , UNITS < 


[IUN) 



C+++ 

C+++ Determine the center locations of each element in Xt-Zt plane 
C+++ 

CXXX WRITE (2,*) ' I J XCE ' , 

CXXX *' YCE ZCE ' 

DO 70 1=1, NZ 

Z=-l .0*ZALENG/2.0+I*ZESPAC 
IF (NZ .EQ. 1) Z=0.0 
DO 60 J=1 ,NX 

X=- 1 . 0*XALENG/2 . 0+J*XESPAC 
IF (NX .EQ. 1) X=0.0 
XCE ( I , J)=X 
YCE(I , J)=0.0 
ZCE(I,J)=Z 

CXXX WRITE (2 ,*) I , J,XCE (I , J) ,YCE(I , J) ,ZCE (I , J) 

60 CONTINUE 
70 CONTINUE 

WRITE (2 , *) 'Subroutine ELMCEN end' 

RETURN 

END 



C 

C 

C 

C 

C 



C+++ 

C+++ 



*********************** 

* SUBROUTINE NECCMD * 

*********************** 

SUBROUTINE NECCMD 
REAL LAMBDA 

COMMON /GENDAT/ C,PI,DTR,Z0 
COMMON /NECDAT/ I UN, I US , FREQ, LAMBDA 
WRITE (2,*) 

WRITE (2 ,*) 'Subroutine NECCMD' 

The following NEC-BSC commands are needed to determine 
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C+++ 

C+++ 

C+++ 



c 
c 
c 
c 
c 

C+++ 

C+++ 

C+++ 



c??? 

c??? 

c??? 

c??? 

c??? 

c??? 

c??? 



the units used by the whole data set 

They should be consistance with the rest of the data 



PRINT 


* 

/ 


'Enter the fo 


PRINT 


* 

/ 


'UN, US, FR' 


PRINT 


* 

/ 


'where UN and 


PRINT 


/ 


' UN = 1 


PRINT 


* 

/ 


' UN = 2 


PRINT 


* 

/ 


' UN = 3 


PRINT 

PRINT 


* 

* 

/ 


' US = 0 


PRINT 


* 

/ 


' US = 1 


PRINT 


* 

/ 


' US = 2 


PRINT 


* 

. / 


' US = 3 



PRINT 

PRINT 

PRINT 



* 

* 

*' 



unit in m' 
unit in ft 1 
unit in in 1 

unit in wavelength 1 
unit in m 1 
unit in ft 1 
unit in in 1 



‘FR is frequency in GHz 1 



WRITE 
WRITE 
WRITE 
WRITE 
RETURN 
END 



C/FREQ 


2,*) 


' I UN = ' , 


2 '1\ 


' FREQ = 


2 ,*) 


' LAMBDA = 


2,* 


'Subroutine ] 



IUS = 

, 1 GHz’ 
)A, 1 m 1 



1 ,IUS 



*********************** 

* SUBROUTINE ROTATE * 

* -k "k * * * * -k * * * * 7c -k * * * * * * -k * * 



This subroutine computes the axes rotation matrix 

SUBROUTINE ROTATE ( LORTHO ) 

REAL ROT (3, 3) 

LOGICAL LORTHO 

COMMON /E ANGLE/ THSZ,PHSZ,THSX,PHSX 
COMMON /GENDAT/ C,PI,DTR,Z0 
COMMON /RTDATA/ ROT , XCARE A , YCAREA , ZCAREA 
LORTHO= . TRUE . 

WRITE (2,*) 

WRITE (2,*) 'Subroutine ROTATE' 

THZ=THSZ*DTR 

THX=THSX*DTR 

PHZ=PHSZ*DTR 

PHX=PHSX*DTR 

ROT ( 1 , 1 ) =COS ( PHX ) *S IN ( THX ) 

ROT (2, 1)=SIN(PHX)*SIN(THX) 

ROT (3 , 1 )=COS ( THX) 

ROT (1 , 3 )=COS (PHZ)*SIN (THZ) 

ROT (2 , 3 )=SIN(PHZ )*SIN(THZ) 

ROT (3 ,3)=COS(THZ) 

Dl=ROT(l,3)*ROT(l,l)+ROT(2,3)*ROT(2,l)+ROT(3,3)*ROT(3,l) 

IF (ABS(Dl) .GT. 0.1) GOTO 10 
ROT (1,1 )=ROT (1,1) -ROT ( 1 , 3 ) *D1 
ROT ( 2 , 1 ) =ROT ( 2 , 1 ) -ROT ( 2 , 3 ) *D1 
ROT (3,1 )=ROT (3,1 ) -ROT(3 , 3 )*D1 

D2=SQRT ( ROT (1,1) ^ROT (1,1) +ROT (2,1) *R0T (2,1) +R0T (3,1) *ROT (3,1)) 
ROT (1,1) =ROT ( 1 , 1 ) / D2 
ROT (2 , 1 )=ROT( 2 , 1 ) /D2 
ROT (3,1) =ROT ( 3 , 1 ) / D2 

ROT (1,2) =ROT (2,3) *ROT (3,1) -ROT (3,3) *ROT (2,1) 

ROT (2,2) =ROT (1,1) *ROT (3,3) -ROT (1,3) *ROT (3,1) 
ROT(3 / 2)=ROT(l,3)*ROT(2,l)-ROT(2,3)*ROT(l,l) 

WRITE (2,*) 'AXIS ROTATION MATRIX' 

WRITE (2,*) ROT (1,1), ROT (1,2), ROT (1,3) 

WRITE (2,*) ROT (2,1), ROT (2,2), ROT (2,3) 

WRITE (2,*) ROT (3,1), ROT (3,2) , ROT (3,3) 
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WRITE (2,*) 'Subroutine ROTATE end' 
RETURN 
10 CONTINUE 

LORTHO= . FALSE . 

WRITE (2 ,*) 'Subroutine ROTATE end' 

RETURN 

END 



* * * * * * -k * * * * * * * * * * * * * * * * 

* SUBROUTINE ROTRAN * 

*********************** 

This subroutine rotates and translate co-ordinates 

SUBROUTINE ROTRAN (X,Y,Z) 

PARAMETER (NELSAX=270) 

REAL X(NELSAX , NELSAX) , Y (NELSAX , NELSAX ) , Z (NELSAX,NELSAX) , ROT(3 , 3) 
COMMON /RTDATA/ ROT , XCAREA , YCAREA , ZCAREA 
COMMON /XZDATA/ NX,NZ , XALENG, ZALENG, A, B , IDIST 
DO 20 1=1, NZ 

DO 10 J=1 ,NX 
X1=X(I,J) 

Y1=Y(I, J) 

Z1=Z(I,J) 

X ( I , J ) =ROT ( 1 , 1 ) *X1 +ROT ( 1 , 2 ) *Y1 +ROT (1,3) *Z1 +XCAREA 
Y ( I , J ) =ROT ( 2 , 1 ) *X1 +ROT (2,2) *Y1+R0T (2,3) *Z1+YCAREA 
Z(I, J)=R0T(3,1)*X1+R0T(3,2)*Y1+R0T(3,3)*Z1+ZCAREA 
10 CONTINUE 
20 CONTINUE 
RETURN 
END 



C 

c 
c 
c 

C+++ 

C+++ 

C+++ 



c 
c 
c 
c 
c 

C+++ 

C+++ 

C+++ 



10 



20 



C+++ 

C+++ 

C+++ 



********************** 

* SUBROUTINE TAPER * * 
********************** 



This subroutine computes the distribution in Xt or Zt axis 

SUBROUTINE TAPER( ITAP , XZ$ , AMP , CENTER) 

PARAMETER (NESAX=45000) 

PARAMETER (NELSAX=270) 

REAL AMP(NELSAX) , CE (NELSAX) , CENTER(NELSAX, NELSAX) 
CHARACTER* 1 XZ$ 

COMMON /GENDAT/ C,PI,DTR,Z0 

COMMON /XZDATA/ NX, NZ, XALENG, ZALENG, A,B, IDIST 
AMPMAX=0 . 0 

IF (XZ$ .EQ. 'X') THEN 
N=NX 

ALENG=XALENG 
DO 10 1=1, NX 

CE ( I )=CENTER (1,1) 

CONTINUE 

ELSE 

N=NZ 

ALENG=ZALENG 
DO 20 1=1, NZ 

CE ( I ) =CENTER (1,1) 

CONTINUE 
END IF 

IF (ITAP .EQ. 1 .OR. ITAP .EQ. 2) GOTO 40 
IF (ITAP .EQ. .3) GOTO 50 
IF (ITAP .EQ. 4) GOTO 60 

Uniform tapering 



DO 30 1=1, N 
AMP ( I )=1 . 0 
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30 CONTINUE 
RETURN 

C+++ 

C+++ Cosine and cosine square tapering 

C+++ 

40 CONTINUE 
J=ITAP 
DO 45 1=1 ,N 

AMP ( I ) =COS ( CE ( I ) /ALENG*PI ) ** J 
IF (AMPMAX .LT. AMP(l)) AMPMAX=AMP ( I ) 
45 CONTINUE 
GOTO 70 

C+++ 

C+++ Triangular tapering 

C+++ 

50 CONTINUE 

DO 55 1=1 ,N 

AMP (I )=l-2 . 0*ABS (CE (I ) )/ALENG 
IF (AMPMAX .LT. AMP(I)) AMPMAX=AMP ( I ) 
55 CONTINUE 
GOTO 70 



C+++ 

C+++ Binomial tapering 
C+++ 

60 CONTINUE 

DO 66 1=0 ,N-1 
' COEFF=l .0 

IF (I .EQ. 0 .OR. I .EQ. N-l) GOTO 64 
DO 62 J=1,I 

COEFF=COEFF*(N-l-J+l)/J 
62 CONTINUE 

64 CONTINUE 

IF (AMPMAX .LT. COEFF) AMPMAX=COEFF 
AMP(l+l)=COEFF 
66 CONTINUE 
C+++ 

C+++ Normalization 



C+++ 

70 CONTINUE 

DO 80 1=1 ,N 

AMP ( I ) =AMP ( I ) / AMPMAX 
80 CONTINUE 
RETURN 
END 



C--- 

C 

C 

C 

C 

C+++ 



C+++ 

C+++ 

C+++ 
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■k-k'k’k'k'k'k-k'k'k^’k'k’k'k'k^'k'k^'k'k-k 



* SUBROUTINE TAPSEL * 
*********************** 



This subroutine selects the Xt and Zt axis tapering type 
and the tapering pedestal height 

SUBROUTINE TAPSEL (IXTAP , IZTAP , AMPPED) 

WRITE (2,*) 

WRITE (2,*) 'Subroutine TAPSEL 1 
CONTINUE 

CALL FRTCMS ( 1 CLRSCRN ' ) 

PRINT *, 'Select one of the tapering type in Xt direction' 
PRINT * 

PRINT *, '0) Uniform' 

PRINT *, 'll Cosine' 

PRINT *, '2) Cosine square' 

PRINT *, '3) Triangular' 

PRINT *, '4) Binomial' 

PRINT * 

PRINT *, 'Enter the number before the tapering type' 

PRINT * 

READ (5, '(II)') IXTAP 
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nnnnn 



IF (IXTAP .LT. 0 .OR. IXTAP ,GT. 4) GOTO 10 
20 CONTINUE 
PRINT * 

PRINT * 

PRINT *, 'Select one of the tapering type in Zt direction 1 
PRINT * 

PRINT *, '0) Uniform' 

PRINT *, '1) Cosine' 

PRINT *, '2) Cosine square' 

PRINT *, '3) Triangular' 

PRINT *, '4) Binomial' 

PRINT * 

PRINT *, 'Enter the number before the tapering type 1 
PRINT * 

READ (5, 1 (II) ' ) IZTAP 

IF (IZTAP .LT. 0 .OR. IZTAP .GT. 4) GOTO 20 
CALL FRTCMS ( ' CLRSCRN ') 

PRINT *, 'Enter pedestal height relative to max amplitude' 
PRINT *, 'Enter 0.0 if no pedestal is needed' 

PRINT * 

READ (5,*) AMPPED 
WRITE (2,*) ' IXTAP = \IXTAP 

WRITE 2,* ' IZTAP = ', IZTAP 

WRITE (2,*) ' AMPPED = ', AMPPED 

WRITE (2,*) 'Subroutine TAPSEL end' 

RETURN 

END 



*********************** 

* SUBROUTINE WRITED * 

■k'k’k'k’k'k’kjs'kjsjs'k'kjzjz'k’kjs'kjs'k'k'k 

C+++ 

C+++ This subroutine writes the output data in the format 

C+++ set by NEC-BSC code 

C+++ 

SUBROUTINE WRITED (XCE , YCE , ZCE , E ) 

PARAMETER (NESAX=45000) 

PARAMETER (NELSAX=270) 

REAL XCE (NELS AX , NELSAX ) ,YCE (NELSAX,NELSAX) ,ZCE(NELSAX,NELSAX) 
REAL E ( NELSAX , NELS AX ), ROT (3, 3) 

COMMON /EANGLE/ THSZ , PHSZ , THSX , PHSX 

COMMON /RTDATA/ ROT , XCAREA , YCAREA , ZCAREA 

COMMON /XZDATA/ NX , NZ , XALENG , ZALENG , A , B , IDIST 

NTOTAL=NX*NZ 

WRITE (1,1000) NTOTAL 

CALL ROTRAN ( XCE , YCE , ZCE ) 

DO 20 1=1, NZ 

DO 10 J=1 ,NX 

WRITE (1,1010) XCE (I, J), YCE (I, J), ZCE (I, J) 

10 CONTINUE 

20 CONTINUE 

WRITE (1,1020) THSZ, PHSZ, THSX, PHSX 

IF (IDIST .EQ. 0) IMS=1 

IF (IDIST .EQ. 1) IMS=3 

WRITE (1,1030) IMS , A , B 

PHASE=0 . 0 

DO 40 1=1, NZ 

DO 30 J=1 ,NX 

WRITE (1,1040) E(I,J), PHASE 
30 CONTINUE 

40 CONTINUE 
1000 FORMAT (16) 

1010 FORMAT (F12. 7, 1H, , F12 . 7. 1H, , F9 .4) 

1020 FORMAT (F6 . 1 , 3(1H , , F6 . 1 ) ) 

1030 FORMAT ( 12, 1H, , F9 . 5, 1H, , F9 .5) 

1040 FORMAT(E12 . 6 , 1H, , F6 . 1 ) 

RETURN 

END 
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APPENDIX B 

RADIATION PATTERN PLOTS 




Figure B.l Element Width = 0.25X, Element Length Variable, 
Z t Spacing = 1 -OX and X t Spacing = 0.5X. 



Ill 




Figure B.2 Element Width = 0.2 5X, Element Length Variable, 
Spacing = 0.5X, Spacing = 0.75X for (a) & (b) 
and Z t Spacing = 0.5X for (c). 
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Figure B.3 Element Width = 0.75X, Element Width Variable, 
X ( Spacing = 0.75X, Z t Spacing = 0.75X. 
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Figure B.4 Element Width = 0.75X, Element Width Variable, 
X t Spacing = 0.5X and Z t Spacing = 0.75X. 




Figure B.5 Element Width = 0.75X, Element Width Variable, 



X t Spacing Variable and Z t Spacing = 0.75X. 
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Figure B.6 Vertical Cut at (p = 0° with Fixed Array Area of 
7.5X x 3.75X. Uniform Tapering in both Directions. 
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Angli 6 (dig) 

(a) Element Size = 0.75^ x 0.37SX 
No, of Elements = 10 x 10 




Angli 6 (dig) 



(c) Element Size = 0.5^. x 0.375^ 
No. of Elements = 15 x 10 




Argi« 8 (dig) 



<b) Element Size = 0.75^ x O.Zsk 
No. of Elements = 10 x 15 




No. of Elements = 15 x 15 



Figure B.7 Vertical Cut at <p = 90° with Fixed Array Area of 
1.5k x 3.75X. Uniform Tapering in both Directions. 
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(a) Element Size = 0.75^ x 0.375X (b) Element Size = 0.75X, x 0.25^. 

No. of Elements - 10 x 10 No. of Elements = 10 x 15 




(c) Element Size = 0.5^ x 0.375X 
No. of Elements = 15 x 10 



(d) Element Size = 0.5A x 0.25 X 
No. of Elements - 15 x 15 



Figure B.8 Vertical Cut at (p = 0° with Fixed Array Area of 
1 . 5 \ x 3.75X. Cosine Tapering in both Directions. 
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Figure B.9 Vertical Cut at <p = 90° with Fixed Array Area of 
7.5X x'3.75X. Cosine Tapering in both Directions. 
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(c) Vertical Cut Near Field Plot 
at (p = 90° 



(d) Horizontal Cut Near Field Plot 
at 0 = 90° 



Figure B.10 Far Field and Near Field Plots with Uniform Tapering. 





(a) Vertical Cut Far Field Plot 
at (p = 90° 



(b) Horizontal Cut Far Field Plot 
at 0 = 90° 
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(c) Vertical Cut Near Field Plot 
at <p = 90° 




Figure B.l 1 Far Field and Near Field Plots with Cosine Tapering. 
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at (p = 90° at 0 = 90° 



Figure B.12 Far Field and Near Field Plots with Cosine Square Tapering 




Figure B.13 Far Field and Near Field Plots with Triangular Tapering. 
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(a) Vertical Cut Far Field Plot 
at (p = 90 0 
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(b) Horizontal Cut Far Field Plot 

at 0 = 90° 
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f Cd#g) 



at (p = 900 



(d) Horizontal Cut Near Field Plot 
at 0 = 900 



Figure B.14 Far Field and Near Field Plots with Binomial Tapering. 




(a) Far Field Horizontal Cut (b) Near Field Horizontal Cut 

at 0 = 90° at 0 = 90 ® 



Figure B.15 Uniformly Tapered Linear TE 10 Array with 
20 Elements of Sizes 0.5X x 0.25X each.- 




la) Far Field Horizontal Cut (b) Near Field Horizontal Cut 

at 0 = 90° at 0 = 90 0 



Figure B.16 Single Uniformly Distributed Slot Antenna with Size 

10. OX x 0.25X. 
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(a) Far Field Horizontal Cut <b) Near Field Horizontal Cut 

at 0 = 90° at 0 = 90° 



Figure B.17 Cosine Tapered Linear TE 10 Array with 
20 Elements of Sizes 0.5X x 0.25X each. 




Figure B.1S Single TE 10 Mode Slot Antenna with Size 10.5X x 0.25X each. 




Figure B.I9 Cosine Tapered Linear Array with 20 Uniformly Distributed 
Elements of Sizes 0.5X x 0.25X each. 




Figure B.20 Far Field Patterns of Simplified Model with 
Uniformly Distributed Elements and Triangular Tapering. 
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(c) Horizontal Cut Pattern (d) Horizontal Cut Pattern 

0 = 134 ° 0 = 1440 



Figure B.21 Near Field Patterns of Simplified Model with 
Uniformly Distributed Elements and Triangular Tapering, r = 3.5m. 




(a) Vertical Cut at <p = 90° (b) Vertical Cut at (p = 0° 




Figure B.22 Near Field Radiation Patterns of Uniform - Uniform Tapered Array. 




Figure B.23 Near Field Radiation Patterns of Cosine - Uniform Tapered Array. 
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(c) Horizontal Cut at 0 = 90° 




Figure B.24 Near Field Radiation Patterns of Uniform - Cosine Tapered Array. 
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Figure B.25 Far Field Radiation Patterns of Array using 
Uniform - Cosine Tapering with Reflector Plate Tilted at 45°. 




Figure B.26 Far Field Radiation Patterns of Array using 
Uniform - Cosine Tapering with Reflector Plate Tilted at 22.5°. 
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Figure B.27 Far Field Radiation Patterns of the Full System Model 
with Reflector Plate Tilted at 45’. 





Figure B.28 Far Field Radiation Patterns of the Full System Model 
with Reflector Plate Tilted at 22.5°. 




figure B.29 Near Field Radiation Patterns of Full Model with Reflector Piate 

Facing <p = 0° and l ilted at 45°. 






Figure B.30 Near Field Radiation Patterns of Full Model with Reflector Plate 

Facing (p = 0° and Tilted at 22.5°. 
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(a) Vertical Cut at (p = 45° 
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(b) Vertical Cut at <p = 135« 




Figure B.3i Near Field Radiation Patterns of Full Model with Reflector Plate 

Facing <p = 45° and Tilted at 45°. 
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Figure B.32 Near Field Radiation Patterns of Full Model with Reflector Plate 

Facing (p = 90° and Tilted at 45°. 
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APPENDIX C 

DERIVATION OF GAIN AND MAXIMUM APERTURE E-FIELD 
RELATIONSHIPS BETWEEN ACTUAL AND SIMPLIFIED MODELS 



The derivation of the gain and maximum aperture E-field relationships between 
the actual linear array using either Uniform or Cosine (TE 10 ) distribution elements and 
those simplified single element models of equivalent length are shown in this appendix. 
Based on the different configurations of the array elements, we subdivide them into 3 
sections, namely, 

• Uniformly tapering TE J0 linear array vs single uniformly distributed slot 
antenna. 

• Cosine tapering TE 10 linear array vs single TE 10 slot antenna. 

• Cosine tapering uniformly distributed linear array vs single TE 10 slot 
antenna. 

For the case when the array is formed by uniformly tapered uniformly distributed 
element, the simplified model of single uniformly distributed aperture will yield exact 
solution and hence no derivation is needed. 



1. UNIFORMLY TAPERED TE10 LINEAR ARRAY VS SINGLE UNIFORMLY 
DISTRIBUTED SLOT ANTENNA 

a. Comparing Maximum Aperture E-Field 

For single Uniform Aperture, 



Where E u * and a ' are the maximum field strength and the length of the uniform 
aperture respectively, b is the width of each element and it is the same for both the 
actual array and the simplified model. Z Q is the intrinsic impedence. 

As for the 20 elements uniformly tapered linear area, we have 



NE r 2 ab 

p _ C C 



in 



4Z n 
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for E c and a c are the maximum field strength and the length of each element of the 
Cosine distributed (TE 10 ) aperture respectively. With the same input power and 
a ' = Na . we have, 

U C’ T 

E c - 2E u • 

Hence, the simplified model gives only half the field strength of the actual 

array. 



b. Comparing Gain 

The boresight far field field strength of an uniformly distributed aperture can 
easily be obtained from the available literature [Refs. 5,6], 



|E',J - 



rX 



where r is the far field distance and X is the wavelength. The boresight far field field 
strength of a TE 10 mode rectangular aperture on the other hand, can be derived from 
the two-dimensional Fourier transform of the aperture field given by [Ref. 5] 



P x - o , P y = J " /2 f V * E exp(jPrV) dV dz,' 

y -b/2 -a /Z 

c 

where r is the unit vector in the direction of propagation and r' is the vector pointing 
towards the aperture. The boresight far field field strength from one element is then 
found to be equal to 



2E ab 

C 

JtrX 

Assuming that the distance r is so far that all fields arriving at that point from 
all the X elements are in phase, we then have 



2XE a b 
IEJ = 



"cb 1 



7t rX 



Hence, the gain ratio between the two models is 



lEcbl 

lE'J 



2XE a 

C C 

nE 'u a 'u 
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Since E y = l/2E c and a u = Na c , assuming both apertures are fed with the same input 
power 




This means that the gain of the actual array is about 2.098 dB higher than that of the 
simplified model. 

2. COSINE TAPERED TE10 LINEAR ARRAY VS SINGLE TE10 SLOT 
ANTENNA 

a. Comparing Maximum Aperture E-Field 

The maximum aperture E-field of a single TE 10 slot antenna with length a' is 

given by, 



e . 2 _ 4 V,„ 

E ‘ ~~^b~ 



while that of a Cosine tapered TE 10 linear array is 



2 . 



4Z 0 P in 



a c b j^ E r 2 < n ) 



where E f (n) is the relative E-field of the n* element with maximum relative E-field 
equal to unity. Since Cosine tapering is used, the relative E-field may be expressed in 
terms of n as follow : 



E f (n) = sin 



nJt 



N+l 

and the total sum over N elements is 




N+l 

~ 2 ~ 



Therefore, the maximum aperture E-field is expressed as, 




8 Z 0 P in 

a c b(N+l) 



The ratio of the maximum aperture E-field for the same input power is 
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E c = V(2a' c ) 

E' c V[a c (N+l)] 

By equating a c = (N+ l)a c , we have 

E c = V2 x E' c . 

The actual array maximum aperture E-field is equal to that of the simplified model 
multiplied by *J2. 



b. Comparing Gain 

As derived before, the boresight far field field strength of the single TE 10 
aperture antenna with length a c ' is given by 



l E cbl 



2 E' a' b 

C C 

TtrX 



For a Cosine tapered TE 10 linear array with far field distance r long enough 
to guarantee in-phase E-field addition, we have 



lEcbl 



2a bE 

C C 

7trX 






2a c bE c sin[it/(N+ 1)] 
TtrX 1 - cos[7t/(N+ 1)] 
2abE 

= — ^ cot{7i/[2(N+l)]} . 



Hence, the gain ratio is, 



l E (±J _ a c E c 
lE'cbl 



, -C0t{7t/[2(N+1)]} . 

a c E c 



From above, w r e have E c = V 2 x E' c , for a' c = (N+ l)a c . 



|E rh l 



cb 1 _ 



V2 



lE'cbl 



N+l 



cot{«/[2(N+l)]} . 



Table 4 shows the gain ratio as N varies. 

From the table, it shows that for N > 8, the boresight gain difference is 
reduced to 1 dB or less. 
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TABLE 4 

GAIN RATIO BETWEEN COSINE TAPERED TE10 LINEAR ARRAY AND SINGLE TE10 

SLOT ANTENNA 


N 


Ratio 


Difference in dB 


2 


0.8165 


-1.7610 


5 


0.8797 


-1.1138 


8 


0.8912 


-1.0009 


10 


0.8942 


-0.9714 


20 


0.8986 


-0.9283 


50 


0.9000 


-0.9148 


100 


0.9002 


-0.9128 


200 


0.9003 


-0.9123 



3. COSINE TAPERED UNIFORMLY DISTRIBUTED LINEAR ARRAY VS 
SINGLE TE10 SLOT ANTENNA 

a. Comparing Maximum Aperture E-Field 

The maximum aperture E-field of a Cosine tapered uniformly distributed linear 

array with element length a ' can be expressed as follows: 




2Z n P- 

0 in 



a„b(N+l) 



and the maximum aperture E-field of a single TE 10 slot antenna is 

, 4Z n P. 

° " a- c b • 

Hence, the maximum aperture E-field ratio is, 



u _ 



Va'_ 



E' c V[2a u (N+l)] 
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By letting a c = a u (N+ 1), 



E u _ ' 



E - c V 2 

b. Comparing Gain 

The boresight far field field strength of a single element uniformly distributed 
aperture is, 

E u a u b 



rX 



as given before. With N Cosine tapered elements, the far field field strength is 



E a b il 

l E J - rr ! i E » 



rX 



E a b 

u u - cot (n/[2(N+l)]} , 



rX 

assuming r is large enough to guarantee in-phase E-field addition. 

The boresight far’ field field strength of the single TEj Q aperture is, 

2E' a' b 
|E' r J 



JtrX 



Hence, the gain ratio is, 

IE .1 kE a 

“ 7 “ = u cot (lt/(2(N+ 1)]} • 
|EJ 2E c a c 

With a' c = a u (N+ 1) and E u = (1/V2)E' C , we have 

it 

o /o 
cb 1 



\M 

|EVJ 2 V 2 (N + 1 ) 



cot{jt/[2(N+ 1)]} . 



Table 5 shows the boresight gain ratio as N varies. It shows that for large N, 
the difference in dB approaches -3.01 dB. 
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TABLE 5 

GAIN RATIO BETWEEN COSINE TAPERED UNIFORNLY DISTRIBUTED LINEAR 
ARRAY AND SINGLE TE10 SLOT ANTENNA 


N 


Ratio 


Difference in dB 


2 


0.6413 


-3.8591 


5 


0.6909 


-3.2120 


10 


0.7023 


-3.0696 


20 


0.7058 


-3.0265 


50 


0.7069 


-3.0130 


100 


0.7070 


-3.0110 


200 


0.7071 


-3.0105 
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APPENDIX D 

INPUT DATA FILE FOR GTD2P51 



CT (FWRP1W1F ) $ 

CE 

UN GEOMETRY UNIT IN M 

1 

US SOURCE UNIT IN WAVELENGTH 

0 

FR FREQUENCY IS 10.0 GHZ 

10.0 

RT ROTATE THE ANTENNA REFLECTOR 

0 , 0,0 

45.0. 90.270 

PG REFLECTOR PLATE - FRONT PLATE 

4,0 

1.35.0. 05,-2.55 
-1.35,0.05,-2.55 
-1.35,0.05,2.55 

1.35.0. 05.2.55 

PG REFLECTOR PLATE - BACK PLATE 

4,0 

1.35, -0.05,-2.55 

1.35, -0.05,2.55 
-1.35,-0.05,2.55 
-1.35,-0.05,-2.55 

PG REFLECTOR PLATE - TOP PLATE 

4,0 

1.35,0.05,2.55 

-1.35,0.05,2.55 

-1.35,-0.05,2.55 

1.35,-0.05,2.55 

PG REFLECTOR PLATE - BOTTOM PLATE 



4,0 

1.35,0.05,-2.55 

1.35, -0.05,-2.55 
-1.35,-0.05,-2.55 
-1.35,0.05,-2.55 

PG REFLECTOR PLATE - LEFT PLATE 

4,0 

1.35.0. 05,-2.55 

1.35.0. 05.2.55 

1.35, -0.05,2.55 

1.35, -0.05,-2.55 

PG REFLECTOR PLATE - RIGHT PLATE 



4,0 

-1.35,0.05,-2.55 

-1.35,-0.05,-2.55 

-1.35,-0.05,2.55 

-1.35,0.05,2.55 

RT ROTATE THE ARRAY ANTENNA 



0 , 0,0 

0,0,90,270 

CG PIVOT FOR REFLECTOR 



-1.46,0,0 
90,0,90,90 
0.14,0.14 
-0.11,90,0.11,90 
CG PIVOT FOR REFLECTOR 



1.46.0. 0 

90.0. 90.90 
0.14,0.14 
-0.11,90,0.11,90 

CG SIDE SUPPORTING CYLINDER 1 
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CG 



CG 



CG 



CG 



CG 



PG 



PG 



PG 



PG 



PG 



SA 



-1.84,0,-1.37 

16.314,0,106.314,0 

0.05,0.05 

-1.27,73.686,1.28,73.686 

SIDE SUPPORTING CYLINDER 2 



1.84,0,-1.37 

16.314,180,106.314,180 

0.05,0.05 

-1.27,73.686,1.28,73.686 

SIDE SUPPORTING CYLINDER 3 
-1.48,0.595,-1.35 
22,-90,112,-90 
0.05,0.05 
-1.34,68,1.35,100 

SIDE SUPPORTING CYLINDER 4 
-1.48,-0.595,-1.35 
22,90,112,90 
0.05,0.05 
-1.34,68,1.35,100 

SIDE SUPPORTING CYLINDER 5 



1.48,0.595,-1.35 

22,-90,112,-90 

0.05,0.05 

-1.34,68,1.35,100 

SIDE SUPPORTING CYLINDER 6 
1.48,-0.595,-1.35 
22,90,112,90 
0.05,0.05 
-1.34,68,1.35,100 

ARRAY ANTENNA GROUND PLATE 

4,0 

2. 25, 1.5, -2. 6 
-2. 25, 1.5, -2. 6 
-2. 25, -1.5, -2. 6 

2. 25, -1.5, -2. 6 

ARRAY ANTENNA SIDE PLATE 1 

4,0 

2. 25, -1.5, -2. 6 
-2. 25, -1.5, -2. 6 
-2. 25, -1.5, -3.1 

2. 25, -1.5, -3.1 

ARRAY ANTENNA SIDE PLATE 2 

4,0 

2. 25, 1.5, -2. 6 

2. 25, -1.5, -2. 6 

2. 25, -1.5, -3.1 
2.25,1.5 -3.1 

ARRAY ANTENNA SIDE PLATE 3 

4,0 

-2. 25, 1.5, -2. 6 

2. 25. 1.5, -2. 6 

2. 25. 1.5, -3.1 
-2.25 1.5 -3.1 

ARRAY' ANTENNA SIDE PLATE 4 



4,0 

-2. 25, -1.5, -2. 6 
-2. 25, 1.5, -2. 6 
-2. 25, 1.5, -3.1 
“2 25 "1 5 “3 1 

ARRAY ELEMENTS CONFIGURATION 
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-0.0000003, 

-0.0000003, 

-0.0000003, 

-0.0000003, 

-0.0000003, 

-0.0000003, 

-0.0000003, 

-0.0000003, 



-0.9893153, 

-0.9743257, 

-0.9593360, 

-0.9443464, 

-0.9293568, 

-0.9143671, 

-0.8993775, 

-0.8843879, 



-2.5999 

-2.5999 

-2.5999 

-2.5999 

-2.5999 

-2.5999 

-2.5999 

-2.5999 
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-0.0000003, 


-0.8693983, 


-2.5999 


-0.0000003, 


-0.8544087, 


-2.5999 


-0.0000003, 


-0.8394190, 


-2.5999 


-0.0000003, 


-0.8244294, 


-2.5999 


-0.0000003, 


-0.8094398, 


-2.5999 


-0.0000002, 


-0.7944502, 


-2.5999 


-0.0000002, 


-0.7794605, 


-2.5999 


-0.0000002, 


-0.7644709, 


-2.5999 


-0.0000002, 


-0.7494813, 


-2.5999 


-0.0000002, 


-0.7344916, 


-2.5999 


-0.0000002, 


-0.7195020, 


-2.5999 


-0.0000002, 


-0.7045124, 


-2.5999 


-0.0000002, 


-0.6895228, 


-2.5999 


-0.0000002, 


-0.6745331, 


-2.5999 


-0.0000002, 


-0.6595435, 


-2.5999 


-0.0000002, 


-0.6445539, 


-2.5999 


-0.0000002, 


-0.6295643, 


-2.5999 


-0.0000002, 


-0.6145747, 


-2.5999 


-0.0000002, 


-0.5995850, 


-2.5999 


-0.0000002, 


-0.5845954, 


-2.5999 


-0.0000002, 


-0.5696058, 


-2.5999 


-0.0000002, 


-0.5546162, 


-2.5999 


-0.0000002, 


.-0.5396265, 


-2.5999 


-0.0000002, 


-0.5246369, 


-2.5999 


-0.0000002, 


-0.5096473, 


-2.5999 


-0.0000002, 


-0.4946576, 


-2.5999 


-0.0000002, 


-0.4796680, 


-2.5999 


-0.0000001, 


-0.4646784, 


-2.5999 


-0.0000001, 


-0.4496888, 


-2.5999 


-0.0000001, 


-0.4346991, 


-2.5999 


-0.0000001, 


-0.4197095, 


-2.5999 


-0.0000001, 


-0.4047199, 


-2.5999 


-0.0000001, 


-0.3897303, 


-2.5999 


-0.0000001, 


-0.3747407, 


-2.5999 


-0.0000001, 


-0.3597510, 


-2.5999 


-0.0000001, 


-0.3447614, 


-2.5999 


-0.0000001, 


-0.3297718, 


-2.5999 


-0.0000001, 


-0.3147821, 


-2.5999 


-0.0000001, 


-0.2997925, 


-2.5999 


-0.0000001, 


-0.2848029, 


-2.5999 


-0.0000001, 


-0.2698132, 


-2.5999 


-0.0000001, 


-0.2548236, 


-2.5999 


-0.0000001, 


-0.2398340, 


-2.5999 


-0.0000001, 


-0.2248444, 


-2.5999 


-0.0000001, 


-0.2098548, 


-2.5999 


-0.0000001, 


-0.1948651, 


-2.5999 


-0.0000001, 


-0.1798755, 


-2.5999 


-0.0000001, 


-0.1648859, 


-2.5999 


0 . 0000000 , 


-0.1498963, 


-2.5999 


0 . 0000000 , 


-0.1349066, 


-2.5999 


0 . 0000000 , 


-0.1199170, 


-2.5999 


0 . 0000000 , 


-0.1049274, 


-2.5999 


0 . 0000000 , 


-0.0899377, 


-2.5999 


0 . 0000000 , 


-0.0749481, 


-2.5999 


0 . 0000000 , 


-0.0599585, 


-2.5999 


0 . 0000000 , 


-0.0449688, 


-2.5999 


0 . 0000000 , 


-0.0299792, 


-2.5999 


0 . 0000000 , 


-0.0149896, 


-2.5999 


0 . 0000000 , 


0 . 0000000 , 


-2.5999 


0 . 0000000 , 


0.0149889, 


-2.5999 


0 . 0000000 , 


0.0299788, 


-2.5999 


0 . 0000000 , 


0.0449686, 


-2.5999 


0 . 0000000 , 


0.0599585, 


-2.5999 


0 . 0000000 , 


0.0749474, 


-2.5999 


0 . 0000000 , 


0.0899372, 


-2.5999 


0 . 0000000 , 


0.1049271, 


-2.5999 


0 . 0000000 , 


0.1199169, 


-2.5999 


0 . 0000000 , 


0.1349058, 


-2.5999 


0 . 0000000 , 


0.1498957, 


-2.5999 
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0.0000001, 


0.1648855, 


-2.5999 


0.0000001, 


0.1798754, 


-2.5999 


0.0000001, 


0.1948643, 


-2.5999 


0.0000001, 


0.2098541, 


-2.5999 


0.0000001, 


0.2248440, 


-2.5999 


0.0000001, 


0.2398338, 


-2.5999 


0.0000001, 


0.2548227, 


-2.5999 


0.0000001, 


0.2698126, 


-2.5999 


0.0000001, 


0.2848024, 


-2.5999 


0.0000001, 


0.2997923, 


-2.5999 


0.0000001, 


0.3147821, 


-2.5999 


0.0000001, 


0.3297710, 


-2.5999 


0.0000001, 


0.3447609, 


-2.5999 


0.0000001, 


0.3597507, 


-2.5999 


0.0000001, 


0.3747406, 


-2.5999 


0.0000001, 


0.3897295, 


-2.5999 


0.0000001, 


0.4047194, 


-2.5999 


0.0000001, 


0.4197092, 


-2.5999 


0.0000001, 


0.4346991, 


-2.5999 


0.0000001, 


0.4496880, 


-2.5999 


0.0000001, 


0.4646778, 


-2.5999 


0.0000002, 


0.4796677, 


-2.5999 


0.0000002, 


0.4946575, 


-2.5999 


0.0000002, 


0.5096464, 


-2.5999 


0.0000002, 


0.5246363, 


-2.5999 


0.0000002, 


0.5396263, 


-2.5999 


0.0000002, 


0.5546160, 


-2.5999 


0.0000002, 


0.5696049, 


-2.5999 


0.0000002, 


0.5845947, 


-2.5999 


0.0000002, 


0.5995846, 


-2.5999 


0.0000002, 


0.6145744, 


-2.5999 


0.0000002, 


0.6295643, 


-2.5999 


0.0000002, 


0.6445532, 


-2.5999 


0.0000002, 


0.6595430, 


-2.5999 


0.0000002, 


J 0.6745329, 


-2.5999 


0.0000002, 


0.6895227, 


-2.5999 


0.0000002, 


0.7045116, 


-2.5999 


0.0000002, 


0.7195015, 


-2.5999 


0.0000002, 


0.7344913, 


-2.5999 


0.0000002, 


0.7494812, 


-2.5999 


0.0000002, 


0.7644701, 


-2.5999 


0.0000002, 


0.7794600, 


-2.5999 


0.0000002, 


0.7944498, 


-2.5999 


0.0000003, 


0.8094397, 


-2.5999 


0.0000003, 


0.8244286, 


-2.5999 


0.0000003, 


0.8394184, 


-2.5999 


0.0000003, 


0.8544083, 


-2.5999 


0.0000003, 


0.8693981, 


-2.5999 


0.0000003, 


0.8843870, 


-2.5999 


0.0000003, 


0.8993769, 


-2.5999 


0.0000003, 


0.9143667, 


-2.5999 


0.0000003, 


0.9293566, 


-2.5999 


0.0000003, 


0.9443464, 


-2.5999 


0.0000003, 


0.9593353, 


-2.5999 


0.0000003, 


0.9743252, 


-2.5999 


0.0000003, 


0.9893150, 


-2.5999 


90.0, 0.0, 


90.0, 90.0 




1, 83.39101, 


0.50000 




0 . 124507E+05 , 


0.0 




0 . 248939E+05 , 


0.0 




0 . 373239E+05 , 


0.0 




0 . 497329E+05 , 


0.0 




0 . 621151E+05 , 


0.0 




0 . 744631E+05 , 


0.0 




0 . 867697E+05 , 


0.0 




0 . 990291E+05 , 


0.0 




0 . 111234E+06 , 


0.0 




0 . 123378E+06 , 


0.0 




0 . 135453E+06 , 


0.0 
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0 . 147455E+06 , 0.0 
0 . 159375E+06 , 0.0 
0 . 171208E+06 , 0.0 
0 . 182946E+06 , 0.0 
0 . 194584E+06 , 0.0 
0 . 206116E+06 , 0.0 
0 . 217533E+06 , 0.0 
0 . 228832E+06 , 0.0 
0.240004E+06, 0.0 
0 . 251045E+06 , 0.0 
0 . 261948E+06 , 0.0 
0 . 272706E+06 , 0.0 
0 . 283315E+06 , 0.0 
0 . 293768E+06 , 0.0 
0 . 304059E+06 , 0.0 
0 . 314184E+06 , 0.0 
0.324135E+06, 0.0 
0 . 333909E+06 , 0.0 
0 . 343499E+06 , 0.0 
0 . 352900E+06 , 0.0 
0 . 362108E+06 , 0.0 
0 . 371116E+06 , 0.0 
0 . 379920E+06 , 0.0 
0 . 388516E+06 , 0.0 
0 . 396898E+06 , 0.0 
0.405062E+06, . 0.0 
0 . 413003E+06 , 0.0 
0 .4207 17E+06 , 0.0 
0 . 428200E+06 , 0.0 
0 . 435448E+06 , 0.0 
0.442456E+06, 0.0 
0 . 449221E+06 , 0.0 
0.455739E+06, 0.0 
0.462007E+06, 0.0 
0 . 468021E+06 , 0.0 
0.473777E+06, 0.0 
0.479273E+06, 0.0 
0 . 484506E+06 , 0.0 
0.489472E+06, 0.0 
0.494170E+06, 0.0 
0.498595E+06, 0.0 
0 . 502747E+06 , 0.0 
0 . 506622E+06 , 0.0 
0 . 510219E+06 , 0.0 
0 . 513536E+06 , 0.0 
0 . 516570E+06 , 0.0 
0 . 519320E+06 , 0.0 
0 . 521785E+06 , 0.0 
0 . 523963E+06 , 0.0 
0 . 525854E+06 , 0.0 
0 . 527455E+06 , 0.0 
0 . 528766E+06 , 0.0 
0 . 529786E+06 , 0.0 
0 . 530516E+06 , 0.0 
0 . 530954E+06 , 0.0 
0 . 531099E+06 , 0.0 
0 . 530954E+06 , 0.0 
0 . 530516E+06 , 0.0 
0 . 529786E+06 , 0.0 
0 . 528766E+06 , 0.0 
0 . 527455E+06 , 0.0 
0 . 525854E+06 , 0.0 
0 . 523963E+06 , 0.0 
0 . 521785E+06 , 0.0 
0 . 519321E+06 , 0.0 
0 . 516570E+06 , 0.0 
0 . 513536E+06 , 0.0 
0 . 510219E+06 , 0.0 
0 . 506623E+06 , 0.0 
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0.502747E+06, 0.0 

0.498595E+06, 0.0 

0.494170E+06, 0.0 

0.489473E+06, 0.0 

0 .484506E+06 , 0.0 

0.479273E+06, 0.0 

0.473777E+06, 0.0 

0.468021E+06, 0.0 

0.462007E+06, 0.0 

0.455739E+06, 0.0 

0 . 449221E+06 , 0.0 

0 . 442456E+06 , 0.0 

0.435448E+06, 0.0 

0 .428200E+06 , 0.0 

0.420717E+06, 0.0 

0.413003E+06, 0.0 

0 . 405062E+06 , 0.0 

0 . 396898E+06 , 0.0 

0.388516E+06, 0.0 

0 . 379920E+06 , 0.0 

0 . 371117E+06 , 0.0 

0 . 362108E+06 , 0.0 

0 . 352901E+06 , 0.0 

0 . 343499E+06 , 0.0 

0 . 333910E+06 , 0.0 

0 . 324136E+06 , 0.0 

0 . 314184E+06 , 0.0 

0 . 304059E+06 , 0.0 

0 . 293768E+06 , 0.0 

0 . 283315E+06 , 0.0 

0 . 272706E+06 , 0.0 

0 . 261948E+06 , 0.0 

0.251045E+06, 0.0 

0 . 240005E+06 , .0.0 

0 . 228832E+06 , 0.0 

0.217534E+06, 0.0 

0 . 206116E+06 , 0.0 

0 . 194585E+06 , 0.0 

0 . 182947E+06 , 0.0 

0 . 171208E+06 , 0.0 

0 . 15937 5E+06 , 0.0 

0 . 147455E+06 , 0.0 

0 . 135454E+06 , 0.0 

0 . 123378E+06 , 0.0 

0 . 111234E+06 , 0.0 

0.990297E+05, 0.0 

0 . 867702E+05 , 0.0 

0 .744631E+05 , 0.0 

0 . 621151E+05 , 0.0 

0 . 497329E+05 , 0.0 

0.373244E+05, 0.0 

0 . 248944E+05 , 0.0 

0 . 124507E+05 , 0.0 

RT RESET TO ORIGINAL GLOBAL COORDINATE 

0,0,0 
0,0,90,0 

PG VEHICLE BODY SIDE PLATE 1 

4,0 

1.75, -1.5, -3. 2 
-2. 25, -1.5, -3. 2 
-2. 25, -1.5, -4. 25 

1.75, -1.5, -4. 25 

PG VEHICLE BODY SIDE PLATE 2 

4,0 

3. 15, -1.4, -3. 2 

1.75, -1.4, -3. 2 

1.75, -1.4, -4.1 
3.15 -1.4 -4.1 

PG ' VEHICLE BODY SIDE PLATE 3 



149 



4,0 

3. 3, -1.1, -3. 2 

3. 15, -1.4, -3. 2 

3. 15, -1.4, -4.1 

3. 3, -1.1, -4.1 

PG VEHICLE BODY SIDE PLATE 4 

4,0 

3. 15. 1.4, -3. 2 

3. 3. 1.1, -3. 2 

3. 3. 1.1, -4.1 

3. 15. 1.4, -4.1 

PG VEHICLE BODY SIDE PLATE 5 

4,0 

1.75.1. 4, -3. 2 

3. 15. 1.4, -3. 2 

3. 15. 1.4, -4.1 

1.75. 1.4, -4.1 

PG VEHICLE BODY SIDE PLATE 6 

4,0 

-2. 25, 1.5, -3. 2 

1.75. 1.5, -3. 2 

1.75.1.5, -4.25 
-2.25,1.5,-4.25 

PG VEHICLE BODY SIDE PLATE 7 

4,0 

1.75, -1.4, -3. 2 

1.75, -1.5, -3. 2 

1.75, -1.5, -4.1 

1.75, -1.4, -4.1 

PG VEHICLE BODY SIDE PLATE 8 

4,0 

1.75. 1.5, -3. 2 

1.75. 1.4, -3. 2 

1.75. 1.4, -4.1 

1.75. 1.5, -4.1 

PG VEHICLE BODY BACK PLATE 

4.0 

-2. 25, 1.5, -3. 2 
-2. 25, -1.5, -3. 2 
-2.25,-115,-4.25 
-2.25,1.5,-4.25 

PG VEHICLE BODY TOP PLATE 

12.0 

2. 8. 1.1, -3. 2 

3. 3. 1.1, -3. 2 

3. 15. 1.4, -3. 2 

1.75. 1.4, -3. 2 

1.75. 1.5, -3. 2 
-2. 25, 1.5, -3. 2 
-2. 25, -1.5, -3. 2 

1.75, -1.5, -3. 2 

1.75, -1.4, -3. 2 

3. 15, -1.4, -3. 2 

3. 3, -1.1, -3. 2 
28-1 1 -32 

PG ' ' VEHICLE* BODY BOTTOM PLATE 1 
4,0 

4. 45, -1.1, -4.1 

3. 3, -1.1, -4.1 

3. 15, -1.4, -4.1 
4.45 -1.4 -4.1 

PG ’ VEHICLE BODY BOTTOM PLATE 2 
4,0 

4. 45. 1.4, -4.1 

3. 15. 1.4, -4.1 

3. 3. 1.1, -4.1 

4. 45. 1.1, -4.1 

PG VEHICLE HEAD TOP PLATE 

4,0 



150 



3. 9.1. 1, -2. 6 

2 . 8 . 1 . 1 , - 2. 6 
2. 8, -1.1, -2. 6 

3. 9, -1.1, -2. 6 

PG VEHICLE HEAD BACK PLATE 

4,0 

2. 8, -1.1, -2. 6 
2. 8, 1.1, -2. 6 

2. 8. 1.1, -3. 2 
2 8 “1 1 -3 2 

PG * ' VEHICLE' HEAD FRONT UPPER PLATE 
4,0 

3. 9. 1.1, -2. 6 

3. 9, -1.1, -2. 6 

4. 45, -1.1, -3. 5 

4. 45. 1.1, -3. 5 

PG VEHICLE HEAD FRONT LOWER PLATE 

4.0 

4. 45. 1.1, -3. 5 

4. 45, -1.1, -3. 5 

4. 45, -1.1, -4.1 

4. 45. 1.1, -4.1 

PG VEHICLE HEAD SIDE PLATE 1 

7.0 

3. 9, -1.1, -2. 6 
2. 8, -1.1, -2. 6 
2. 8, -1.1, -3. 2 

3. 3, -1.1, -3. 2 

3. 3, -1.1, -4.1 

4. 45, -1.1, -4.1 

4. 45, -1.1, -3. 5 

PG VEHICLE HEAD SIDE PLATE 2 

7,0 

2 . 8 . 1 . 1 , - 2. 6 

3. 9. 1.1, -2. 6 

4. 45. 1.1, -3. 5 

4. 45. 1.1, -4.1 

3. 3. 1.1, -4.1 

3. 3. 1.1, -3. 2 
2811 ~3 2 

LP * 'data output enabled 

T 

RD FAR FIELD RANGE 

10000 

PF FAR FIELD RADIATION PATTERN PLOT 1 

0.0,0.0,90.0,0.0 
F , 0 . 0 

45.0,0.25,361 
XQ EXECUTE 

PF FAR FIELD RADIATION PATTERN PLOT 2 

0.0,0.0,90.0,0.0 
T , 90 . 0 

-45.0,0.25,361 
XQ EXECUTE 

EN END OF INPUT DATA 
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